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I.  SUMMARY  OF  RESEARCH  ACCOMPLISHMENTS 

This  report  summarizes  the  research  activites  carried  out  at  the  Insti¬ 
tute  for  Physical  Science  and  Technology,  University  of  Maryland,  under  the 
auspices  of  the  Office  of  Naval  Research  Grant  N00014-82-K-0208  from  March  1, 
1982  to  March  31,  1984.  We  are  pleased  to  report  that  significant  progress 
has  continued  in  the  study  of  the  acceleration  of  electrons  and  ions  by  wave- 
particle  interactions.  Several  theoretical  and  numerical  investigations  of 
these  acceleration  processes  with  applications  to  space  physics  and  laboratory 
plasmas  have  been  conducted.  Copies  of  the  six  completed  papers  are  attached. 
The  ONR  grant  was  acknowledged  appropriately.  Several  other  manuscripts  are 
presently  in  preparation.  Also  included  for  your  information  are  eleven 
abstracts  of  talks  presented  at  scientific  meetings.  Dr.  Gaffey  chaired  two 
sessions  at  the  Spring  1984  Meeting  of  the  American  Geophysical  Union. 

1 .  Overview 

The  research  program  is  devoted  to  analyzing  the  processes  respon¬ 
sible  for  the  acceleration  of  ions  and  electrons  in  plasmas,  to  studying 
instabilities  resulting  from  energetic  ion  beams  streaming  across  a  magnetic 
field,  or  other  non-equilibrium  features,  such  as  temperature  anisotropy  or  a 
loss-cone  distribution  function,  and  to  calculating  the  spontaneous  emission 
of  radiation  at  the  fundamental  and  harmonics  of  the  cyclotron  frequency. 

The  program  is  divided  into  six  closely  related  projects.  The 
first  and  second  parts  deal  with  electromagnetic  and  electrostatic  instabili¬ 
ties  resulting  from  an  energetic  ion  beam  streaming  across  the  magnetic  field 
and  electron  t_emperature  anisotropy  in  the  earth^s  bow  shock  and  in  magnetic 
mirrors.  The  third  part  treats  the  acceleration  of  electrons  by  plasma  turbu¬ 
lence  and  an  electric  field  parallel  to  the  magnetic  field.  The  fourth  part 
considers  the  electron-cyclotron  maser  instability  in  a  high  temperature 
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plasma  with  a  loss-cone  distribution  function*  In  the  fifth  part,  the  spec¬ 
tral  emissivity  of  spontaneous  synchrotron  radiation  at  the  fundamental  and 
harmonics  of  the  cyclotron  frequency  is  calculated  using  the  complete  relativ¬ 
istic  resonance  condition*  In  the  sixth  part,  the  stochastic  acceleration  of 
ions  by  a  large  amplitude  ion  cyclotron  wave  is  discussed. 

The  problems  discussed  above  are  all  interrelated.  For  example, 
the  kinetic  cross-field  streaming  instability  considered  in  sections  2  and  3 
can  generate  the  plasma  turbulence  which  results  in  the  electron  acceleration 
considered  in  section  4*  Moreover,  the  acceleration  of  the  electrons  along 
the  magnetic  field  can  result  in  a  loss-cone  distribution  function,  which  can 
be  unstable  to  the  electron  cyclotron  maser  instability  discussed  in  sec¬ 
tion  5.  The  expression  for  the  spontaneous  cyclotron  emissivity  obtained  in 
section  6  has  been  used  to  calculate  the  radiation  spectrum  from  energetic 
trapped  electrons  with  a  loss-cone  distribution  function.  Finally,  the 
electron  current  along  the  magnetic  field  discussed  in  section  4  can  excite 
the  ion  cyclotron  waves  which  produce  the  stochastic  ion  acceleration  consid¬ 
ered  in  section  7. 

We  now  briefly  describe  the  key  results  and  major  accomplishments 
achieved  thus  far  in  our  research  program.  A  detailed  report  on  each  problem 
is  given  in  Chapter  III. 

2.  Effect  of  Electron  Thermal  Anisotropy  on  the  Kinetic  Cross-Field 
Streaming  Instability 

We  have  investigated  both  numerically  and  analytically  the  kinetic 
cross— field  streaming  instability,  which  is  excited  by  a  relative  electron— ion 
cross-field  drift  and  enhanced  by  an  electron  temperature  anisotropy  T  > 
Te i| "  Tiis  investigation  is  motivated  by  the  research  on  collisionless  shock 
waves  and  applied  to  the  earth's  bow  shocki’^  The  unstable  waves  have  been- 
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shown  to  be  obliquely  propagating  whistlers,  and  the  instability  is  not  sup¬ 
pressed  by  electromagnetic  effects  when  the  drift  velocity  VQ  exceeds  the 
Alfven  velocity  v^  and  the  plasma  parameter  3  is  of  order  unity*  The  effect 
of  increasing  Vq/va  is  to  shift  the  direction  of  propagation  more  nearly 
parallel  to  the  magnetic  field.  The  effect  of  increasing  the  electron  temper¬ 
ature  ratio  TeJ/Te|  is  to  increase  the  growth  rate  and  to  shift  the  direction 
of  propagation  more  nearly  perpendicular  to  the  magnetic  field.  The  mode  has 
mixed  electrostatic  and  electromagnetic  polarization.  However,  for  3-1  and 
for  nearly  perpendicular  propagation,  the  mode  is  predominantly  electrostatic; 
whereas,  for  more  oblique  propagation,  the  electromagnetic  and  electrostatic 
components  are  nearly  equal.  Finally,  for  small  3,  the  mode  becomes  the  well- 
known  modified  two— stream  instability? 

In  summary,  we  have  studied  the  plasma  instability  resulting  from  a 
cross— field  drift  and  an  electron  temperature  anisotropy,  T  /T  >  \  # 

find  that  both  effects  significantly  affect  the  instability.  In  particular, 
the  present  study  has  led  to  the  following  major  conclusions: 

(1)  A  relative  electron-ion  cross— field  drift  can  cause  insta¬ 
bility,  even  if  VQ/vA  >  1. 

(2)  An  electron  temperature  anisotropy  T  /Tej  >  1  significantly 
enhances  the  peak  growth  rate. 

(3)  The  unstable  waves  have  both  electrostatic  and  electromag¬ 
netic  components  in  general.  For  3  -  1  and  nearly  parallel  propagation,  the 

electrostatic  and  electromagnetic  components  are  comparable  and  the  polariza- 
tion  is  mixed. 

(4)  The  instability  is  highly  kinetic  for  high  0  -  1.  From 

numerical  solutions  of  the  dispersion  equation,  we  find  that 
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over  the  entire  range  of  angles  of  propagation.  Moreover,  electromagnetic 
effects,  high-order  electron  cyclotron  harmonics,  and  ion  terms  all  affect  the 
dispersion  equation  significantly  and,  therefore,  it  is  difficult  to  treat  the 
dispersion  equation  accurately  by  analytic  techniques. 

3.  Instabilities  Excited  by  an  Energetic  Ion  Beam  and  Electron  Temper¬ 
ature  Anisotropy  in  Tandem  Mirrors 

Tandem  Mirrors  are  magnetic  confinement  devices  in  which  the  ions 
in  the  central  (solenoidal)  cell  are  prevented  from  leaking  out  the  end  by 
means  of  an  electrostatic  potential  maintained  by  a  denser  plasma  in  mirror- 


end  cells.  The  basic  concept  has  been  successfully  verified  in  the  Tandem 
Mirror  Experment  (TMX).  In  order  to  generate  the  barrier,  the  density  in  this 
region  is  reduced  by  pumping  out  some  of  the  ions.  This  is  achieved  by  injec¬ 
ting  neutral  beams  at  oblique  angles  to  create  "sloshing  ion  distributions," 
which  have  a  density  minimum  at  the  center  of  the  barrier.  Furthermore,  a  hot 
electron  population  is  also  created  in  the  barrier  by  means  of  ECRH  to  reduce 
the  power  transfer  between  the  plug  and  central-cell  electrons.  The  thermal 
barrier  concept  is  presently  being  tested  in  the  upgraded  TMX-U  experiment^ 

While  the  inclusion  of  such  complexity  into  the  end  cell  of  the 
Tandem  Mirror  may  improve  the  overall  confinement  significantly,  it  is  implic¬ 
itly  assumed  that  such  configurations  can  be  stably  maintained.  However,  the 
introduction  of  (neutral)  ion  beams  and  hot  anisotropic  electrons  in  the  bar¬ 
rier  region  suggests  a  variety  of  possible  microinstabilities.  Injection  of 
neutral  beams  to  form  the  sloshing  ion  distribution  allows  the  possibility  of 
various  streaming  instabilities8’ 7  The  formation  of  hot  anisotropic  electrons 
may  also  give  rise  to  unstable  modes,  such  as  the  whistler8-11  and  the  elec¬ 
tron  cyclotron  maser  instabilities}2-21 
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In  a  recent  paper,  we  have  investigated  such  a  combination  of  ions 
drifting  across  a  magnetic  field  in  which  the  background  electron  temperature 
was  anisotropic  in  the  environment  of  the  earth's  bow  shock?^  In  that  situa¬ 
tion,  the  ion  streaming  results  from  ion  reflection  at  supercritical  shocks 
and  the  electron  anisotropy  arises  from  adiabatic  compression  at  the  shock 

surface.  We  have  found  that  the  well-known  (low-beta)  modified  two-stream 
23  l 

instability  is  modified  by  finite-beta  effects1  and  further  enhanced  by  the 

?? 

temperature  anisotropy. 

Using  the  same  basic  model  (local  linear  theory,  homogeneous 
plasma,  unmagnetized  ions,  etc.),  we  extend  the  theory  of  cross-field  ion 
streaming  instabilities  to  the  parameter  regime  of  the  end  cells  of  the 
upgraded  Tandem  Mirror  Experiment  with  neutral  beam  injection.  The  cross¬ 
field  ion  drift  as  well  as  the  electron  thermal  anisotropy  T  >  T  provides 
the  free— energy  which  drives  various  instabilities.  Three  instabilities,  a 
nearly  perpendicular  propagating  modified  two-stream  instability,  an  obliquely 
propagating  ion— ion  streaming  instability,  and  an  obliquely  propagating  elec¬ 
tromagnetic  lower-hybrid  instability,  have  been  identified.  The  first  two 
waves  are  electrostatic  and  have  the  largest  growth  rate.  For  the  actual 
operation  conditions  of  the  Tandem  Mirror  Experiment,  it  is  found  that  the 
ion-ion  streaming  instability  has  the  largest  growth  rate.  When  more  ener¬ 
getic  neutral  beams  become  available,  the  modified  two-stream  instability  may 
play  the  dominant  role  in  the  stability  of  these  devices. 

4.  Acceleration  of  Electrons  by  Plasma  Turbulence  and  a  Quasistatic 
Parallel  Electric  Field 

Recently,  some  laboratory  experiments  have  demonstrated  that  the 
combined  effect  of  a  parallel  electric  field  and  plasma  turbulence  is  effec¬ 
tive  in  producing  energetic  electrons?4  This  phenomenon  has  been  observed  in 
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a  tokamak,  and  we  expect  that  the  physical  mechanism  should  be  operative  in 
the  auroral  region  also.  It  is  based  on  this  notion  that  we  were  motivated  to 
proceed  with  the  present  study. 

The  acceleration  of  electrons  along  the  magnetic  field  by  large 
amplitude  plasma  turbulence  and  a  quasistatic  parallel  electric  field  has  been 
investigated.  The  plasma  amplifies  perturbations  excited  by  a  cross-field 
current  due  to  a  relative  electron-ion  drift.  The  fluctuating  electric  field 
is  primarily  electrostatic  and  propagates  nearly  perpendicularly  to  the  mag¬ 
netic  field.  The  instability  is  assumed  to  saturate  rapidly.  The  time  evolu¬ 
tion  of  the  velocity  and  energy  of  an  individual  test  electron  is  calculated 
from  the  moment  equations.  A  Greenes  function  is  constructed  for  the  quasi— 
linear  diffusion  equation,  with  a  power  law  spectral  energy  density?^  The 
time  evolution  of  an  initial  Maxwellian  distribution  is  given;  and  the  average 
electron  velocity,  energy,  and  velocity  spread  are  obtained  for  various  cases 
by  taking  moments  of  the  distribution  function.  We  have  found  that  the  com¬ 
bined  effect  of  a  parallel  electric  field  and  electrostatic  plasma  turbulence 
is  especially  effective  in  producing  energetic  electrons.  In  particular,  the 
present  study  has  led  to  the  following  major  conclusions: 

(1)  A  larger  electric  field  produces  greater  acceleration;  how¬ 
ever,  acceleration  by  an  electric  field  alone  does  not  increase  the  velocity 
spread  of  the  electron  distribution  function. 

(2)  A.  larger  spectral  index,  i.e.,  a  narrower  spectrum  of  waves, 
produces  greater  turbulent  acceleration.  Turbulent  acceleration  increases 

both  the  velocity  spread  and  the  average  energy  of  the  electron  distribution 
function. 

(3)  The  turbulent  acceleration  mechanism  is  more  effective  at 
higher  altitudes,  whereas  the  electric  field  acceleration  is  independent  of 
altitude . 
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(4)  For  some  values  of  the  parameters  which  are  typical  of  the 
auroral  zones  the  turbulent  acceleration  and  the  electric  field  acceleration 
are  comparable. 

(5)  When  both  acceleration  mechanisms  are  operable  the  accelera¬ 
tion  is  greater  than  when  only  one  acceleration  mechanism  is  operable. 

Although  the  results  have  been  applied  to  the  auroral  regions,  they 
also  have  important  implications  for  laboratory  experiments  as  well,  for 
example,  radio  frequency  current  drive  in  tokamaks. 

5 .  Electron-Cyclotron  Maser  Instability  in  Hot  Electron  Mirror 
Machines 

It  has  now  been  firmly  established  that  energetic  electrons  with  a 
loss— cone  distribution  function  are  unstable  to  the  electron— cyclotron  maser 
instability.  The  basic  physical  mechanism  of  the  instability  is  a  resonant 
wave— particle  interaction  via  a  Doppler— shifted  cyclotron  resonance  pro— 
cess»  The  relativistic  resonance  condition  has  been  shown  to  represent  an 
ellipse  in  velocity-space,  which  is  qualitatively  different  from  the  usual 
nonrelativistic  resonance  condition,  which  represents  a  line  in  velocity- 
space.  If  the  portion  of  the  distribution  function  covered  by  the  resonant 
ellipse  has  a  positive  derivative,  then  the  electron-cyclotron  maser  insta¬ 
bility  will  develop.  This  is  a  sufficient  condition  for  instability.  The 
free-energy  which  drives  the  instability  is  the  velocity— space  anisotropy  of 
the  distribution  function.  The  result  of  the  instability  is  the  transfer  of 
energy  from  the  energetic  electrons  to  the  electromagnetic  waves  via  a  rela¬ 
tivistic  Doppler— shif ted  cyclotron  resonance  process. 

This  instability  has  received  considerable  attention  in  the  past 
few  years  by  the  space  physics  community  as  a  possible  source  mechanism  of  the 
auroral  kilometric  radiation!^  ^  More  recently,  it  has  been  suggested  that 
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this  instability  may  also  occur  in  some  laboratory  experiments,  such  as  the 
Tandem  Mirror  Machine  and  the  EBT  device.  In  a  recent  paper  by  Lau  and  Chu}9 
a  special  case  of  the  cyclotron  maser  instability  associated  with  a  loss-cone 
distribution  function  is  discussed  for  hot  electron  mirror  experiments.  Their 
analysis  is  restricted  to  the  case  that  the  radiation  propagates  parallel  to 
the  applied  magnetic  field  Bq  (k^  =  0)  and,  moreover,  emphasis  is  placed  on 
the  limit  kj(  *  0.  Actually,  before  Ref.  19  was  published,  the  theory  of  such 
an  instability  had  been  extensively  developed  by  numerous  authors , 20 , 26 
(mostly  in  the  space  physics  community),  although  the  research  has  been  moti¬ 
vated  by  the  investigation  of  the  auroral  kilometric  radiation. 

In  view  of  the  fact  that  strong  bursts  of  cyclotron  radiation  have 
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been  observed  in  hot  electron  mirror  machines,  it  is  important  and  inter¬ 

esting  to  study  the  cyclotron  maser  instability  for  such  experiments  more 
fully. 

We  have  generalized  the  work  of  Lau  and  Chu  to  treat  the  case  of  a 
plasma  consisting  of  a  background  of  low-energy  electrons  as  well  as  a  popula¬ 
tion  of  suprathermal  electrons,  and  our  analysis  is  carried  out  for  arbitrary 
direction  of  propagation.  The  major  conclusions  of  our  preliminary  study  are: 

(1)  The  cyclotron  maser  instability  can  occur  at  all  angles  of 
propagation  over  a  wide  range  of  parameters. 

(2)  In  general,  the  instability  at  large  0  (i.e.,  0  ~  90°)  is 
more  serious  than  that  at  0  =  0°. 

(3)  The  peak  of  the  growth  rate  maximized  over  k  increases  with 
the  temperature  T  and,  for  large  T,  occurs  at  a  finite  k  rather  than  k  =  0. 

(4)  The  range  of  unstable  values  of  w  /fi  broadens  as  T 

e  e 

increases . 
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(5)  The  instability  can  be  supressed  by  the  presence  of  a  suffi¬ 
ciently  large  population  of  cold  electrons  if  T  is  not  too  high. 

In  passing,  we  remark  that,  besides  the  maser  instability,  a 
whistler  instability  may  also  occur  in  the  mirror  experiments  as  pointed  out 
by  Gladdi1  Experimentally,  it  should  not  be  difficult  to  distinguish  the 
maser  instability  from  the  whistler  instability,  even  if  both  should  be  opera¬ 
tive,  because  the  extraordinary  mode  can  be  unstable  at  all  angles  of  propaga¬ 
tion,  whereas  the  whistler  mode  is  most  unstable  for  parallel  or  nearly  paral¬ 
lel  propagation. 

Finally,  we  hope  that  more  observations  will  be  carried  out  in  the 
near  future  so  that  a  comparison  between  experimental  results  and  theoretical 
predictions  may  become  possible. 

6 .  Calculation  of  Spontaneous  Cyclotron  Emlssivlty  Using  the  Complete 
Relativistic  Resonance  Condition 

General  expression  for  the  emissivity  of  cyclotron  radiation  from 
temperature  plasmas  including  collective  effects  have  been  derived  and 
discussed  in  two  recent  articles  by  Freund  and  Wu^  and  by  Audenaerde^  (rele¬ 
vant  bibliographies  concerning  earlier  publications  are  cited  in  these  arti¬ 
cles)  .  Although  sufficient  generality  has  been  retained  in  these  results, 
practical  evaluations  of  these  expresions  are,  however,  difficult.  As  a 
result,  emissivity  has  been  discussed  in  two  distinct  regimes:  (1)  N  cos  0  » 
v/c  and  (2)  N  cos  9  «  v/c,  where  N  is  the  index  refraction  of  the  radiation, 
0  is  the  angle  between  the  wavevector  and  the  ambient  magnetic  field,  v 
denotes  a  characteristic  velocity  of  the  electrons,  and  c  is  the  speed  of 
light.  In  the  first  regime,  the  Doppler  effect  prevails  over  the  relativistic 
effect  in  the  wave— particle  resonance  condition  and  vice  versa  for  the  second 
regime.  It  is  evident  that  when  the  refractive  index  is  close  to  unity  the 
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later  regime  only  occurs  in  the  very  small  range  of  6  (i.e.,  8  -  90°)  for 
weakly  relativistic  electron  energies  (<  10  keV) .  Obviously  the  aforemen¬ 
tioned  approximations  become  invalid  in  the  intermediate  regime  where 
N  cos  8  -  v/c.  The  range  of  9  which  defines  the  intermediate  regime  increases 

with  electron  energy,  and  results  in  a  serious  limitation  of  the  previous 
calculations?®’®9 

Spurred  by  efforts  to  explain  the  auroral  kilometric  radiation®® 
considerable  progress  has  been  made  in  treating  the  general  resonance  condi¬ 
tion,  and  in  obtaining  general  expressions  for  the  growth  rates  of  the  escape 
modes  (ordinary  and  extraordinary)  in  the  presence  of  a  population  of  supra- 
thermal  electrons  with  a  loss-cone  distribution  functioni2-^®’ ’ ®®  The  suc¬ 
cess  of  these  efforts  provides  the  motivation  for  the  present  work  in  which  we 

28  29 

extend  the  previous  »  work  on  the  emissivity  to  treat  arbitrary  angles  of 
propagation  with  no  approximation  imposed  on  the  resonance  condition*  In  this 
work  we  specifically  consider  the  emissions  of  the  escape  modes  due  to  ener¬ 
getic  electrons  with  a  loss-cone  distribution.  The  applications  of  the 
results  of  our  calculation  to  the  study  of  both  the  auroral  kilometric  radia¬ 
tion  and  to  magnetic  mirror  machine  experiments,  however,  will  be  presented  in 
separate  articles . 

In  summary,  the  spontaneous  cyclotron  emissivity  has  been  calcu¬ 
lated  using  the  complete  relativistic  resonance  condition.  Unlike  prior  cal- 
culations  of  synchrotron  radiation®8’29  the  result  is  valid  over  the  entire 
angular  range  of  propagation  without  regard  to  the  bulk  energy  of  the  electron 
population.  We  are  specifically  interested  in  the  evaluation  of  the  radiation 
spectrum  from  an  energetic  population  of  trapped  electrons  having  a  loss-cone 
distribution  in  a  low-density  plasma.  The  particular  feature  of  such  a  plasma 
(i.e.,  we  <  is  that  a  relatively  high  radiation  intensity  is  found  in  the 
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X— mode  at  frequencies  just  above  cutoff  due  to  the  fundamental  gyroresonance 
which,  since  this  is  an  escape  mode,  can  readily  propagate  out  of  the 
plasma.  While  cyclotron  damping  due  to  the  backgroun  plasma  can  be  expected 
to  be  negligible  for  sufficiently  low-bulk  temperatures  of  background  plasma, 
amplification  of  the  radiation  can  occur  due  to  the  anisotropic  nature  of  the 
suprathermal  electron  distribution.  The  detailed  nature  of  this  instability 
has  been  amply  discussed  in  the  literature  dealing  with  auroral  kilometric 

ration^“18 

Finally,  it  should  be  observed  that,  while  the  overall  level  of 
spontaneous  synchrotron  radiation  is  greatest  for  the  X-mode  above  cutoff,  the 
detailed  angular  spectrum  is  more  complex.  In  particular,  the  X-mode  radia¬ 
tion  intensity  at  the  higher  gyroharmonics  actually  exceeds  that  near  the 
fundamental  (  oj  >  ,  oj  >  J  for  angles  of  propagation  nearly  perpendicu¬ 

lar  to  the  ambient  magnetic  field. 

^ •  Preferential  Stochas tic  Acceleration  of  Heavy  Ions  by  a  Coherent 

Electrostatic  Hydrogen  Cyclotron  Wave 

In  a  hydrogen  dominated  multi-ion  plasma,  the  minority  ion  species 
with  large  mass  to  charge  ratios  are  preferentially  accelerated  by  a  coherent 
electrostatic  hydrogen  cyclotron  wave  propagating  nearly  perpendicularly  to 
the  magnetic  field.  The  equations  of  motion  have  been  studied  analytically 
and  numerically  using  an  appropriate  Hamiltonian  formalism.  Surfaces  of  sec¬ 
tions  have  been  plotted  numerically  to  visualize  the  solutions  and  to  make 
comparisons  between  the  numerical  results  and  the  analytical  predictions. 
Stochasticity  has  been  found  to  onset  rather  abruptly  when  the  wave  amplitude 
exceeded  a  critical  threshold.  The  stochastic  region  has  an  upper  bound  as 
well  as  a  lower  bound  on  the  ion  kinetic  energy.  The  maximum  ion  kinetic 
energy  is  proportional  to  the  ion  mass  to  the  5/3  power.  The  agreement 
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between  the  theoretically  predicted  values  and  the  numerically  calculated 

values  is  quite  good  in  most  cases.  In  general,  the  numerical  v  .  is 

mm 

slightly  less  than  the  theoretical  value,  and  the  numerical  vmax  is  slightly 
larger  than  the  theoretical  value.  The  effect  of  increasing  the  wave 
amplitude  is  to  decrease  vmin  and  to  increase  vmax,  thus  broadening  the  sto¬ 
chastic  region.  A  diffusion  equation  which  describes  the  stochastic  accelera¬ 
tion  of  the  ions  has  been  derived  and  the  resulting  ion  distribution  function 
has  been  determined.  The  importance  of  this  process  to  the  acceleration  of  0+ 
ions  in  magnetosphere,  as  observed  by  satellites,  has  been  discussed.  We 
expect  similar  behavior  during  stochastic  acceleration  by  other  types  of 
waves.  Hence,  these  results  may  have  important  implications  for  radio  fre¬ 
quency  heating  experiments  in  tokamaks  and  other  devices. 

33 

In  a  recent  paper,  Singh  et  al.  make  a  comparison  between  the  sto¬ 
chastic  acceleration  of  ionospheric  ions  by  electrostatic  hydrogen  cyclotron 
waves  observed  in  their  numerical  simulations  and  the  predicted  acceleration 
of  these  ions  based  on  the  theory  of  stochastic  acceleration  by  Papadopoulos 
et__al;4  The  assertion  made  by  Singh  et  al??  that  the  heating  of  the  helium 
ions  observed  in  their  simulations  is  appreciably  larger  than  the  heating  pre¬ 
dicted  by  the  theory  of  stochastic  acceleration  by  Papadopoulos  et  al?^  is  the 
result  of  an  algebraic  error.  In  fact,  upon  correcting  the  error,  we  find 
agreement  beween  the  values  of  the  minimum  and  maximum  velocities  of  the  sto¬ 
chastic  region  calculated  from  the  theory  of  Karney33  and  Papadopoulos  et  al?^ 
agree  quite  well  with  the  numerical  simulations  of  Singh  et  al?3  and  Gaffey36 
in  most  cases. 
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An  expression  for  the  spectral  emissivity  of  spontaneous  synchrotron  radiation  for  a  plasma 
which  consists  of  both  thermal  and  suprathermal  electron  components  is  derived  using  the 
complete  relativistic  cyclotron  resonance  condition.  The  expression  is  valid  over  all  angles  of 
propagation.  The  result  is  applied  to  the  study  of  the  emission  of  radiation  from  an  energetic 
population  of  electrons  with  a  loss-cone  distribution  in  a  relatively  low-density  plasma  (i.e.,  the 
electron  plasma  frequency  is  less  than  the  cyclotron  frequency). 


PACS  numbers:  52.25.Ps,  52.60.  +  h 
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I.  INTRODUCTION  !  ! 


General  expressions  for  the  emissivity  of  cyclotron  ra¬ 
diation  from  high-temperature  plasmas  including  collective 
effects  are  derived  and  discussed  in  two  recent  articles  by 
Freund  and  Wul  and  Audenaerde2  (relevant  bibliographies 
concerning  earlier  publications  are  cited  in  these  articles). 
Although  sufficient  generality  has  been  retained  in  these  re¬ 
sults,  practical  evaluations  of  these  expressions  are  difficult. 
As  a  result,  the  emissivity  has  been  discussed  in  two  distinct 
regimes:  (1)  A  cos  <9>Z;/cand  (2);Vcos  0<£v/c,  where N is  the 
index  refraction  of  the  radiation,  6  is  the  angle  between  the 
wave  vector  and  the  ambient  magnetic  field,  u  denotes  a 
characteristic  velocity  of  the  electrons,  and  c  is  the  speed  of 
light.  In  the  first  regime,  the  Doppler  effect  prevails  over  the 
relativistic  effect  in  the  wave-particle  resonance  condition 
and  vice  versa  for  the  second  regime.  It  is  evident  that  when 
the  refractive  index  is  close  to  unity  the  latter  regime  only 
occurs  in  a  very  small  range  of  0  (i.e.,  0^ 90°)  for  weakly 
relativistic  electron  energies  (  <  10  keV).  Obviously  the  afore¬ 
mentioned  approximations  become  invalid  in  the  intermedi¬ 
ate  region  where  N  cos  0c^v/c.  The  range  of  6  which  defines 
the  intermediate  regime  increases  with  electron  energy  and 
results  in  serious  drawback  of  the  calculations  presented  in 
Refs.  1  and  2. 
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Spurred  by  efforts  to  explain  the  auroral  kilometric  ra¬ 
diation, 3  considerable  progress  has  been  made  in  treating  the 
general  resonance  condition  and  in  obtaining  general  expres¬ 
sions  for  the  growth  rates  of  the  escape  modes  (ordinary  and 
extraordinary)  in  the  presence  of  a  population  of  supra¬ 
thermal  electrons  with  a  loss-cone  type  of  distribution  func¬ 
tion.4"12  The  success  of  these  efforts  provides  the  motivation 
for  the  present  work  in  which  we  extend  the  previous1,2  work 
on  the  emissivity  to  treat  arbitrary  angles  of  propagation 
with  no  approximation  imposed  on  the  resonance  condition. 

In  this  paper  we  specifically  consider  the  emissions  of  the 
escape  modes  due  to  energetic  electrons  with  a  loss-cone  dis¬ 
tribution.  The  applications  of  the  results  of  our  calculation  , 
are  to  the  study  of  both  the  auroral  kilometric  radiation  and  1 
to  magnetic  mirror  machine  experiments.  In  the  case  of  the 
auroral  kilometric  radiation  (AKR),  the  most  commonly  ac¬ 
cepted  theory  involves  the  induced  emission  of  extraordin¬ 
ary  mode  radiation  due  to  the  presence  of  a  suprathermal 
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electron  component  with  a  loss-cone  distribution  function. 
Indeed,  observations  of  energetic  electron  distributions  in 
the  auroral  regions  confirm  the  existence  of  such  energetic 
populations. 13  However,  while  growth  rates  for  such  an  in¬ 
stability  have  been  calculated,  the  basic  level  of  spontaneous 
radiation  which  determines  the  initial  condition  has  not  been 
discussed.  The  present  work  addresses  the  issue  of  the  spon¬ 
taneous  emissivity  at  both  the  fundamental  and  higher  har¬ 
monics.  The  topic  is  highly  relevant  to  the  question  of  why 
the  observed  second  harmonic  AKR  is  often  more  effective 
than  the  theory  predicts,  since  the  growth  rate  is  one  order  of 
magnitude  smaller  than  that  for  the  fundamental  AKR.  We 
will  show  later  that  for  nearly  perpendicular  propagation  the 
spontaneous  emissions  of  the  X  mode  at  the  second  harmon¬ 
ic  is  much  stronger  than  that  at  the  fundamental.  This  find¬ 
ing  may  explain  the  question  mentioned  earlier. 

In  addition  to  the  theory  of  AKR,  the  discussion  of  the 
spontaneous  emission  presented  in  this  paper  is  also  relevant 
to  the  magnetic  mirror  machine  experiments  for  controlled 
fusion  research.  In  particular,  in  the  planned  MFTF-B  ex¬ 
periment  at  the  Lawrence  Livermore  National  Laboratory  it 
is  anticipated  that  the  suprathermal  electron  energies  will  be 
of  the  order  of  100  keV.  In  each  of  these  cases,  the  relativistic 
effects  on  the  electron  cyclotron  resonance  condition  is  im¬ 
portant  and  will  be  included  in  the  present  analysis  in  a  self- 
consistent  manner.  Finally,  we  note  that  the  relativistic  ef¬ 
fects  on  the  resonance  condition  have  recently  been  shown  to 
be  important  to  the  application  of  electron  cyclotron  current 
drive  in  tokamaks  even  for  weakly  relativistic  electron  ener¬ 
gies.  14 

The  organization  of  the  paper  is  as  follows.  The  general 
formulation  of  the  emissivity  is  given  in  Sec.  II,  along  with  a 
discussion  of  the  complete  resonance  condition.  The  emissi¬ 
vity  is  evaluated  numerically  in  Sec.  Ill  for  the  specific  case 
of  a  loss-cone  distribution  using  the  fully  general  resonance 
condition.  A  summary  and  discussion  are  given  in  Sec.  IV. 

II.  GENERAL  CONSIDERATIONS 

The  physical  configuration  under  consideration  is  that 
of  a  uniformly  magnetized  plasma  (B0  =  B0ez)  which  con¬ 
sists  of  a  thermal  background  as  well  as  a  population  of 
suprathermal  electrons.  Since  the  suprathermal  electrons 
are  assumed  to  have  a  much  lower  density  than  the  thermal 
species,  the  dielectric  properties  of  the  plasma  are  deter- 
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mined  primarily  by  the  thermal  background,  and  the  supra- 
thermal  electrons  may  be  considered  to  act  solely  as  a  source 
of  emission  for  the  normal  modes  of  the  plasma.  Within  the 
context  of  this  model,  the  emissivity  rj  (defined  as  the  power 
radiated  per  unit  volume  per  unit  frequency  and  per  unit 
solid  angle  subtended  by  the  wave  vector  k)  can  be  expressed 
in  the  form 

, .  to.,9 ) .  £i.  “LiV.  P±.«El  (j>uFM 

2-  C  2p  J 

xm  1  [y± 

X5{co  —  mf2c/y  —  t;,, ),  (1) 

where  (<y,k)  are  the  frequency  and  wave  vector  of  the  radi¬ 
ation,  k  =ke. ,  0  ^cos “ 1  (k^  /k  \  k,  —  k , ,  ns  denotes  the 
suprathermal  electron  density,  u==  p/m  =  yy  is  the  momen¬ 
tum  per  unit  mass,  j'ssfl  +  u2/c2)ul  =  (1  -  v2/c2)~l/2  is  the 
relativistic  factor,  Fs  is  the  suprathermal  electron  distribu¬ 
tion,  I2e==\eBr/mec\  is  the  electron  cyclotron  frequency,  vL 
and  u are  the  components  of  the  velocity  in  directions  per¬ 
pendicular  and  parallel  to  B0,  /?iJf  s=t;lfi|/c,  and  Jm  and 
are  the  regular  Bessel  function  of  the  first  kind  of  order  m 
and  its  derivative,  respectively.  In  addition, 

y  co  [(1— a~)(l  —  N  ±  J3  cos  0)—N2±  sin2  <9  ] 

Fte  AX  sin  6  (sin2  6  p ) 

(2) 

-  P  sin  6 

=  =(1  -  .v=  /?,  cos  9 )'  (3) 

/r^sin4  9  +■  4 \or/{2  2 )(1  —  a2)2  cos2  9,  (4) 

where  a-=cj] /or,  coe  is  the  background  electron  plasma  fre¬ 
quency,  S  -  ck  / a >  is  the  refractive  index,  and  the  subscript 
±  refers  to  the  ordinary  and  extraordinary  modes,  respec¬ 
tively.  The  indices  of  refraction  for  these  wave  modes  are 
given  by  the  Appleton-Hartree  dispersion  relation15 

•v:r  =  1  —  2a-(l  —  or-)/[2(l  —  «-)  —  (sin-  ^  pj ] , 

(5) 

where /?:=/?;/jy2. 

The  resonance  condition  described  by  the  delta  func¬ 
tion  in  Eq.  ( 1)  may  be  expressed  in  the  form 

ul/am  +(«',-  Um)2/bl„  =  1,  (6) 

where 
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c*  1  —  N 2  cos 
and 

bi=a2m/(  1  —  N 2  cos2  0 ). 

In  the  regime  applicable  to  the  escape  modes  (iV 2  cos2  9  <  1 ) 
the  reality  of  the  solutions  requires  that 

m~fl \/io2 >  1  —  N2  cos2  9,  (10) 

and  the  resonance  condition  describes  an  ellipse  in  (uL,u  ) 
space,  centered  at  (0,wm ),  with  extrema  at  ( am  ,um )  and  (0 ,um 

cno2 


+  )•  This  is  illustrated  schematically  in  Fig.  1(a).  Observe 

that  in  this  case  a~m  and  b  are  both  positive.  We  restrict  our 
consideration  to  that  section  of  the  ellipse  for  which  ut  >  0. 
From  the  resonance  condition  we  find  that  resonance  occurs 
for  uL  =  um,  where 
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However,  the  range  of  u}}  over  which  this  is  valid  is  not  infi¬ 
nite  and  must  be  restricted  to  within  the  extrema  of  the 
ellipse.  The  resulting  integral  over  «,,,  therefore,  is  defined 
by  the  range  where 

wlr  ±bm.  (12) 

In  the  opposite  regime,  in  which  N2  cos2  $>  1,  Eq.  (6) 
may  be  rewritten  in  the  form 

(u\\  -“m)2/bl  -u2/\ai\  =  1.  (13) 

where  we  observe  that  a2m  <  0  and  b  2m  >  0.  This  case  is  sche¬ 
matically  illustrated  in  Fig.  Kb)  in  which  it  is  evident  that  the 
resonance  curve  is  a  hyperbola  with  u ,  intercepts  at  um 
±  Hence,  resonance  is  possible  only  for  w,,  <i?m  —  bm 
and  u;]  >um  -f  bm. 

The  normal  modes  described  by  the  Appleton-Hartree 
dispersion  equation  are  shown  schematically  in  Fig.  2  for 
0  <  #<t/2.  The  ordinary  mode  (i.e.,  IV 2  =  N\  )  exhibits  a 
cutoff  at  co  =  coe  and  a  resonance  at  co  =  &>lh  {0 ),  where 
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FIG.  1.  Schematic  representations  of  (a)  elliptical  and  (b)  hyperbolic  reso¬ 
nance  curves. 
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FIG.  2.  Schematic  illustration  of  the 


magneto-ionic  wave  modes. 
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is  a  resonance  frequency  for  which  aIh(d]=£a,  cos  9  when 
coe<ne,  and  co]h(9)czf2e  cos  9  when  coe>fie,  It  should  be 
remarked  that  this  expression  becomes  invalid  in  the  limit  as 
9  —  77/2  since  ojlh  (9 )  —  0  and  ion  contributions  become  im¬ 
portant.  The  branch  above  the  cutoff  frequency  is  the  escape 
mode  f  Ar  2  <  1 )  and  is  usually  referred  to  as  the  O  mode,  which 
below  the  resonance  the  mode  is  called  the  electron  whistler 
(EW)  branch.  Similar  results  apply  to  the  extraordinary 
mode  (i.e.,  *V:  =  .V:_  )  which  displays  a  resonance  at 
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and  a  cutoff  at 

+  [  1  “f  Mco)/0  ^  1  2},  (16)  i 

The  branch  below  the  resonance  is  called  either  the  Z  mode 
or  the  slow  extraordinary  mode.  The  escape  mode  is  above 
the  cutoff  and  is  termed  either  the  fast  extraordinary  mode 
or  simply  the  X  mode.  We  shall  confine  our  attention  in  this 
work  to  the  escape  modes  for  which  N2  cos2  9  <  1  and  the 
resonance  curve  is  an  ellipse.  However,  it  should  be  noted  j 
that  while  .V :  >  1  over  much  of  the  Z  and  EW  modes,  it  is  still 
possible  to  describe  the  resonance  curve  as  an  ellipse  at  suffi¬ 
ciently  large  angles  of  propagation. 

The  conditions  under  which  the  complete  relativistic 
form  of  the  resonance  condition  must  be  employed  involve  a 
com plex  relationship  between  the  wave  frequency  and  angle  j 
of  propagation  and  the  bulk  energy  of  the  electrons.  In  gen¬ 
eral,  the  relativistic  corrections  become  important  whenever 
either  j 


iV  cos  9  5  (v)/2 c, 
or 

(tr)/2c-  £  1 1  —  m/2t/a>\. 


From  these  inequalities  [(17)  and  (18)]  it  is  evident  that 
relativistic  corrections  to  the  resonance  condition  become 
important  whenever  (1)  the  wave  propagation  direction  is 
nearly  normal  to  the  ambient  magnetic  field,  (2)  the  index  of 
refraction  of  the  wave  is  much  less  than  unity  (i.e.,  near  cut¬ 
off),  or  (3)  the  wave  frequency  is  near  the  cyclotron  harmon¬ 
ic. 

Finally,  it  should  be  remarked  that,  while  we  are  pri¬ 
marily  concerned  with  the  spontaneous  radiation  in  this  pa¬ 
per,  the  relativistic  effects  on  the  resonance  condition  also 
apply  to  the  study  of  the  growth  or  damping  of  the  coherent 
emissions. 

Utilizing  the  resonance  condition  (11),  we  obtain  the 
following  expression  for  the  emissivity  of  a  gyrotropic  distri¬ 
bution  Fs  (tfjpwJ, 

V  ±{co,0)  =  e2nscocN  EslIHL  6  f  P"  du„F 

[rn,Y±  JJb±)  +  (um/c)J’Jb^)]2 

XH[m2n]/co2  +  N2±  cos-0-  1),  (19) 

where  the  relativistic  factor  ym  =  (mOJco)  +  [ut/c)N  _ 
cos  9,  H  is  the  Heaviside  function,  and  b  _  =  mN  _ 
sin  9(um/ymc )  (1  —  N±  /?  cos  9)~'. 


III.  NUMERICAL  ANALYSIS 

Because  of  the  complexity  of  this  expression  for  the 
emissivity,  there  is  no  recourse  at  this  stage  but  to  turn  to 
numerical  integration  for  specific  choices  of  the  supra- 
thermal  electron  distribution  function.  We  choose  a  loss- 
cone  type  of  distribution. 

X exp(  -  Bj/ajj)  exp(  -  u2/a:),  (20) 

where  a;U  describes  the  parallel  and  perpendicular  momen¬ 
tum  spreads,  respectively.  Such  a  distribution  function 
might  be  used  to  model  the  hot  electrons  in  the  end  cells  of  a 
tandem  mirror  machine  or  a  population  of  energetic  elec¬ 
trons  trapped  in  a  dipole-like  magnetic  field,  such  as  might 
be  found  at  relatively  low  amplitudes  in  planetary  magnetic 
fields  or  the  solar  corona. 

Equations  (19)  and  (20)  have  been  evaluated  numerical¬ 
ly  for  the  escape  modes  in  a  low-density  plasma  [coe  </2,), 
and  the  results  are  shown  in  Figs.  3-5  for  the  choice  of  pa¬ 
rameters  such  that  coe  =0.1  f2e,  and  ai/c  =  a/c  =  0.14 
(corresponding  to  thermal  energies  of  5  keV).  Note  that  for 
this  choice  of  parameters  the  X-mode  cutoff  frequency  oc¬ 
curs  at  co a  ~  1 .0 1  nt ,  and  the  stop  band  between  &>uh  [6 )  and 
the  cutoff  frequency  ts  narrow.  Significant  levels  of  sponta¬ 
neous  emission  are  possible,  therefore,  because  of  the  funda¬ 
mental  X-mode  gyroresonance  (i.e.,  m  =  1)  at  frequencies 
oo>coM.  Such  an  interaction  (at  frequencies  co>flf),  how¬ 
ever,  occurs  by  means  of  the  contribution  of  Doppler  shift  to 
the  gyroresonance  and  consequently  vanishes  when 
k-B0  =  0  (i.e.,  9  =  90").  A  spectrum  is  shown  in  Fig.  3  in 
which  we  plot  vs  co/flt  at  9  =  60*  (where  17  is  in 

units  of  erg-sr  '-cm-3),  and  the  solid  (dashed)  line  repre¬ 
sents  the  X(O)  mode.  Evidently  the  X-mode  emission  near 
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FIG.  3.  Spontaneous  synchrotron  spectra  for  the  O  and  X  modes  at  $  =  60\ 
cj  / 1}  =0.!.  and  or  /c  =  a  /c  —  0. 14.  The  solid  (dashed)  line  represents 
the.XiO  mode. 

the  fundamental  gyroharmonic  is  peaked  at co—  1.016/2,.  In 
addition,  the  peak  decreases  in  intensity  by  about  an  order  of 
magnitude  and  broadens  in  width  fas  measured  by  the  full- 
width-half-maximum)  at  successively  higher  harmonics. 


FIG  4.  Angular  variations  in  peak  radiation  intensity  of  the  O  and  X  modes 
-  -  for<y,./#,  =  0. 1.  and  a./c  =  a./c  =  0.14. 
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FIG.  5.  Angular  variations  in  the  frequency  corresponding  to  the  peak 
emissivity  of  the  O  and  X  modes.  The  solid  (dashed!  line  corresponds  to  the 
XlO)  mode. 


Similar  results  are  found  for  the  O  mode  which,  at  each  har¬ 
monic,  is  found  to  peak  at  approximately  65%  of  the  X- 
mode  maximum  at  this  angle.  It  should  also  be  observed  that 
we  have  not  shown  the  entire  O-mode  peak  near  the  funda¬ 
mental  gyroresonance.  The  reason  for  this  is  that  as  waves 
with  frequencies  below  the  local  cyclotron  frequency  propa¬ 
gate  out  of  the  plasma  (i.e.,  through  regions  of  decreasing 
magnetic  field),  they  will  encounter  a  layer  in  which  the  wave 
frequency  is  sufficiently  close  to  the  local  gyrofrequency  that 
strong  cyclotron  damping  due  to  the  thermal  electrons  will 
occur.  As  a  consequence,  the  observed  radiation  spectrum  at 
these  frequencies  will  bear  little  resemblance  to  the  sponta¬ 
neous  emissivity.  The  X-mode  radiation  near  the  fundamen¬ 
tal  |<y  2  ),  as  well  as  both  X-mode  and  O-mode  radiation, 

will  in  the  course  of  propagation  out  of  the  plasma  encounter 
,  only  higher  harmonic  cyclotron  resonance  layers  for  which 
;  the  cyclotron  resonant  damping  is  relatively  weak.  Finally,  it 
!  is  implicitly  assumed  that  the  temperature  of  the  back¬ 
ground  electrons  is  sufficiently  low  that  they  do  not  partici- 
j  pate  in  the  X-mode  resonance  above  the  cutoff,  and  cylotron 
1  damping  of  these  waves  is  negligible.  This  is  true,  for  exam- 
(  pie.  for  thermal  energies  less  than  about  100  eV. 

The  variation  in  the  peak  of  the  emissivity  versus  angle 
of  propagation  is  shown  in  Fig.  4  for  the  X  mode  at  the  first 
three  harmonics  and  for  the  O  mode  at  the  second  and  third 
harmonics.  The  frequency  at  which  the  peak  occurs  is  shown 
in  Fig.  5  where  the  solid  (dashed)  line  refers  to  the  X(O) 
mode.  It  is  evident  from  the  figures  that  the  fundamental  X 
mode  emission  occurs  over  an  extremely  broad  angular 
range  (<9  S  70°)  with  a  small  peak  at  57°,  and  that  the  emis¬ 
sion  frequency  varies  little  with  angle.  In  contrast,  there  is  a 
very  strong  angular  dependence  in  the  emission  of  both  the 
_  X  and  O  modes  at  the  higher  harmonics.  In  particular,  the 
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X-mode  emissivity  is  maximum  for  perpendicular  propaga¬ 
tion  while  the  O  mode  peaks  for  propagation  at  oblique  an¬ 
gles  |(9~68°  for  the  parameters  considered). 

IV.  CONCLUSIONS 

In  this  work  the  spontaneous  cyclotron  emissivity  has 
been  calculated  using  the  complete  relativistic  resonance 
condition.  Unlike  prior  calculations  of  synchrotron  radi¬ 
ation,1-  the  result  is  valid  over  the  entire  angular  range  of 
propagation  without  regard  to  the  bulk  energy  of  the  elec¬ 
tron  population.  We  are  specifically  interested  in  the  evalua¬ 
tion  of  the  radiation  spectrum  from  an  energetic  population 
of  trapped  electrons  (i.e.,  having  a  loss-cone  distribution)  in  a 
low-density  plasma,  because  the  study  is  motivated  by  the 
research  of  the  auroral  kilometric  radiation  and  possible  ap¬ 
plications  to  the  mirror  machine  experiments.  The  particu¬ 
lar  feature  of  such  a  plasma  (i.e.,  <u,  <  )  is  that  a  relatively 

high  radiation  intensity  is  found  in  the  X  mode  at  frequen¬ 
cies  just  above  cutoff  due  to  the  fundamental  gyroresonance 
which,  since  this  is  an  escape  mode,  can  readily  propagate 
out  of  the  plasma.  While  cyclotron  damping  due  to  the  back-  i 
ground  plasma  can  be  expected  to  be  negligible  for  sufficient¬ 
ly  low  bulk  temperatures  of  the  background  plasma,  amplifi¬ 
cation  ot  the  radiation  can  occur  due  to  the  anisotropic 
nature  of  the  suprathermal  electron  distribution.  The  de¬ 
tailed  nature  of  this  instability  has  been  amply  discussed  in 
the  literature  dealing  with  the  auroral  kilometric  radi¬ 
ation.4-10  Detailed  discussion  of  the  present  results  and  their 
application  to  AKR  will  be  elaborated  in  a  forthcoming  arti-  ! 
cle.  | 

Finally,  it  should  be  observed  that  while  the  overall  lev¬ 
el  of  spontaneous  synchrotron  radiation  is  greatest  for  the  X  ' 
mode  above  cutoff  the  detailed  angular  spectrum  is  more 


complex.  In  particular,  the  X-mode  radiation  intensity  at  the 
higher  gyroharmomcs  actually  exceeds  that  near  the  funda¬ 
mental  \a}'zri',coSQ)M)  for  angles  of  propagation  near  the 
perpendicular  to  the  ambient  magnetic  field. 
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ABSTRACT 

1  ie  electron-cyclotron  maser  instability  for  energetic  electrons  with  a 
loss  cone  distribution  function  is  studied  tor  hot  electron  mirror  experi¬ 
ments.  The  instability  can  occur  at  all  angles  cf  propagation  for  a  wide 
range  of  parameters.  The  growth  rate  is  significantly  reduced  by  the  presence 
of  a  population  of  cold  electrons,  and  the  instability  can  be  suppressed  if 
the  density  of  the  cold  electrons  is  sufficiently  large  and  the  temperature  of 
the  energetic  electrons  is  not  too  high. 


In  a  recent  paper  by  Lau  and  Chu}  a  special  case  of  the  cyclotron  maser 
instability  associated  with  a  loss-cone  distribution  function  is  discussed  for 
hot  electron  mirror  experiments.  Their  analysis  is  restricted  to  the  case 
that  the  radiation  propagates  parallel  to  the  applied  magnetic  field  3Q  (k  = 
0)  and  moreover  emphasis  is  placed  on  the  limit  k  *  0.  Actually,  before 
Ref.  1  was  published,  the  theory  of  such  an  instability  had  been  extensively 
developed  by  numerous  authors2  9  (mostly  in  the  space  physics  community), 
although  the  research  has  been  motivated  by  the  investigation  of  the  auroral 
kiloraetric  radiation.  (Early  publications  which  discussed  the  physics  of  the 
maser  instability  have  been  cited  in  Ref.  1  and  are  not  repeated  here.) 

In  view  of  the  fact  that  strong  bursts  of  cyclotron  radiation  have  been 
observed  in  hot  electron  mirror  machines}0’11  it  is  pertinent  and  interesting 

to  study  the  cyclotron  maser  instability  for  such  experiments  more  fully  than 
was  done  in  Ref,  1, 


In  the  following,  we  report  the  results  of  our  investigation.  We  have 
carried  out  a  fairly  general  analysis  in  which  not  only  radiation  with  k ,  *  0 
and  k|(  *  0  is  considered,  but  also  a  population  of  cold  electrons  is 

included.  The  purpose  of  including  the  latter  effect  is  to  examine  the  possi¬ 
bility  of  whether  the  instability  can  be  stabilized  by  these  electrons. 

To  proceed  with  the  discussion,  we  consider  an  unperturbed  distribution 

i  9 

of  the  form 


Fe  (vV 
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exp 


(  -\)  > 
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where  u  is  the  electron  momentum  per  unit  mass;  a  denotes  the  thermal  spread; 
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is  the  component  of  u  perpendicular  to  the  applied  magnetic  field  Bq;  Z  is 
an  index  which  measures  the  steepness  of  the  loss“cone  feature;  and  ne  is  the 
density  of  the  energetic  electrons.  The  dispersion  equation  may  be  written  as 


det  |  Dij  (k,to)  |  =  0 


(2) 


where 
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and  the  superscript  s  on  Qij  denotes  the  particle  species.  In  the  present 
discussion,  ions  are  ignored.  We  shall  consider  that,  in  general,  there 
exists  a  population  of  cold  electrons  in  addition  to  the  main  component  which 
has  a  thermal  energy  of  several  tens  of  keV.  Since  we  are  interested  in  waves 
with  refractive  indices  N  <  1  and  frequency  to  close  to  the  electron  cyclotron 
frequency  in  calculating  Q?  .  (for  the  energetic  electrons)  we  may  ignore 
terms  associated  with  high  harmonic  numbers  (i.e.,  |nj  >  1)^  and  thus  obtain 
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For  the  cold  electrons,  Q?j  may  simply  be  written  as 


> 


2 

(o  Q 

_ e  e 

2  2 
o)(oj  -  0  ) 
e 


? 

n  oo 
c  e 

n  2 

t  oo 


o  . 


(5) 


(6) 


2  9 

In  Eqs.  (3),  (5),  and  (6),  to  =  4irnte-/m;  nt  is  the  total  electron 
density;  and  nQ  is  the  density  of  the  cold  electrons.  Eq.  (2)  has  been  solved 
numerically  and  the  results  are  presented  in  Figs.  1-5.  The  real  part  of  the 
frequency  is  close  to  the  electron  cyclotron  frequency  in  most  cases  and  the 
graphs  of  the  real  frequency  are  very  similar  to  those  in  Ref.  1.  Henceforth, 
we  discuss  only  the  imaginary  part  of  the  frequency.  The  growth  rate 
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normalized  by  the  cyclotron  frequency  is  plotted  in  Fig.  1  as  a  function  of 

9 

the  normalized  wavenumber  kc/ft  for  T  (=  ma  /2)  =  50  keV,  to  /SI  =  0.5  n  =  0 

and  £  -  1.  Four  values  of  0,  0°,  30°,  60°,  and  90°,  have  been  considered.  We 

note  that,  although  the  maximum  growth  rate  occurs  at  the  cutoff,  k  =  0,  in 

this  case,  the  range  of  wavenumbers  over  which  the  growth  rate  is  large 

increases  as  the  angle  of  propagation  0  increases  and  for  perpendicular 

propagation  the  growth  rate  is  nearly  constant  for  0<k<0.6$l/c. 

~  e 

In  Fig.  2,  the  normalized  growth  rate  is  again  plotted  as  a  function  of 
the  wavenumber  for  =  0.5;  nc  =  0;  £  =1;  T  =  30,  50,  100  keV;  and,  for 

perpendicular  propagation,  0  =  90°.  It  is  seen  that  the  growth  rate  and  the 
range  of  unstable  wavenumbers  increase  as  the  temperature  increases.  More¬ 
over,  for  T  =  100  keV,  the  growth  rate  peaks  at  k  =  0.63  ^0/c,  which  is  a 

significant  fraction  of  the  free— space  wavenumber. 

In  Fig.  3,  the  growth  rate  maximized  over  the  wavenumber  is 

plotted  versus  0  for  four  values  of  u^/S^,  0.1,  0.3,  0.5,  and  0.7.  It  is  seen 
that  as  “e/0e  increases,  ui5inax  becomes  progressively  more  independent  of  0. 
The  maximized  growth  rate  u-j_jraax  is  also  plotted  versus  ue/^g  in  Fig.  4  for 
three  different  energies,  T  =  30  keV,  50  keV,  and  100  keV.  For  each  energy, 

two  indices,  £  =  1  and  £  =  2,  are  studied  for  the  case  0  =  90°.  we  see  that 

the  peak  of  wi>max  increases  with  increasing  energy  and  with  increasing  loss- 

cone  index  and,  furthermore,  it  shifts  to  a  higher  value  of  oj  /Q  . 

e  e 

In  Fig.  5  we  plot  the  growth  rate  ok  as  a  function  of  k  for  three  values 

of  nc/nt,  0,  0.1,  and  0.3.  In  this  case  we  take  T  =  50  keV,  qj  /Q  =0.5 

e  '  e  ’ 

£  -  1,  and  9  =  0°.  The  case  ng  =  0  was  studied  by  Lau  and  Chui  We  notice 
that  as  nc/nt  increases  the  instability  diminishes  gradually.  For  nc/nt  >0.3 
the  unstable  waves  are  practically  suppressed.  However,  this  conclusion  is 


-6- 


only  valid  for  the  case  0  =  0°.  Unstable  waves  can  still  occur  at  other 

angles  as  shown  in  Fig.  6,  in  which  the  case  0  =  90°  is  investigated. 

In  Fig.  6(a)  we  consider  T  =  50  keV,  w  /ft  =  0.5,  and  1  =  1.  It  is  seen 

e  e 

that  the  maximum  growth  rate  occurs  at  a  certain  finite  wavenumber  k  rather 
than  at  k  =  0.  Even  when  nc/nj.  =  0.5,  the  instability  is  not  completely  sup- 
pressed.  The  importance  of  finite  wavenumber  k  is  enhanced  when  the  electron 
energy  T  is  increased.  This  fact  can  be  seen  in  Fig.  6(b)  in  which  T  = 

100  keV  is  considered.  In  this  case,  even  for  nc  =»  0,  the  maximum  growth 
rate  tu^  does  not  occur  at  k  =  0.  Moreover,  for  each  value  of  nc/nt,  the 

maximum  growth  rate  is  higher  than  that  for  the  case  T  =  50  keV. 

In  summary,  the  major  conclusions  of  our  preliminary  study  are:  (1)  the 

cyclotron  maser  instability  can  occur  at  all  angles  of  propagation  over  a  wide 

range  of  parameters;  (2)  in  general,  the  instability  at  large  0  (i.e.,  0  - 

90°)  is  more  serious  than  that  at  0  =  0°;  (3)  the  peak  of  the  growth  rate 

maximized  over  k  increases  with  the  temperature  T  and,  for  large  T,  occurs  at 

a  finite  k  rather  than  k  =  0;  (4)  the  range  of  unstable  values  of  co  /ft 

e  e 

broadens  as  T  increases;  and  (5)  the  instability  can  be  suppressed  by  the 

presence  of  a  sufficiently  large  population  of  cold  electrons  if  T  is  not  too 
high. 

Readers  should  be  cautioned  that  the  above  conclusions  are  for  the 
extraordinary  mode  only.  Generally,  the  growth  rate  for  the  ordinary  mode 
radiation  is  much  smaller.  Nevertheless,  the  ordinary  mode  cannot  be  stabi¬ 
lized  easily  by  the  presence  of  a  small  population  of  cold  electrons.  This 
and  some  other  relevant  results  will  be  discussed  more  fully  in  a  subsequent 
article . 

In  passing,  we  remark  that,  besides  the  maser  instability,  a  whistler 
instability  may  also  occur  in  the  mirror  experiments  as  pointed  out  by 
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maser 


13 

Gladd.  Experimentally,  it  should  not  be  difficult  to  distinguish  the 
instability  from  the  whistler  instability,  even  if  both  should  be  operative, 
because  the  extraordinary  mode  can  be  unstable  at  all  angles  of  propagation, 
whereas  the  whistler  mode  is  most  unstable  for  parallel  or  nearly  parallel 
propagation. 

Finally,  we  hope  that  more  observations  will'  be  carried  out  in  the  near 
future  so  that  comparison  between  experimental  results  and  theoretical  predic¬ 
tions  may  become  possible. 


ACKNOWLEDGMENTS 

The  authors  thank  R.  Ellis  for  informative  discussions.  The  present 
research  was  supported  in  part  by  the  Office  of  Naval  Research  under  grants 
N00014-82-K-0208  and  N000 1 4-8 4-C -0255  and  in  part  by  the  National  Aeronautics 
and  Space  Administration  under  grant  NGL21-002-005 .  One  of  the  authors 
(H.K.W.)  is  a  NAS/NRC  Research  Associate. 


-8- 


REFERENCES 


1.  Y.  Y.  Lau  and  K.  R.  Chu,  Phys.  Rev.  Lett.  _50,  243  (1983). 

2.  C.  S.  Wu  and  L.  C.  Lee,  Astrophys.  J.  230 ,  621  (1979). 

3.  L.  C.  Lee  and  C.  S.  Wu,  Phys.  Fluids  23,  1348  (1980). 

4.  C.  S.  Wu,  C.  S.  Lin,  H.  K.  Wong,  S.  T.  Tsai,  and  R.  L.  Zhou,  Phys.  Fluids 

2A,  2191  (1981). 

5.  H.  K.  Wong.,  C.  S.  Wu,  F.  J.  Ke,  R.  S.  Schneider,  and  L.  F.  Ziebell, 
J.  Plasma  Phys.  28,  503  (1982). 

6.  D.  B.  Melrose,  R.  G.  Hewitt,  and  K.  G.  Ronnmark,  J.  Geophys.  Res.  87_, 
5140  (1982). 

7.  R.  G.  Hewitt,  D.  B.  Melrose,  and  K.  G.  Ronnmark,  Aust.  J.  Phys.  35,  447 
(1982). 

8.  C.  S.  Wu,  H.  K.  Wong,  D.  J.  Gorney,  and  L.  C.  Lee,  J.  Geophys.  Res.  _87^, 
4476  (1982). 

9.  P.  B.  Dusenberry  and  L.  R.  Lyons,  J.  Geophys.  Res.  S7_,  7467  (1982). 

10.  T.  A.  Casper,  Y.  J.  Chen,  R.  Ellis,  R.  James,  and  C.  Lasnier,  LLNL  Report 

UCID  19783  (unpublished.  April  1983). 

11.  N.  T.  Glass,  in  Proc.  Workshop  on  Hot  Electron  Ring  Physics,  San  Diego, 
CA,  1981,  edited  by  N.  A.  Uckan  (Oak  Ridge  National  Laboratory, 
Oak  Ridge,  TN,  1982),  CQNF-811203,  Vol.  1,  p.  247. 

12.  R.  A.  Dory,  G.  E.  Guest,  and  E.  G.  Harris,  Phys.  Rev.  Lett.  14,  131 

(1965). 

13.  N.  T.  Gladd,  Phys.  Fluids  26,  974  (1983). 


FIGURE  CAPTIONS 


Figure  1 


Figure  2 


Figure  3 


Figure  4 


Figure  5 


Figure  6. 


.  Normalized  growth  rate  “± /^e  as  a  function  of  the  normalized  wave- 
number  kc/ft  for  T  «  50  keV,  u  fa  =  0. 5,  %  =  1 ,  n  =  0,  and  0  = 
0°,  30°,  60°,  90°. 

.  Normalized  growth  rate  uk/S^  as  a  function  of  the  normalized  wave- 
number  kc/fte  for  we/ne  =  0.5,  1=1,  0  -  90°,  nQ  =  0,  and  T  =  30, 
50,  100  keV. 

.  Normalized  growth  rate  ,  max/^e  maxilTlized  over  wavenumber  as  a 
function  of  the  propagation  angle  9  for  T  =  50  keV,  1  =  2,  nc  =  0, 
and  03e/ne  =  0.1,  0.3,  0.5,  0.7. 

.  Maximum  growth  rate  uijr!lax/^e  as  a  function  of  the  plasma  density 
parameter  for  9  =  90°;  %  =»  1,  2;  and  three  values  of  the 

energetic  electron  temperature:  (a)'  T  =■  30  keV,  (b)  T  =  50  keV, 

and  (c)  T  =  100  keV. 

,  Growth  rate  to .  /ft  as  a  function  of  the  wavenumber  kc/n  for  T  = 
x  e  e 

50  keV,  to^/fi^  =  0.5,  £  =  1,  9  =  0°,  and  three  values  of  the  cold  to 
total  electron  density  ratio  n  /nt  =  0,  0.1,  0.3. 

Growth  rate  u.  /S2  as  a  function  of  the  wavenumber  kc /S2  for  to  /fi  = 
x  c  e  e  e 

0.5;  i  =  1;  9  =90°;  ac/n^  =  0,  0.1,  0.3,  0.5;  and  two  values  of  the 
temperature:  (a)  T  =  50  keV  and  (b)  T  =  100  keV. 
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ABSTRACT 


The  theory  of  cross-field  ion  streaming  instabilities  is  applied  to  the 
parameter  regime  of  the  end  cells  of  the  upgraded  Tandem  Mirror  Experiment 
with  neutral  beam  injection.  The  cross-field  ion  drift  as  well  as  the  elec¬ 
tron  thermal  anisotropy  Tgi  >  T  provides  the  free-energy  which  drives  vari¬ 
ous  instabilities.  Three  instabilities,  a  nearly  perpendicular  propagating 
modified  two-stream  instability,  an  obliquely  propagating  ion-ion  streaming 
instability,  and  an  obliquely  propagating  electromagnetic  lower-hybrid  insta¬ 
bility,  have  been  identified.  The  first  two  waves  are  electrostatic  and  have 
the  largest  growth  rate.  For  the  actual  operation  conditions  of  the  Tandem 
Mirror  machines,  it  is  found  that  the  ion— ion  streaming  instability  has  the 
largest  growth  rate.  When  more  energetic  neutral  beams  become  available,  the 
two-stream  instability  may  play  the  dominant  role  in  the  stability  of  these 
devices . 
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I. 


INTRODUCTION 


Tandem  Mirrors  are  magnetic  confinement  devices  in  which  the  ions  in  a 

central  (solenoidal)  cell  are  prevented  from  leaking  out  the  end  by  means  of 

an  electrostatic  potential  maintained  by  a  denser  plasma  in  mirror  end 

1  2 

cells.5  The  basic  concept  has  been  successfully  verified  in  the  Tandem 

3 

Mirror  "Experiment  (TMX).  However,  it  was  also  apparent  from  the  initial  work 
that  the  simple  Tandem  Mirror  could  not  be  easily  scaled  to  a  reactor-size 
device,  because  of  the  severe  technological  demands  of  producing  and  main¬ 
taining  the  dense  mirror  plasma.  Alternative  solutions  were  thereby  investi- 

A  C  £ 

gated,  leading  to  the  development  of  the  concept  of  a  thermal  barrier.’ 

The  thermal  barrier  is  located  between  the  central  cell  and  the  plug  and 
allows  the  electrons  in  these  two  regions  to  be  maintained  at  different  tem¬ 
peratures,  thus  producing  the  confining  potential  with  a  smaller  plug  dens¬ 
ity.  In  order  to  generate  the  barrier,  the  density  in  this  region  is  reduced 
by  pumping  out  some  of  the  ions^  This  is  achieved  by  injecting  neutral  beams 
at  oblique  angles  to  create  "sloshing  ion  distributions,"  which  have  a  density 
minimum  at  the  center  of  the  barrier.  Futhermore,  a  hot  electron  population 
is  also  created  in  the  barrier  by  means  of  ECRH  to  reduce  the  power  transfer 
between  the  plug  and  central  cell  electrons;  the  hot  electrons  also  act  as  an 
MHD  anchor^’ ^  The  thermal  barrier  concept  is  presently  being  tested  in  the 
upgraded  TMX-U  experiment^  and  has  been  incorporated  into  various  reactor 
des igns ^ 

While  the  inclusion  of  such  complexity  into  the  end  cell  of  the  Tandem 
Mirror  may  improve  the  overall  confinement  significantly,  it  is  implicitly 
assumed  that  such  configurations  can  be  stably  maintained.  However,  the 
introduction  of  (neutral)  ion  beams  and  hot  anisotropic  electrons  in  the  bar¬ 
rier  region  suggests  a  variety  of  possible  microinstabilities ^  Injection  of 
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neutral  beams  to  form  the  sloshing  ion  distribution  allows  the  possibility  of 
various  streaming  instabilities i°» 1 1  The  formation  of  hot  anisotropic  elec- 
trons  will  also  give  rise  to  unstable  modes,  such  as  the  whistler^”^  and  the 
electron  cyclotron  maser  instabilities J6-20  Instabilities  resulting  from  a 
combination  of  these  two  effects  are  also  possible,  but  as  yet  have  not  been 
examined.  In  a  recent  paper,  we  have  investigated  such  a  combination  of  ions 
drifting  across  a  magnetic  field  in  which  the  background  electron  temperature 
was  anisotropic  in  the  environment  of  the  earth's  bow  shock21  In  that  situa¬ 
tion  the  ion  streaming  results  from  ion  reflection  at  supercritical  shocks  and 

the  electron  anisotropy  arises  from  adiabatic  compression  at  the  shock  sur- 
22 

face.  We  have  found  that  the  well-known  (low-beta)  modified  two-stream 
instability23  is  modified  by  finite  beta  effects24  and  further  enhanced  by  the 
temperature  anisotropy.  In  another  paper  the  nonlinear  development  of  the 
instability  in  a  finite  beta  plasma  with  isotropic  electrons  has  been  stud¬ 
ied?5 

Using  the  same  basic  model  (local  linear  theory,  homogeneous  plasma, 
unmagnetized  ions,  etc.),  we  extend  the  calculations  in  this  paper  to  the 
parameter  regime  of  the  Tandem  Mirror.  We  find  that  three  instabilities  can 
occur:  an  oblique  ion-ion  mode  which  has  the  largest  growth  rate  whenever  the 

ratio  of  the  drift  to  the  thermal  background  ion  velocities  is  not  too  large 
(typically,  when  ^0/v^0  <  4);  a  modified  two-stream  mode  which  has  the  largest 
growth  rate  when  v0/v£0  >  4  and  is  most  unstable  for  wave  Dropagation  nearly 
perpendicular  to  the  magnetic  field  and  an  electromagnetic  lower-hybrid, 
instability  with  a  peak  growth  rate  for  wave  propagation  between  60°  and  75° 
relative  to  §.  T  le  first  two  modes  are  electrostatic  and  their  real  fre¬ 
quencies  are  insensitive  to  the  electron  thermal  anisotropy.  The  modified 
two-stream  mode  and  the  electromagnetic  lower-hybrid  mode  are  the  low-beta 
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limits  of  the  kinetic  cross-field  streaming  modes  with  electron  temperature 
anisotropy  discussed  by  Tsai  et  al?* 

The  plan  of  this  paper  is  as  follows.  In  Section  II  the  basic  model  and 
assumptions  are  reviewed  and  the  dispersion  tensor  elements  are  given.  Numer¬ 
ical  solutions  over  a  range  of  parameters  relevant  to  Tandem  Mirrors  are  shown 
in  Section  III.  In  Section  IV  the  results  are  discussed,  and  the  summary  and 
conclusions  are  given  in  Section  V. 
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II.  BASIC  CONSIDERATIONS 


A.  Physical  Model  and  Assumptions 

The  main  purpose  of  this  work  is  to  apply  and  to  extend  the  theory 

of  cross-field  ion  streaming  instabilities  with  electron  thermal  anlso- 

21  24 

tropy,  9  to  the  Tandem  Mirror  Machines.  In  previous  studies,  the  relative 
electron  ion  cross— field  drift  was  the  result  of  ions  reflected  by  a  high  Mach 
number  quasiperpendicular  shock  wave?2  The  stability  analysis  was  done  in  a 
plane  where  inhomogeneities  could  be  neglected  and  only  one  species  of  ions 
was  Included  in  the  theory.  Although  the  source  of  the  cross-field  ion-drift 
and  the  plasma  parameters  in  the  present  work  are  different  from  those  in 
previous  papers,  the  same  physical  model  and  assumptions  are  adequate  for 
describing  ooth  problems.  In  the  following,  we  discuss  the  physical  model  and 
assumptions  for  the  problem  under  consideration. 

The  drifting  current  in  the  mirror  devices  comes  from  the  injection 
of  a  neutral  beam  into  the  plasma.  Under  the  present  operation  conditions, 
the  Mach  number  is  low  and  an  Ion-ion  instability  may  play  an  important 

role.  Thus,  we  have  to  include  both  ion  species  in  our  theory.  In  the  back¬ 

ground  plasma  frame  we  have  electrons  and  ions  at  rest  with  respect  to  drift¬ 
ing  ions  of  the  beam.  We  Introduce  a  coordinate  system  in  such  a  way  that  the 
z  direction  is  parallel  to  the  ambient  magnetic  field  5  and  the  beam 

injection  velocity  ^  and  the  wave  vector  &  lie  in  the  x-z  plane.  We  use  a 

linear  j-ocal  dispersion  relation  for  a  homogeneous  plasma.  The  effect  of 
electron  temperature  anisotropy  is  again  taken  into  account,  because  in  mag¬ 
netic  mirrors  the  electrons  are  preferentially  heated  perpendicularly  to  $  . 

o 

We  assume  that  the  unstable  waves  have  frequencies  much  higher  than 
the  ion  gyrof requency ,  so  that  the  two  ion  species  may  be  treated  as  unmagne¬ 
tized.  We  take  the  electrons  to  be  magnetized  and  include  electron  thermal 
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anisotropy  with  Tgi  >  Tg||  in  the  theory,  where  the  perpendicular  and  parallel 
directions  are  relative  to  the  magnetic  field.  We  will  return  to  these 
assumptions  in  Section  IV. 

B.  Dispersion  Tensor 

We  consider  an  infinite,  homogeneous,  collisionless  plasma, 
immersed  in  a  uniform  static  magnetic  field  5^,  consisting  of  three  compo¬ 
nents:  Maxwellian  electrons  at  rest  with  temperatures  T  ,  and  T  relative  to 

el  ell 

Bo  and  density  ng;  Maxwellian  ions  at  rest  with  temperature  Tio  and  density 
nio*  anfl  a  Ilian  ion  beam,  with  a  temperature  T^g,  density  n^g ,  streaming 

with  velocity  at  an  angle  a  with  respect  to  the  magnetic  field.  Hereafter, 

,  2  2 

K .  V,,  . 

we  consider  waves  with  the  following  properties:  y^  g  ^ 

°,S  2^2 

i2i<<  “  <<:  where  ^  and  are,  respectively,  the  ion  cyclotron 

f requency  and  ion  thermal  velocity  (background,  beam).  The  component  of  the 
wavevector  k  perpendicular  to  is  denoted  by  k,  =  k  sin  9.  In  such  a 

plasma,  the  ions  are  effectively  unmagnetized.  Moreover,  because  the  ions  are 

massive,  we  can  neglect  their  contributions  to  the  electromagnetic  perturba- 
tions . 


In  the  following,  we  write  the  elements  of  the  dispersion  tensor21 
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Plasma  Parameters 


Our  study  is  applied  to  the  end  cells  of  the  Tandem  Mirror 
device.  The  parameters  chosen  are  intended  to  reproduce  the  current  and 
future  operation  conditions  in  the  TMX-U  machine  and  are  the  following? 


Magnetic  Field:  Bq 

=  5  kG 

Electrons:  T  , 

el 

=  50-85  keV 

T 

ell 

=  17  keV 

ne 

=  5.5  x  1012  era-3 

a 

e 

-  8.8  x  10^®  rad/sec 

0) 

e 

=  1.3  x  10^  rad/sec 

*el- 

=  0.45-0.75 

0ell 

=  0.15 

Background  Ions :  T^o 

-  1.1  keV 

n, 

io 

=  2.7  x  1012  cm-3 

^i 

=4.8  xlO^  rad/sec 

Wio 

=  2.2  x  109  rad/sec 

eio 

=  0.005 

Sloshing  Ions : 


=0.7  keV 


T 

is  II  ,1 

Es  =  1/2mis  V2  =  10-30  keV 

=  2.7  x  1012  cm”2 

-  4.8  x  107  rad/sec 

=  2.2  x  109  rad/sec 

0.  =  0.003 

is 


III.  NUMERICAL  RESULTS 


The  numerical  analysis  of  the  unstable  waves  propagating  in  the  plasma 

described  in  the  previous  section  are  shown  in  Figs.  1  through  14.  In  these 

curves,  the  growth  rate  is  maximized  over  wavenumber  and  the  real  (to  )  and 

imaginary  (y)  parts  of  the  frequency  are  normalized  to  the  frequency  to  = 
i/2  LH 

“i ^  *  The  numerical  solutions  were  obtained  from  the  complete  dispersion 

equation  (1),  including  harmonics  up  to  n  =  ±  20.  To  facilitate  our  discus¬ 
sion,  we  classify  the  waves  in  two  modes,  called  mode  I  and  mode  II.  In  the 
following  we  describe  the  main  characteristics  of  each  mode  separately. 

For  mode  I,  the  assumption  of  magnetized  electrons  and  unmagnetized  ions 
.lS  satisfied.  For  this  mode,  the  values  of  ck/rn^  are  of  order  unity.  We  find 
that  in  the  region  of  maximum  growth  rate,  1  <  ck/oJe  <  2  for  the  whole  set  of 
solutions.  The  maximum  of  the  peak  growth  rate  is  situated  between  58°  and 
72°  relative  to  the  magnetic  field  direction.  In  Fig.  1  we  have  plotted 

“r^LH  and  Y//uLH  as  a  function  of  the  angle  of  propagation  9,  for  T  ,/T 

el  ell 

5  and  4,  Tls/Tlo  -  0.6,  T^/T^  -  15,  Vq/v1o  -  3,  a  -  45°,  and  n^/n.  - 

nis  ne  curves  show  that  the  growth  rate  and  the  real  part  of  the 

frequency  decrease  very  rapidly  with  decreasing  anisotropy  and  that  the  angle 

at  which  the  peak  of  the  maximum  growth  occurs  increases  as  T  /T 

el  el! 

decreases.  The  dependence  of  the  real  and  imaginary  parts  of  the  frequency  on 

the  drift  velocity  is  shown  in  Figs.  2  and  3.  In  these  curves  T  f /T  ,  =  5 

el  e  II  * 

nio/ne  “  nisy/ne  “  a  45°,  and  T  /T  =  15.  Figure  2  shows  to  /ooT „  and 

eii  10  r  LH 

y/a)LH  versus  9  for  Vo/v.Q  =  3,  4,  and  5.  Increasing  the  drift  velocity, 
increases  y  and  cor  and  also  the  angle  at  which  the  peak  of  the  growth  rate 
occurs.  In  Fig.  3  the  variation  with  drift  velocity  of  the  normalized  growth 
rate,  maximized  over  wavenumber  and  wave  propagation  direction  and  the  corres¬ 
ponding  real  part  of  the  frequency  are  shown.  The  maximum  value  of  y/u 

LH 
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7  and  9  =  64°. 


occurs  for  VQ/vio  =  7  and  9  =  64°.  For  VQ/vio  =  10,  we  find  9  =  62°  for  the 
angle  of  propagation  of  the  most  unstable  wave.  In  Fig.  4  we  analyze  the 

instability  as  a  function  of  the  beam  injection  direction.  For  these  curves, 
Tej/Tell  =  5’  Vc/vio  =  and  we  vary  the  anSle  a  e com  45°  to  90°.  The  remain¬ 
ing  parameters  are  the  same  as  in  Fig.  2.  The  growth  rate  peaks  at  o  -  62° 

and  it  does  not  change  significantly  when  a  increases  up  to  75°.  However, 
when  a  >  75°  the  growth  rate  decreases  more  noticeably.  At  a  =  90°,  y/ui  is 
nearly  95%  of  the  value  at  a  =  75°.  The  way  the  instability  is  affected  by 
the  variation  of  the  densities  of  the  background  and  beam  ions  is  shown  in 

5,  where  the  real  and  imaginary  parts  of  the  wave  frequency  are  plotted 
as  a  function  of  n^g/ne  for  the  values  of  wavenumber  and  propagation  angle 
which  maximize  y.  The  paramenters  are  the  same  as  in  Fig.  2,  except  for 
Vvio  =  ^  and  f°r  the  values  of  ion  densities,  which  now  vary  from  nj^s/ne  = 
0.25  up  to  0.60.  The  background  ion  density  can  be  calculated  for  each  curve 
from  the  condition  of  the  quasineutrality  n_^0/ne  =  1  -  nis/ng.  As  the  ion 
beam  density  increases,  the  growth  rate  increases  and  the  real  frequency 
decreases.  .■Je  also  find  that  the  angle  at  which  the  growth  rate  peaks 
decreases  with  increasing  nlg/ne.  In  Fig.  6  the  real  frequency  U)  (k)  corres¬ 
ponding  to  the  maximum  growth  rate  is  plotted  as  a  function  of  the  wavenum¬ 


ber  k.  Electromagnetic  effects  are  important  because  ck/a^  <  2.  In  Fig.  7  we 
show  the  real  part  of  the  arguments  of  the  plasma  dispersion  functions  for  the 
electrons  (zQ  =  wr/k,ve|),  background  ions  (?iQ  =  wr/kviQ),  and  ion  beam 
(^is^  =  l“r  “  ki‘^ol/kvis)-  For  seller  values  of  the  drift  velocity  and 
large  angles  of  propagation  the  instability  is  highly  kinetic,  since  z±o  and 
ki  1  are  approximately  unity. 

For  mode  II,  the  assumption  of  unmagnetized  ions  is  again  valid.  How¬ 
ever,  the  numerical  results  show  that  hio>uis  »  >  1.  So  the  electrons  are 
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now  only  weakly  magnetized  compared  with  the  ions  species.  This  mode  exhibits 
two  peaks  for  the  growth  rate  maximized  over  wavenumber.  One  is  located 

nearly  perpendicular  to  the  magnetic  field  (9  =  89.5°)  and  nearly  disappears 
for  v0/vi0  <  A.  The  second  peak  is  situated  at  an  angle  between  52°  and  70°, 
approximately,  and  has  a  larger  peak  growth  rate  than  the  nearly  perpendicular 
peak  when  v0/v£0  <  A.  In  Fig.  8  the  normalized  imaginary  and  real  parts  of 
the  frequency  are  given  as  a  function  of  the  wave  propagation  direction  for 
Tel/Teil  =  5,  A,  and  3.  The  other  parameters  are  the  same  as  in  Fig.  1.  The 
growth  rate  peaks  at  0  =  52°  and  increases  as  the  anisotropy  increases.  The 
real  part  of  the  frequency  is  very  insensitive  to  the  ratio  T  ,  /T  „  .  In 
Figs.  9  and  10  we  repeat  for  mode  II  the  study  done  In  Figs.  2  and  3  for 
mode  I,  taking  T^/T^  =  A.  The  peak  at  9  =  89.5°  decreases  very  rapidly  with 
decreasing  v0/v^0.  The  second  peak  moves  toward  the  perpendicular  direction 
to  the  magnetic  field  as  the  drift  velocity  increases.  In  Fig.  10  both  peaks 
are  analyzed  as  a  function  of  v0/vi0*  The  peak  at  9  =  89.5°  (curve  A)  is 

maximized  at  v0/v^0  =  6.8  and  the  less  oblique  peak  (curve  B)  at  v0/v^0  = 
A.  5.  When  v0/vio  -  A. 85,  both  peaks  have  the  same  amplitude.  Above  this 

value  of  drift  velocity,  the  nearly  perpendicular  wave  has  the  larger  growth 
rate  and  below  ,  the  less  oblique  peak  has  the  larger  one.  In  Fig.  11  the 

dependence  of  the  unstable  waves  on  the  beam  direction  is  studied,  taking 
Vo/vio  =  ^  with  the  reamining  parameters  the  same  as  in  Fig.  A.  The  less 
oblique  peak  (curve  B)  is  maximized  when  70°  <  a  <  75°.  The  growth  rate  of 
the  perpendicular  peak  (curve  A)  increases  with  increasing  beam  injection 
angle.  We  also  find  the  the  angle  9  of  the  maximum  growth  rate  moves  toward 
the  direction  of  the  injected  beam.  The  suggestion  of  an  ion-ion  interaction 
as  the  instability  mechanism  for  the  unstable  wave  with  a  maximum  growth  peak 
less  perpendicular  to  is  shown  in  Fig.  12,  where  the  real  frequency  and 
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growth  rate  maximized  over  wavenumber  and  propagation  angle  9  are  given  as  a 

function  of  n.g/ne,  for  T^/T^  =  4,  VQ/vio  =  3,  and  a  =  45°.  The  maximum 

growth  rate  occurs  for  equal  background  and  ion  beam  densities.  We  also  find 

that  the  angle  at  which  the  growth  rate  peaks  is  insensitive  to  the  variation 

of  nig/ne.  In  Fig.  13  the  same  analysis  done  in  Fig.  6  for  mode  I  has  been 

made  for  mode  II.  We  can  see  that  the  electromagnetic  corrections  can  be 

neglected  in  mode  II  since  ck/we  »  1.  Finally,  in  Fig.  14  we  plot  the  real 

oart  of  the  arguments  of  the  plasma  dispersion  function  for  the  three  plasma 

components.  The  perpendicular  mode  is  kinetic  in  nature,  since  oj  /k  v  1 

r  II  e  II 

when  0  =  89.5°. 

Summarizing,  we  have  observed  three  kinds  of  instabilities.  The  first, 

described  as  mode  I,  is  an  oblique  wave  which  has  a  relatively  small  growth 

rate  compared  to  the  other  two.  It  is  an  electromagnetic  wave  with  fre— 

quency  and  is  the  low  8  limit  of  the  smaller  growth  rate  peaks  in 

2 1 

Figs.  2  and  3  of  Tsai  et  al .  for  8  =  1.  It  is  instructive  to  emphasize  that 
this  mode  is  distinct  from  the  kinetic  cross-field  streaming  instability 
(KCFSI ) ,  becoming  stable  when  Tg^  =  Tg||  and  remaining  electromagnetic  even 
when  8  is  small.  The  basic  sources  of  free  energy  arise  from  the  electron 
temperature  anisotropy  and  ion  beam  streaming  energy.  In  mode  II,  there  is  a 
nearly  perpendicular  wave  which,  as  will  be  shown  in  the  next  section,  is  a 
modified  two-stream  instability.  As  discussed  by  Wu  et_  al?t  this  mode  is  the 
low  8  limit  of  the  KCFSI  and  corresponds  to  the  nearly  perpendicular  growth 
rate  peaks  in  Figs.  2  and  3  of  Tsai  et  al?^  The  third  wave,  described  in 
mode  II,  is  also  an  oblique  wave  and  has  the  largest  growth  rate  for  the 
actual  operating  conditions  of  the  Tandem  Mirror  Machines.  It  is  an  ion-ion 
streaming  wave.  Both  waves  in  Mode  II  are  electrostatic.  In  Section  IV,  we 
present  a  discussion  of  these  waves. 
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IV.  DISCUSSION 


Numerical  solution  of  the  complete  dispersion  equation  (1)  for  the 
parameters  of  the  end  cells  of  the  upgraded  Tandem  Mirror  Experiment  with  neu¬ 
tral  beam  injection  have  been  presented  in  the  previous  section.  Three  kinds 
of  unstable  waves  have  been  identified.  One  is  electromagnetic  and  two  are 
electrostatic.  In  order  to  give  a  better  interpretation  of  these  waves  we 
will  next  present  an  analytical  analysis  of  each  mode,  using  appropriate 

approximations  to  the  dispersion  equation,  based  on  the  results  of  the  previ¬ 
ous  section. 

We  start  with  mode  I.  For  convenience  we  normalize  the  frequency  to  the 


usual  whistler  frequency  ajQ=  cV^cosO  /<u*.  Because  the  electrons  are 


strongly  magnetized,  we  retain  only  the  terms  associated  with  harmonics  n  = 
0,  ±1  in  (1).  Introducing  the  index  of  refraction  N  =  ck/w,  using  the  small 
argument  expansion  of  the  Bessel  functions  and  the  large  argument  expansion  of 
the  electron  plasma  dispersion  function  for  n  =  ±  1,  the  elements  of  the  dis- 
persion  matrix  can  be  written  approximately  as 


co  oo 
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2  ^ 

Dxz  =  DZX  =  n2{  sin0cos0+ j(3ei-0e||  )sin0cos0-ig  tan0+  ^  __2_  e  ^ ^inOcos©  } 


o)  c  k 
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J. 

Dyz  =  ”  Dzy  =  ”  iN  1  JT  [1  +  zoZ(zo)]  sine  +  if  tan6  }  ; 


o  e  II 


where  is  given  in  Eq.  (2), 
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For  tor  »  | y  I  ,  the  real  and  imaginary  parts  of  the  frequency  can  be  cal- 
culated  by  the  standard  approach, 


Re  D  (5,uO  «  0 

(5) 

Im  D  (i£,a)  ) 

y  =  - - £ _ 

r~-  Re  D  (£,«  ) 

9<or  r 

where  D  denotes  the  determinant  of  . . 

For  the  real  part  of  the  frequency,  we  can  neglect  the  term  (<o2/u>2) 

2  2  2  2  ox 

(c  k  Aoe)  cos  0  in  Dxx  and  Dyy,  because  to2  «  m2  and  ck/o^  £  2.  Using  this, 

we  obtain  an  equation  similar  to  Tsai  et_al_2J  equation  (B6), 
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Equation  (6)  reduces  to  the  usual  whistler  dispersion  equation  when 

Tel  “  "ell  and  0  =  °°*  1,16  waves  we  ate  discussing  satisfy  neither  of  these 

conditions  and  the  numerical  solution  gives  wr/«  £  10-3  at  the  peak  growth 
rate.  The  expressions  for  the  growth  rate  are  very  similar  to  those  given  in 
Ref.  (21).  The  significant  differences  are  that  in  the  present  case,  6  «  1, 
two  ion  species  have  been  considered  and  ^  typically  makes  an  angle  a  =  45° 
relative  to  B,  whereas  in  the  case  treated  in  Ref.  (21),  0-1,  only  one 
species  of  ions  is  considered,  and  a  =  90°. 

Next,  we  analyze  Mode  II.  Since  these  are  low  frequency  waves  and  the 
normalized  wavenumber  ck/wfi  is  large  compared  to  unity,  we  treat  them  in  the 
electrostatic  approximation.  Neglecting  electromagnetic  corrections,  the  dis- 
persion  relation  becomes 
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Substituting  (1)  into  (8),  we  obtain 
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This  expression  can  be  simplified  by  using  the  large  argument  expansion 
of  the  plasma  dispersion  function  for  the  electron  cyclotron  harmonics  with 
n  *  0  and  the  large  argument  expansion  of  the  Bessel  function  IQ.  The  use  of 
these  expansions  is  justified  by  the  numerical  solutions.  For  instance,  from 
the  previous  section  it  is  possible  to  show  that  typically  p  =  10,  while 

Ui(s,o)  >  102*  Hence’ 

,  2.2 
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where  the  prime  denotes  the  derivative  of  Z  with  respect  to  the  argument. 

For  the  obliquely  propagating  waves,  zQ  is  small,  as  shown  in  Fig.  14 
and  the  dispersion  relation  (10)  can  be  approximated  by 


*  2“ 

D  (k,u>)  «  1  +  ——  -  - 
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el 


is 


(w  -  k  •  t'  )2 


(11) 


This  is  essentially  the  dispersion  equation  for  the  ion-ion  streaming  insta¬ 
bility.  >  The  solution  of  Eq.  (11)  ts  independent  of  the  electron  mass 

and  the  parallel  electron  temperature.  The  dependence  on  the  perpendicular 
electron  temperature  is  weak  because  the  second  term  in  (11)  is  small.  Fur¬ 
thermore,  the  only  angular  dependence  is  on  the  absolute  value  of  the  differ¬ 
ence  between  0  and  a,  rather  than  on  0  and  a  separately. 

it  can  be  shown  that  the  largest  value  of  the  growth  rate  occurs  for 
Wio  Uis  -  wi*  i,e,»  for  ni0  =  nis  =  n0/2  in  agreement  with  Fig.  12.  For 
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this  case,  eq.  (11)  is  easily  solved  and  gives  unstable  waves  if  |  £•$  |  < 

The  maximum  growth  rate  occurs  for  I  I  =  3  ^  m  ,  in  which  case 

co  =  a)i  (3  '2  +  i)/2. 

For  the  nearly  perpendicular  propagating  mode,  9  i  90°  and  zq  becomes 
nearly  unity,  as  can  be  seen  in  Fig.  14.  In  this  case  the  dispersion  equa¬ 
tion  (10)  becomes  similar  to  the  dispersion  equation  for  the  modified  two- 
22  24 

stream  mode  *  in  the  electrostatic  limit.  It  Is  again  independent  of  the 

parallel  electron  temperature  and  depends  only  weakly  on  the  perpendicular 

electron  temperature.  However,  a  significant  difference  Is  that  there  is  a 

dependence  on  the  electron  mass  in  this  case. 

*^e  emphasize  that  this  modified  two— stream  mode  can  occur  only  in  a  very 

limited  range  of  angles  9  <  90°.  As  9  decreases,  the  parallel  component  of 

the  wavevector  increases  and  zQ  becomes  smaller  than  unity,  as  shown  in 

Fig.  14.  The  applicable  dispersion  equation  then  becomes  (11),  which  gives 

the  ion-ion  streaming  mode.  On  the  other  hand,  if  0  -  9  ,  where 
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the  electron-ion  beam  interaction  is  comparable  to  the  ion-ion  interaction  and 
the  mode  becomes  a  modified  two-stream  mode. 

Figure  10  can  now  be  easily  interpreted.  The  nearly  perpendicular  mode 
(which  has  zQ  =  1)  is  more  strongly  electron  Landau  damped  than  the  oblique 
mode  (which  has  zQ  «  1),  while  the  difference  in  the  background  ion  Landau 
damping  on  either  of  the  modes  is  not  so  pronunciated  for  a  given  value  of 
^o^vio*  nn  t^ie  other  band,  the  electron  Landau  damping  is  rather  insensitive 
to  VQ •  vi0>  ln  contrast  to  what  occurs  for  the  background  ion  damping  and  ion 
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beam  growth  rate,  as  can  be  inferred  from  Fig.  14.  Therefore,  for  plasma 
parameters  such  that  at  smaller  values  of  VQ/v^0  the  electron  Landau  damping 
(Te)  is  larger  than  the  background  ion  damping  (y^  ) ,  we  expect  that  at  a 
sufficient  large  value  of  the  drift  velocity,  y  >  y  ln  this  way,  the 
nearly  perpendicular  mode,  which  is  governed  by  an  electron- ion  beam  inter¬ 
action  and  has  the  damping  mechanism  provided  essentially  by  the  electrons, 
will  have  the  greater  growth  rate  when  Vo/vio  is  sufficiently  large. 
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V. 


SUMMARY  AND  CONCLUSIONS 


In  this  paper,  we  have  studied  analytically  and  numerically  three  insta¬ 
bilities  associated  with  the  injection  of  a  neutral  beam  into  the  end  cells  of 
the  Tandem  Mirror  Machine.  These  unstable  waves,  which  are  excited  by  the 
relative  electron-ion  cross-field  drift  and  by  the  electron  temperature  aniso¬ 
tropy,  may  also  occur  in  the  perpendicular  bow  shock.  However,  since  the  two 
problems  deal  with  different  situations,  a  particular  unstable  wave  may  have 
distinct  relative  importance  and  the  growth  rate  may  be  maximized  for  differ¬ 
ent  values  of  the  parameters  appropriate  to  the  space  and  laboratory  plasmas, 
even  though  they  result  from  the  same  source  of  free  energy.  Also,  in  the 
present  work,  more  parameters  are  free  to  vary,  for  example,  the  beam 
injection  angle  a,  which  provide  more  ways  to  optimize  a  particular  insta¬ 
bility. 

For  propagation  oblique  to  the  magnetic  field  $,  we  find  two  unstable 
waves.  One  has  been  identified  as  an  electromagnetic  lower-hybrid  instability 
(EMLHI )  and  the  second,  as  an  electrostatic  ion-ion  streaming  instability 
(IISI). 

It  is  found  that  the  EMLHI  has  the  following  properties: 

(1)  The  growth  rate  y  is  maximized  for  60°  <  0  <  75°. 

(2)  y  and  increase  with  increasing  Tel/Teu •  The  mode  becomes  stable 

when  the  electron  temperature  anisotropy  is  reduced,  i.e.,  T  =  T  . 

ei  ell 

(3)  y  is  nearly  constant  with  respect  to  the  angle  of  beam  injection 
for  45°  <  a  <  75°  and  decreases  more  rapidly  for  larger  values  of  a. 

(4)  y  increases  and  decreases  with  increasing  nig  (density  of  beam 

ions ) . 
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(5)  The  instability  is  highly  kinetic  for  smaller  drift  velocities  and 
large  angles  of  propagation. 

The  IISI  has  the  characteristics  that: 

(1)  The  growth  rate  y  is  maximized  for  50°  <  9  <  70°  # 

(2)  The  growth  rate  is  the  largest  compared  to  the  other  two  inodes, 
when  VQ/vio  <  4. 

(3)  y  exhibits  a  peak  for  70°  <  a  <  75°. 

(4)  y  is  maximized  for  n^g  =  n^Q. 

For  propagation  nearly  perpendicular  to  we  find  an  electrostatic 

modified  two-stream  mode  (MTSI)  with  the  properties  that: 

(1)  The  growth  rate  y  is  maximized  for  9  =  89.5°. 

(2)  y  increases,  but  wr  remains  nearly  constant,  as  Tel/Te[j  increases. 

(3)  xhe  growth  rate  is  the  largest  compared  to  the  other  to  modes  when 

Vno  > 4- 

(4)  y  increases  with  increasing  vQ/vio. 

(5)  The  mode  is  kinetic  in  nature  because  to  /k..  v  ..  ~  1  when  9  -  89.5°. 

r  1!  e  I! 

The  parameters  have  been  chosen  to  describe  typical  operation  conditions 
of  the  TMX-L  device.  Since  presently  VQ/v^0  **  we  expect  the  ion-ion 

streaming  mode  to  play  the  dominant  role  in  the  stability  of  the  Tandem  Mirror 
Machines.  However,  new  experiments  are  likely  to  have  larger  values  of  the 
ratio  Vo/vio,  and  we  expect  that  the  nearly  perpendicular  two— stream  mode  may 
have  the  largest  growth  rate  in  the  future  experiments. 


-21- 


ACKNOWLEDGEMENTS 


The  authors  thank  Dr.  M.  Tanaka  for  helpful  discussions  about  the  compu¬ 
ter  codes  used  in  the  numerical  computation.  The  present  research  was  sup¬ 
ported  in  part  by  the  U.  S.  Office  of  Naval  Research  under  grants  N00014-82-K- 
0208  and  N00014— 84— C— 0255 ,  in  part  by  the  National  Aeronautics  and  Space 
Administration  under  grant  NGL2 1-002-005,  and  in  part  by  the  Conselho  Nacional 
de  Desenvolviraento  Cientifico  e  Tecnologico  of  Brazil,  CNPq . 


-22- 


REFERENCES 


1.  G.  I.  Dimov,  V.  V.  Zakaidakov,  and  M.  E.  Kishinevsky,  Fiz.  Plasiny  2,  597 
(1976)  [Sov.  J.  Plasma  Phys.  2,  326  (1976)],  In  Proceedings  of  the  Sixth 
International  Conference  on  Plasma  Physics  and  Controlled  Nuclear  Fusion 
Research ,  Berchtesgaden,  West  Germany,  1976  (International  Atomic  Energy 
Agency,  Vienna,  1977).  Paper  No.  C4. 

2.  T.  K.  Fowler  and  B.  G.  Logan,  Comments  Plasma  Phys.  Controlled  Fusion  2, 
167  (1977). 

3.  D.  L.  Correll  _et__al_. ,  Nucl.  Fusion  _22,  223  (1982). 

4.  D.  E.  Baldwin  and  B.  G.  Logan,  Phys.  Rev.  Lett.  43,  1318  (1979). 

5.  J.  Kesner,  R.  S.  Post,  D.  K.  Smith,  D.  E.  Baldwin  and  Y.  C.  Lee,  Nucl. 

Fusion  _22,  577  (1982). 

6.  J.  Kesner,  R.  S.  Post,  B.  D.  McVey  and  D.  K.  Smith,  Nucl.  Fusion  22j  549 
(1982). 

7.  T.  A.  Casper,  Y.-J.  Chen,  R.  Ellis,  R.  James  and  C.  Lasnier,  Lawrence 

Livermore  National  Laboratory,  Report  No.  UCID-19783,  April,  1983  (unpub¬ 
lished)  . 

8.  R.  W.  Conn,  V.  K.  Dhir,  N.  M.  Ghoniem,  D.  M.  Goebel,  S.  P.  Grotz, 
F.  Kantrowitz,  N.  S.  Kim,  T.  K.  Mau,  G.  W.  Shuy  and  M.  Z.  Youssef,  Nucl. 
Tech.  /Fusion  2_,  563  (1982). 

9.  MARS  Study  Group,  Lawrence  Livermore  National  Laboratory,  Report 

No.  TB-047  [GPO  No.  689-504],  March,  1983  (unpublished) 

10.  P.  L.  Auer,  R.  W.  Kilb  and  W.  F.  Crevier,  J.  Geophys.  Res.  76_,  2927 

(1971). 

11.  K.  Papadopoulos,  R.  C.  Davidson,  J.  M.  Dawson,  I.  Haber,  D.  A.  Hammer, 
N.  A.  Krali  and  R.  Shanny,  Phys.  Fluids  JAj  849  (1971). 

12.  C.  F.  Kennel  and  H.  E.  Petschek,  J.  Geophys.  Res.  71,  1  (1966). 


-23- 


13. 


J.  E.  Scherer  and  A.  W.  Trivelpiece,  Phys.  Fluids  10 ,  591  (1967). 

14.  A.  B.  Mihkailovskii,  Theory  of  Plasma  Instabilities.  Vpl.  _1_  (Consultants 
Bureau,  New  York,  NY,  1974),  p.  235. 

15.  N.  T.  Gladd,  Phys.  Fluids  _26,  974  (1983). 

16.  L.  C.  Lee  and  C.  S.  Wu,  Phys.  Fluids  _23_,  1348  (1980). 

17.  R.  G.  Hewitt,  D.  B.  Melrose  and  K.  G.  Ronnmark,  Aust.  J.  Phys.  35,  447 
(1982). 

18.  Y.  Y.  Lau  and  K.  R.  Chu,  Phys.  Rev.  Lett.  50,  243  (1983). 

19.  H.  K.  Wong,  C.  S.  Wu,  F.  J.  Ke,  R.  S.  Schneider  and  L.  F.  Ziebell, 

J.  Plasma  Phys.  28 ,  503  (1982). 

20.  H.  K.  Wong,  C.  S.  Wu  and  J.  D.  Gaffey,  Jr.,  Phys.  Rev.  Lett.  _52_  (1984), 
submitted. 

21.  S.  T.  Tsai,  M.  Tanaka,  J.  D.  Gaffey,  Jr.,  E.  H.  da  Jornada,  C.  S.  Wu  and 
L.  F.  Ziebell,  J.  Plasma  Phys.  _3l_  (1984),  in  press. 

22.  C.  S.  Wu,  D.  Winske,  Y.  M.  Zhou,  S.  T.  Tsai,  P.  Rodriguez,  M.  Tanaka, 

K.  Papadopoulos,  K.  Akimoto,  C.  S.  Liu,  M.  M.  Leroy  and  C.  C.  Goodrich, 
Space  Sci .  Rev.  37 ,  63  (1984). 

23.  J.  B.  McBride,  E.  Ott,  J.  P.  Boris  and  J.  H.  Orens ,  Phys.  Fluids  J_5_,  2367 
(1972). 

24.  C.  S.  Wu,  Y.  M.  Zhou,  S.  T.  Tsai,  S.  C.  Guo,  D.  Winske  and  K.  Papa- 
dopoulos,  Phys.  Fluids  _26,  1259  (1983). 

25.  D.  Winske,  M.  Tanaka  and  C.  S.  Wu,  Phys.  Fluids  _27_  (1984),  to  be 
submitted. 


-24- 


FIGURE  CAPTIONS 


Figure  1.  Mode  I  growth  rate  (dashed  curve)  maximized  over  wavenumber  and 
real  part  of  the  frequency  (solid  curve)  as  a  function  of  propaga¬ 
tion  angle  9  for  two  values  of  electron  temperature  anisotropy 

Tex/Te.  ‘  4'5>  Vvlo  ’  ’  l5'.  «  ‘  43°  >  Tlt/Tlo  -  0.6, 

and  n,  ^  =  n„* „  =  0.5  n  . 


Ho 


is 


Figure  2.  Mode  I  growth  rate  (dashed  curve)  maximized  over  wavenumber  and 
real  part  of  the  frequency  (solid  curve)  as  a  function  of  propaga— 
tion  angle  9  for  several  values  of  the  relative  electron-ion  beam 
drif t  velocity  normalized  to  the  background  ion  thermal  velocity 

Vvio  ■  3-4-5>  Tei/Tel  '  5;  nlo  ■  “Is  ■  °-5  ne:  «  '  45°: 

Tls/Tlo  '  0-6-  Te,/Tlo  '  15‘ 

Figure  3.  Mode  I  growth  rate  (dashed  curve)  maximized  over  wavenumber  and 

propagation  angle  9  and  the  corresponding  real  frequency  (solid 

curve)  versus  relative  electron-ion  beam  drift  velocity  V  /v.  . 

J  o  io 

The  remaining  parameters  are  the  same  as  in  Fig.  2. 


Figure  4.  Mode  I  growth  rate  (dashed  curve)  maximized  over  wavenumber  and 

propagation  angle  9  and  the  corresponding  real  frequency  (solid 

curve)  versus  beam  in  lection  angle  a  for  V  /v,  =  4.  T  /T  =5 

o  io  el  e  tl  9 

nio  =  nis  =  °-5  ne>  and  Tell/Tio  =  15 * 


Figure  5.  Mode  I  growth  rate  (dashed  curve)  maximized  over  wavenumber  and 


propagation  angle  9  and  the  corresponding  real  frequency  (solid 


curve)  as  a  function  of  the  ion  beam  density  n^/n  for  V  /v .  = 

S  0  O  I.  O 


A,  T  , /T 


5.  T 


15 
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Figure  6. 

Mode  I  real  part  of  the  frequency  versus  wavenumber  ck/u  cor- 

e 

responding  to  maximum  growth  rate  for  vQ/vio  =  3,5;  T  /T  „  = 
5  (dashed  curve),  4  (solid  curve);  nio  =  nis  =  0.5  ne;  and 

Te,l/Tio  =  15. 

Figure  7. 

Mode  I  values  of  I  <i»r  -  |  /kvig ,  cor/kvio,  and  «  /k,v  ,  cor¬ 

responding  to  maximum  growth  versus  angle  of  propagation  9  for 

Vvlo  ■  3'5;  Tei/Te«  ■  5i  "to  -  "is  -  0-5  V  «  ■  45°i 
T.l,/Tlo  ■  15‘ 

Figure  8. 

Mode  II  growth  rate  (dashed  curve)  maximized  over  wavenumber  and 

the  real  frequency  (solid  curve)  as  a  function  of  propagation 
angle  9  for  T^/T<(  -  3,4,5;  V0/vlo  -  3,  T,,  /Tlo  -  15;  „lo  -  nts  - 
0.5  ne;  and  a  =  45°. 

Figure  9. 

Mode  II  growth  rate  (dashed  curve)  maximized  over  wavenumber  and 

the  real  frequency  (solid  curve)  as  a  function  of  propagation 

angle  6  for  Vvio  '  3>4-5;  Tei/Te.  “  4i  Tel 'Tto  '  15i  "i„  *  "is  ' 
0.5  ne;  and  a  =  45°. 

Figure  10. 

Mode  II  nearly  perpendicular  (A)  and  oblique  (3)  growth  rates 

(dashed  curves)  maximized  over  propagation  angle  0  and  wavenumber 

and  the  real  part  of  the  frequency  (solid  curves)  versus  drift 
velocity  VD/vl0  for  T^/T,,  -  4,  T#,/Tlo  -  15,  nlo  -  „ls  -  0.5 

and  a  =  45°. 

Figure  11. 

Mode  II  nearly  perpendicular  (A)  and  oblique  (B)  growth  rates 

(dashed  curves)  maximized  over  wavenumber  and  propagation  angle  0 

and  the  corresponding  real  frequency  (solid  curves)  versus  beam 
injection  angle  a  for  VQ/vio  =  3,  T^/T^  =  5,  /1±q  =  15,  and 

nio  =  nis  =  °-5  ne* 

-26- 


Figure  12.  Mode  II  maximum  oblique  growth  rate  (dashed  curve)  and  real 
frequency  (solid  curve)  as  a  function  of  ion  beam  density  nig/ne 
for  v0/v£q  -  3,  Tej_/Tg||  =  4,  Te||/T£0  =  15,  and  a  =  45°. 


Figure  13.  Mode  II  real  part  of  the  frequency  versus  wavenumber  ck/u  for 

e 

several  values  of  the  drift  velocity  V  /v,  =  3  4  5:  T  /T 

o  lo  *  el  ell 

4  (dashed  curve),  5  (solid  curve);  T  „ /T.  =  15;  n.  *  n,  = 

ell  10  10  Is 

0.5  ne;  and  a  =  45°. 


Figure  14.  Mode  II  values  of  |  u>r  -  |/kvig,  wr/kvio,  and  w^/k^v^  cor¬ 

responding  to  maximum  growth  rate  versus  angle  of  propagation  0 

for  Vvio  =  4»5;  Tel/Tell  =  4;  Tell/Tio  =  155  and  °  =  45°- 
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ABSTRACT 


The  investigation  of  the  kinetic  cross-field  streaming 
instability,  which  is  motivated  by  the  research  of  collision¬ 
less  shock  waves  and  studied  by  Wu  et  al .  (1983),  is  discussed 
more  fully  in  the  present  work.  Since  in  the  ramp  region  of  a 
quasiperpendicular  shock  electrons  can  be  preferentially  heated 
in  the  direction  transverse  to  the  ambient  magnetic  field,  the 
inclusion  of  the  effect  of  the  thermal  anisotropy  on  the  insta¬ 
bility  associated  with  a  shock  is  both  desirable  and  necessary. 
The  present  study  has  found  that  can  significantly 

enhance  the  peak  growth  rate  of  the  cross-field  streaming  insta¬ 
bility  when  the  electron  beta  is  sufficiently  high.  Furthermore, 
the  present  analysis  also  improves  the  analytical  and  numerical 
solutions  obtained  from  the  earlier  theory  by  Wu  et  al.  (1983). 
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I.  INTRODUCTION 


In  a  previous  article  (Wu  et  al . ,  1983)  it  has  been  reported  that 
the  relative  electron-ion  cross-field  drift  in  a  perpendicular  shock 
wave  can  result  in  a  kinetic  instability  which  has  the  following  essen¬ 
tial  features:  (1)  the  characteristic  real  frequency  of  the  unstable  waves 
is  of  the  order  of  the  lower-hybrid  frequency,  (2)  the  unstable  waves  are  basic 
ally  whistlers  and  have  mixed  electrostatic  and  electromagnetic  polarization, 

(3)  the  instability  is  highly  kinetic,  and  (4)  unlike  the  usual  modified 
two-stream  instability  (McBride  et  al.,  1972;  McBride  and  Ott,  1972) 
this  instability  is  not  suppressed  by  electromagnetic  effects  when  the 
relative  electron-ion  drift  exceeds  the  Alfven  speed  of  the  plasma.  Sine.e 
the  free  energy  which  drives  the  instability  is  associated  with  the 


relative  cross-field  drift,  the  unstable  whistlers  must  inherently  have 

oblique  or  nearly  perpendicular  propagation.  It  is  anticipated  that  the 

growth  rate  should  diminish  for  waves  with  wave  vectors  nearly  parallel 

to  the  ambient  magnetic  field  BQ  .  The  fact  that  whistlers  can  be  am¬ 
plified  by  streaming  ions  was  first  pointed  out  by  Boss  (1970)  . 

The  instability  theory  is  motivated  by  the  research  of  collisionless 
shock  waves.  The  basic  physical  situation  is  described  and  discussed  in 
a  review  article  by  Wu  et  al .  (1984).  In  view  of  the  fact  that  in  a 
shock  front  electrons  can  be  heated  preferentially  in  the  direction 
transverse  to  the  ambient  magnetic  field,  it  is  anticipated  that  in 
this  region  the  perpendicular  electron  temperature,  Tg ,  ,  may  exceed 
the  parallel  electron  temperature,  ,  as  evidenced  by  in  situ  space¬ 

craft  observations  (Scudder  et  al. ,  1983).  Consequently  such  a  thermal 
anisotropy  can  provide  an  additional  source  of  free  energy.  The  per¬ 
tinent  question  is  how  does  the  anisotropy  affect  the  cross-field  streaming 
instability .  With  this  in  mind  we  feel  that  the  preceding  analysis  by 
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Wu  et  al.  (1983).  needs  to  be  generalized  to  include  the  effect  of 
Tej.  >  Te||  SO  that  the  stability  theory  can  be  impr0ved. 

In  passing  we  remark  that  in  the  preceding  work  by  Wu  et  al.  (1983) 

we  have  imposed  an  approximation.  That  is,  in  the  dispersion  equation 

we  have  only  retained  electron  terms  associated  with  n  =  0  (where  n 

is  the  cyclotron  harmonic  number) .  The  justification  is  that  we  are 

only  interested  in  waves  with  frequencies  oj  and  with  wavelengths 

2  2  2 

such  that  kj.veJL/fie  <<:  1  •  However,  numerical  solutions  of  the  disper¬ 
sion  equation  show  that  the  unstable  waves  can  have  wave— numbers  k.<  w  /c 

e 

This  implies  that  for  the  case  8  ^  1  we  have  k,v  ,  .0,  <  1  ,  which 

is  inconsistent  with  the  postulate.  As  a  result,  it  appears  that  in  order 
to  improve  the  numerical  solutions  we  need  to  consider  terms  associated 
with  n  ^  0  .  This  task  is  carried  out  in  the  present  paper  in  addition 
to  the  investigation  of  the  effect  of  the  thermal  anisotropy  on  the  in¬ 
stability  . 

The  organization  of  the  paper  is  described  as  follows:  in  Sec.  II 
the  basic  considerations  are  introduced.  In  this  section  the  physical 
model  and  assumptions  are  described,  the  dispersion  tensor  is  discussed  and 
the  plasma  parameters  for  the  problem  of  interest  are  listed.  In  Sec. 

Ill  a  numerical  analysis  of  the  complete  dispersion  equation  in  presented 
and  the  treatment  of  the  kinetic  cross-field  streaming  instability  of 
Wu  et_  al.  (1983)  is  generalized  to  include  the  effects  of  electron  tem¬ 
perature  anisotropy.  The  discussion  is  given  in  Sec.  IV.  Finally,  the 
summary  and  conclusions  are  provided  in  Sec.  V.  An  alternative  form  of 
the  dispersion  tensor  is  derived  in  Appendix  A  and  approximate  analytical 
expressions  are  obtained  in  Appendix  B. 


II.  BASIC  CONSIDERATIONS 


a.  Physical  Model  and  Assumptions 

In  the  present  study  the  physical  model  under  consideration  is  the 
same  as  that  described  by  Wu  et  al .  (1983,  1984)  .  in  addition  to  the 
cross-field  ion  drift,  the  electron  thermal  anisotropy  is  included 
in  the  discussion  of  this  paper.  It  may  be  helpful  to  comment  on 
the  assumptions  and  to  describe  the  model  briefly  so  that  the  sub- 

sequent  discussion  is  self-contained.  We  consider  a  one  dimensional 
perpendicular  shock.  To  facilitate  our  discussion  »e  introduce  the  fol¬ 
lowing  coordinate  system.  We  choose  the  t  axis  to  be  parallel  to  the  am-  \ 
blent  magnetic  field  BQ  .  All  macroscopic  quantities  vary  in  the  x  direc 

tion  normal  to  the  shock.  It  is  anticipated  that  a  cross-field  cur¬ 
rent  can  occur  in  the  shock  front  region  because  of  two  principal 
causes:  first,  the  presence  of  the  electrostatic  potential  barrier  and 
the  diamagnetic  effects  on  the  electrons;  and,  second,  the  reflection 
of  a  small  population  of  ions.  In  general,  the  relative  drift  between 
electrons  and  ions  can  occur  in  both  the  x  and  y  directions.  The  rela¬ 
tive  drift  in  the  y  direction  may  be  attributed  to  both  effects  just 
mentioned,  whereas  the  relative  drift  in  the  x  direction  is  only  caused 
by  the  transmitted  and  reflected  ions.  The  important  point  is  that 
there  exists  a  cross-field  current  in  the  y  direction,  but  the  plasma 
remains  current  free  in  the  x  direction  within  the  context  of  the  pres¬ 
ent  model,  as  shown  in  Fig.  1  of  Wu  et  al .  (1983).  Thus,  in  principle, 
the  plasma  density  inhomogeneity  may  be  significant  for  the  stability 
of  waves  propagating  in  the  y-z  plane  (Wu  et  al.,  1984;  Zhou  et  al., 

1983).  In  the  subsequent  discussion,  the  effects  of  plasma  inhomogene- 
lty  are  neglected.  Obviously  for  waves  propagating  in  the  x-z  plane. 
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this  assumption  is  generally  justified  (Wu  et  al . ,  1983;  1984).  More¬ 
over.  if  the  Eq  X  Bq  drift  prevails  over  the  diamagnetic  drift,  the 
gradient  of  density  and  temperature  may  be  of  secondary  importance  in 
the  stability  analysis,  even  for  waves  propagating  in  the  y-s  plane. 

We  postulate  that  the  unstable  waves  have  fluencies  and  growth  rates 
-ch  higher  than  the  ion  gyrof reguency  so  that  ions  may  be  treated  as  un¬ 
magnetized.  This  assumption  will  be  checked  a  posteriori.  In  the  elec¬ 
tron  frame  the  reflected  ions  and  the  transmitted  ions  have  opposite  drift- 
velocities  (Leroy  et  al . ,  1981;  1982).  Thus,  it  is  expected  that  there 
may  exist  two  types  of  unstable  waves;  one  is  excited  by  the  reflected 
ions  and  the  other  by  the  transmitted  ions.  Hence,  without  loss  of  gene¬ 
rality,  we  only  need  to  consider  one  ion  species  in  the  stability  analysis) 

because  only  one  ion  species  can  be  in  resonance  with  the  wave,  since 
the  transmitted  and  reflected  ions  have  oppositely  directed  drift  relo¬ 
cation  in  the  electron  frame.  Of  course,  this  leaves  out  an  ion-ion 

instability,  which  may  not  occur  in  a  high  Mach  number  shock  wave  any- 
way  (Wu  et  al . ,  1984). 


b .  Comments  on  High  Bata  Effects 

Whistlers  with  oblique  propagation  have  been  discussed  in  the  litera¬ 
ture.  However,  in  most  discussions  it  is  implicitly  assumed  that  the 
electron  beta  of  the  plasma  is  very  low  (i.e.  «  1)  .  Because  the 

subtlety  of  this  point  may  not  be  appreciated  at  first,  an  explanation 
appears  to  be  not  only  instructive  but  also  rather  desirable. 

It  is  usually  argued  that  to  study  whistlers  with  oblique  propagation, 
one  may  approximate  and  simplify  the  dispersion  tensor  by  neglecting  the 

electric  field  perturbation  along  the  magnetic  field,  <5E  .  [Readers 

z 

are  referred  to  the  discussion  given  on  pages  97  and  103  of  the  review 
article  by  Mikhailovskii  (1975)].  The  point  is  that  for  whistlers  the 
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2  2 

dielectric  tensor  element  e  ,  which  is  of  order  u)  /co  where 

zz  e  * 

2  2 

03^  =  47Tne  /m  is  large  in  comparison  with  the  other  elements  of  the 

2  2  2 

dispersion  tensor,  which  are  of  order  k  c  /to  .  This  is  so  for  whistlers 
because  the  condition 


2sn12  2 
0)  >>  k  c 

e 


2  2 

is  satisfied.  However,  it  should  be  pointed  out  that  E  *v.  0 (co  /a)  )  is 

Z  Z  6 

true  only  when  to/k,j  >  v^  ,  where  vg  is  the  electron  thermal  speed  and 
k|j  denotes  the  component  of  wave  vector  k  parallel  to  the  ambient 
magnetic  field.  Considering  the  usual  whistler  dispersion  relation 


oo 


2  2 

— ^  cos  8  , 

00 

e 


where  ^eB^/m^ci  is  the  electron  cyclotron  frequency  and  cos  0  = 
k-B^  /kB^  ,  we  find 


03 

kuv 


e 


ck  C^e  ck  1 


(2) 


For  0  -  1  ,  the  parallel  phase  velocity  co/k^  is  apparently  smaller  than 

the  thermal  speed  v^  .  In  this  case  the  aforementioned  argument  becomes 

2  2 

invalid  because  £  is  no  longer  of  0(toe  /to  )  .  In  fact,  one  can 
readily  show  that  in  this  case 


€ 

ZZ 


1 2^2 

k„A 

II  e 


2to 


j^e 


i  2  2 

k  jV 

II  e 


(3) 


where  X^  is  the  electron  Debye  length.  From  Eqs.  (1)  and  (3)  we  find 
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2.  2 

c  k 


=  2 


2n2 
c  SI 

_ e 

2  2 
v  to 

e  e 


2 


(4) 


This  result  implies  that  if  3  ~  1  ,  e  is  also  of  0(c2k2/(o2) 

6  Z  Z 

Evidently  the  approximation  of  neglecting  5E^  is  unjustified  and  the 
complete  3x3  dispersion  matrix  must  be  considered  for  a  high-beta 

4 

plasma.  This  is  the  task  which  we  carry  out  in  the  present  investigation. 


c  •  Dispersion  Tensor  D  ^  for  an  Anisotropic  Plasma 

To  proceed  with  our  discussion  we  first  obtain  the  expressions  for 

the  dispersion  tensor  elements.  In  doing  so  we  choose  to  work  in  the 

electron  frame  of  reference,  in  which  the  ions  are  considered  to  have  a 

cross-field  drift  V  .  For  convenience  we  assume  k*V  «  k  V  where  * 

-  -0  -  0 

the  wavevector k  lies  in  the  x-z  plane.  Because  ions  are  massive  their 
contributions  to  the  electromagnetic  perturbations  are  neglected.  More¬ 
over,  as  mentioned  earlier,  they  are  treated  as  unmagnetized.  In  the  follow 
ing,  the  ions  and  electrons  are  described  by  the  model  distribution  functions 


n . 

Fi  =  ~VT1  exp 


TT 


V  . 
1 


_1L 

2 

v. 

l 


o  -vny 

x  0 

2 
v . 

l 


(5) 


F  = 
e 


n 

“372  ~  exp 

TT  V  t,V 

eli  ex 


Vi, 


V 


2  1 


v. 


v  J 
ex 


(6) 


where  v±  ,  v^  and  veJ_  are  the  ion  and  the  parallel  and  perpendicular 
electron  thermal  velocities,  v2  =  2T^  /m^.  and  vd  represents  the  *E  x  B 

and  diagmagnetic  drifts.  Waves  propagating  in  the  x-z  plane  are  not 

affected  by  v , . 

d 
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The  expression  for  the  dispersion  tensor  D. ,  of  interest  for  the 
present  case  can  be  found  in  the  literature  (Stix,  1962;  Tsai  et  al., 
1981)  and  are  evaluated  in  Appendix  A  for  the  distribution  functions 
(5)  and  (6)  . 

Hereafter ,  we  shall  restrict  our  discussion  to  waves  which  have  the 
following  properties: 

t  2  2 

kivei 

- o—  <  1  and  <<  oo  <<  £2  /• 

20  -  1  e 

e 


0)  ±  0 

e 

k„  v 


en 


>  1 


(8) 


and 


(9) 


(10) 


Based  on  these  assumptions  we  can  estimate  the  order  of  magnitude  of 

various  terms  in  the  dispersion  equation.  We  shall  return  to  this  point  in 
Appendix  B. 


d .  Basic  Plasma  Parameters 

Since  the  present  study  is  mainly  motivated  by  the  research  of  the 
earth's  bow  shock,  the  physical  parameters  under  consideration  should 
be  compatible  with  those  of  the  solar  wind,  particularly  at  the  bow 
shock  region.  Thus,  we  shall  consider: 


10 


8e  =  6i  *  1 


0) 


n 


=  150 


Bn  s  7  Y 


1  £ 

T 

e|» 


£  2 


5va  <  V0  <  20va 


n  -  10  cm 
e 


-3 


(VA  the  Alfven  speed) 

(ne  is  the  electron  density)  . 


The  magnitude  of  the  cross-field  drift  VQ  depends  upon  its  nature  and 
direction.  The  range  given  here  is  based  on  recent  in  situ  observations  I 

by  spacecraft  (Paschmann  et  al . ,  1982)  and  computer  simulation  (Leroy  et  al.. 


1981  and  1982). 
(1984). 


Discussion  along  this 


line  is  also  presented  in  Wu  et  al. 
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III.  NUMERICAL  ANALYSIS  OF  THE  COMPLETE  DISPERSION  EQUATION 

In  this  section  we  present  a  numerical  treatment  of  the 
complete  dispersion  equation,  which  is  given  in  Appendix  A.  In 

Fig.  1  the  growth  rate,  maximized  over  wavenumber,  and  the  real  part  of 

the  frequency  are  plotted  as  a  function  of  the  angle  of  propagation  9  for 

several  values  of  the  ion  drift  velocity.  These  curves  have  general 

features  which  are  similar  to  those  obtained  by  Wu  et  al.  (1983)  for 

the  same  values  of  the  parameters',  6  =  B,=1,  T  /T  =1  and 

e  i  ex  en  * 

“e/fie  =  150 '  Xt  is  seen  that  increasing  the  ion  drift  velocity  to  several 
times  the  Alfven  velocity  does  not  suppress  the  instability,  but  rather 
shifts  the  direction  of  propagation  more  nearly  parallel  to  the  magnetic 
field.  This  result  may  be  attributed  to  retaining  the  electromagnetic  * 
terms  in  the  dispersion  equation.  A  significant  difference  between  Fig. 

1  and  Figs.  2  and  3  of  Wu  et  al .  (1983)  is  that  the  peak  growth  rate  in 
Fig.  1  in  this  paper  is  a  slowly  decreasing  function  of  the  drift  velo¬ 
city  for  VQ/vA  >  7.5  whereas  in  the  figures  of  Wu  et  al .  it  is  an  in¬ 
creasing  function  of  VQ/vA  up  to  Vq/va  =  20.  This  discrepancy  occurs 
because  in  the  present  calculation  we  have  retained  the  sum  of  electron 
cyclotron  harmonic  terms  with  n  f  0  .  In  general,  we  find  that  the  terms 
associated  with  the  cyclotron  harmonics  n  =  ±  1  in  the  dispersion 
equation  must  be  retained.  These  terms  give  rise  to  a  significant  modi¬ 
fication  of  the  growth  rate,  as  shown.  Moreover,  terms  with  |n|  ^2 
result  in  a  correction  of  order  20  %. 

The  combined  effects  of  electron  temperature  anisotropy  and  ion 
drift  are  illustrated  in  Fig.  2.  This  plot  is  the  same  as  Fig.  1,  except 

that  the  electron  temperature  ratio  T  /T  =  1.5  in  this  case.  It  is 

el  elt 

seen  that  the  electron  temperature  anisotropy  reduces  the  real  frequency. 
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However,  the  temperature  anisotropy  significantly  increases  the  peak 
growth  rate  and  shifts  the  corresponding  direction  of  propagation  to 
slightly  larger  angles.  These  effects  are  illustrated  more  clearly 
in  Fig.  3.  Here  the  growth  rate  is  plotted  as  a  function  of  0  for 
Vq/v^  =  10  and  T^/T  =  1.0,  1.5  and  2.0.  Decreasing  3  at  constant 
v0/vi  greatly  increases  the  growth  rate  (Wu  et  al . ,  1983);  however, 
increasing  the  temperature  ratio  Tg|/T  results  in  only  a  very  small 
increase  in  the  growth  rate  when  8  <<  1  . 

In  Fig.  4  the  growth  rate  and  real  frequency  are  plotted  as  fun¬ 

ctions  of  the  ion  drift  velocity  for  the  same  parameters  as  in  Figs. 

1  and  2  Both  the  growth  rate  and  the  real  frequency  increase  rapidly  with 
increasing  drift  velocity,  reach  a  peak,  and  then  decrease  rather  slowly 
after  the  peak.  It  is  particularly  important  to  note  that  the  peak 

value  of  the  growth  rate  occurs  for  v0^VA  ~  7‘5'  Increasing  the  tem¬ 

perature  ratio  Te ,  /Te„  increases  the  growth  rate  for  all  values  of  VQ/vA, 

but  increasing  the  temperature  ratio  does  not  change  the  value  Vn/v  PS  7.5 
where  the  peak  growth  rate  occurs.  It  is  also  seen  that  the  assumptions 
u)r  >  y  >  are  satisfied. 

The  real  frequency  03^  (k)  corresponding  to  the  maximum  grow7th 

rate  is  plotted  as  a  function  of  the  wavevector  k  in  Fig.  5.  It  is 

seen  that  the  assumption  ck/u)^  <<  1  is  net  satisfied  in  many  cases  and 

2  2  2 

the  terms  of  order  c  k  /co  must  be  retained  in  the  dispersion  equation, 

as  discussed  in  Wu  et  al.  (1983)  and  Appendix  B. 

The  arguments  of  the  plasma  dispersion  function  =  03  /k,r v 

0  r  II  ell 

and  jc^l  =  |  (wr  -  kjVg^kVil  corresponding  to  maximum  growth  rate  are 
plotted  in  Fig.  6  as  a  function  of  the  angle  of  propagation.  The  fol¬ 
lowing  two  conditions. 


are  satisfied  over  the  entire  range  of  9  .  This  is  a  further  indi¬ 
cation  of  the  kinetic  nature  of  the  instability  studied  in  this  paper 
For  the  maximum  growth  rate 


k  flv 


He|| 


0) 


r 


«  0,5  and 


and  the  assumption  (9),  Wr/knve||<  1  is  only  marginally  satisfied.  Be¬ 
cause  of  the  highly  kinetic  nature  of  the  instability  and  the 


importance  of  the  terms  cV/uJ  ,  it  is  not  possible  to  treat 


this  instability  accurately  in  a  simple  way.  However,  the  discussion 
in  Sec.  IV  helps  to  illustrate  the  physics  of  this  instability.  £ 

Another  important  characteristic  of  the  instability  is  that  it  has 
mixed  polarization,  with  electromagnetic  as  well  as  electrostatic  com¬ 
ponents.  The  ratio  of  electromagnetic  to  electrostatic  electric  field 
energy  density  is  plotted  in  Fig.  7  as  a  function  of  6  for  V_/v  =  5,10, 

A 

15,20  and  T^/T^j  =  1.0  and  1.5.  For  smaller  values  of  Vq/v^  >  the 
mode  propagates  more  nearly  perpendicular  to  the  magnetic  field,  and  the 
polarization  is  predominantly  electrostatic.  For  larger  values  of  VQ/vA 
the  mode  propagates  more  obliquely,  and  the  electromagnetic  component  of 
the  electric  field  becomes  comparable  to  the  electrostatic  component.  The 
ratio  of  the  electromagnetic  to  electrostatic  field  energy  density  in¬ 
creases  with  increasing  drift  velocity.  Increasing  the  temperature 
ratio  T  /T  increases  | 6E  |  / | 6E  |  slightly.  Although  the 
ratio  of  electromagnetic  to  electrostatic  electric  field  energy  density 
is  less  than  one  in  all  cases,  the  ratio  of  fluctuating  magnetic  field 
energy  density  to  electrostatic  electric  field  energy  density  can  be 
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very  large,  since 
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from  Ampere Ts  law  and 


c2k2/  |  oo  J  2 


is  very  large  because  the  phase 


velocity  of  the  waves  is  small  (Lemons  and  Gary  1977)  . 
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IV.  DISCUSSION 


The  numerical  analysis  in  the  previous  section  shows  that  an  elec¬ 
tron  temperature  anistoropy  t  >  T  results  in  a  substantial  increase 
of  the  peak  growth  rate  and  a  decrease  of  the  real  frequency  for  large 
angles  of  propagation.  However,  the  fundamental  nature  of  the  unstable 
waves  is  essentially  the  same  as  in  the  isotropic  case;  that  is,  they 
are  basically  obliquely  propagating  whistlers.  This  suggests  that  the 
graphs  in  Sec.  Ill  may  be  qualitatively  interpreted  using  the  theory  for 
the  isotropic  case. 

Before  proceeding,  a  significant  point  should  be  also  mentioned. 
That  is,  the  ion  dynamics  influence  the  wave  dispersion  even  when  the  r 
relative  electron-ion  cross-field  drift  is  small  and  negligible.  This 
can  occur  in  the  case  of  the  solar  wind,  for  which  the  ion  plasma  fre¬ 
quency  is  comparable  to  the  electron  cyclotron  frequency.  In  the  fol¬ 
lowing  we  write  the  dispersion  equation,  including  the  ion  dynamics  (Wu 
et  al.  1983) ,  as 


where 
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A 
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Since  the  curves  in  the  previous  section  were  calculated  for  the 


maximum  growth  rate,  we  use  the  ion  resonant  condition  oo^  =  k,VQ  to  obtain 


OJ 

r 


ULH 


ckV0  .  fl 
— - sm  8 

00  VA 

e  A 


(12) 


which  is  an  increasing  function  of  the  wavenumber,  of  the  drift  velocity 
and  of  the  angle  of  propagation.  This  equation  describes  how  the  real 
frequency  changes  as  a  function  of  (Vp/v^sin  0  ,  for  the  value  of.  k 
that  maximizes  the  growth  rate.  If  we  now  substitute  the  resonant  con¬ 
dition  into  Eq.  (11),  neglecting  the  effect  of  the  ions  for  simplicity, 
we  obtain 

Vn  m.  h  . 

—  sm  0  =  ( — )  —  cos  0  . 

v.  m  co 

A  e  e 

The  set  of  equations  (12)  and  (13)  are  convenient  to  describe  the 

nature  of  Figs .  1  to  6.  From  Eq.  (13)  we  infer  that  (V^/v^sin  0  and 

Vq/va  >  for  the  value  of  k  that  maximizes  the  growth  rate ,  are  both 

decreasing  functions  of  the  propagation  angle.  Therefore,  the  real 

frequency  is  an  increasing  function  of  the  drift  velocity.  In  Fig.  5, 

the  larger  the  ratio  Vq/v^  5  digger  va-lue  °f  the  frequency. 

Since  a)  /a)  increases  linearly  with  k  ,  as  shown  in  Eq .  (12),  our 

qualitative  estimate  is  in  good  agreement  with  Fig.  5,  as  well  as  with 

Fig.  4  when  Vq/v^  is  not  to°  large.  Furthermore,  from  Eq.  (13)  we 

also  deduce  that,  at  sufficiently  small  propagation  angles,  the  drift 

velocity  assumes  large  values  and  the  increase  of  the  parameter  V./v 

0  A 

dominates  the  increase  of  (Vq/v  )siii  0  as  we  decrease  the  angle  of 
propagation.  Therefore,  when  the  drift  velocity  is  large,  a  small 
decrease  in  0  can  produce  a  significant  variation  in  Vn/v  ,  without 

U  CL 


changing  the  value  of  (V0/vA)sin  0  .  In  this  way,  the  flat  region  of 
the  real  frequency  for  large  values  of  the  drift  velocity  shown  in  Fig. 

4  can  also  be  qualitatively  described. 

In  Figs.  1,  2  and  3,  the  real  frequency  is  plotted  as  a  function 
of  the  propagation  angle,  for  several  values  of  the  drift  velocity.  Since 
each  point  on  these  curves  has  a  different  value  of  k  ,  which  approxi¬ 
mately  satisfies  Eq.  (12),  it  is  convenient  to  combine  Eqs.(12)  and  (13) 
to  eliminate  k  .  This  gives 

i  t,2 

Wr  me  **  VQ 

aT“  =  —  sin  9  tan  9  *  <14> 

LH  l  v. 

A 

which  increases  very  rapidly  with  angle  at  a  fixed  V^/v^  »  as  shown 
in  Figs.  1,  2  and  3,  and  increases  with  Vq/v^  at  fixed  angle,  as  shown 
in  Figs.  1  and  2. 

In  Fig.  6  it  is  seen  that  |o>r  -  k^V^  |  /kv\  is  of  order 

unity,  even  though  co^/kv^  and  V^/v.  are  large.  This  shows  that  the 

ion  resonant  condition  is  nearly  satisfied  in  all  cases  ,  which  justifies 

the  use  of  Eq.  (12)  in  this  discussion.  The  quantity  w  /k„v  „  is  also 

r  it  ell 

of  order  unity  when  the  drift  velocity  or  the  propagation  angle  is  large. 
The  results  |  -  k^v^  | /kv_^  >  1  and  G0^/k||V^  <  1  clearly  indicate 

the  kinetic  nature  of  the  cross-field  streaming  instability.  Similar 
results  were  obtained  in  the  previous  study  of  Wu  et  al.  (1983).  Based 
on  the  arguments  presented  at  the  beginning  of  this  section,  we  could  have 
anticipated  that  the  inclusion  of  the  electron  thermal  anisotropy  results 
only  in  minor  modifications  of  the  curves,  as  shown  in  Fig.  6. 

Although  it  is  difficult  to  obtain  quantitatively  accurate  expres¬ 
sions  for  the  growth  rates  because  they  are  sensitive  to  small  changes 
of  the  variables,  the  qualitative  shape  of  the  graphs  of  the  growth 


rates  can  be  understood  easily.  For  very  small  angles,  the  mode  is 
stable  because  the  nearly  parallel  propagating  wave  is  unable  to  inter¬ 
act  efficiently  with  the  perpendicular  ion  drift.  As  the  angle  increases, 
the  ion  growth  initially  increases  faster  than  the  electron  damping  and 
an  instability  develops.  Eventually,  the  growth  rate  curve  turns  over 
and  decreases  with  increasing  angle  when  the  electron  damping  increases 
faster  than  the  ion  growth.  Finally,  the  mode  becomes  stable  at  an 
angle  less  than  90°.  These  effects  are  shown  clearly  in  Fig.  1,  2  and  3. 
Furthermore,  because  the  qualitative  analysis  for  the  maximum  growth  rate 
has  demonstrated  that  the  ratio  v0^VA  is  a  decreasin§  function  of  6  > 
we  expect  the  peak  growth  rate  to  move  toward  the  region  of  increasing 
propagation  angle,  as  the  drift  velocity  decreases,  in  agreement  with 
Figs.  1  and  2.  • 

In  Fig.  4  the  growth  rates  increase  initially  with  the  drift  velocity 

because  the  free  energy  which  feeds  the  instability  is  proportional  to 

2 

h  huVq  •  The  growth  rates  reach  a  maximum  at  V  /v^  -7.5  ,  and  then  de¬ 
crease  with  increasing  drift  velocity  because  the  angle  at  which  the 
growth  rate  peaks  decreases ,  which  increases  k^  and  the  electron  Landau 
damping . 

In  short,  this  discussion  qualitatively  explains  the  essential 
features  of  the  numerical  analysis  of  the  kinetic  cross-field  streaming 
instability  presented  in  Sec.  III.  The  most  important  results  are  that 
the  unstable  waves  have  been  identified  as  whistlers,  and  that  the  growth 
rate  is  largest  for  oblique  angles  of  propagation. 


V.  SUMMARY  AND  CONCLUSIONS 


We  have  investigated  both  analytically  and  numerically  the 

kinetic  instability  produced  by  a  relative  electron-ion  cross-field 

drift  and  an  electron  temperature  anisotropy  T  /T  >1  The 

ei  ell  ' 

unstable  waves  have  been  shown  to  be  whistlers  as  discussed  in  Sec.  IV. 

The  instability  is  not  suppressed  by  electromagnetic  effects  when  V  /v  > 

0  A 

and  8^1.  This  agrees  with  the  results  of  Lemons  and  Gary  (1977) .  As 
can  be  seen  in  the  numerical  results  in  Sec.  m,the  effect  of  increasing 
V0^VA  is  to  shift  the  direction  of  propagation  more  nearly  parallel  to 
the  magnetic  field.  The  effect  of  increasing  the  electron  temperature 
ratio  Tei/T  is  to  increase  the  growth  rate  and  to  shift  the  direction 
of  propagation  more  nearly  perpendicular  to  the  magnetic  field.  The  mods 
has  mixed  electrostatic  and  electromagnetic  polarization.  However,  for 


8^1  and  for  nearly  perpendicular  propagation  the  mode  is  predominantly 
electrostatic,  whereas  for  more  oblique  propagation,  the  electromagnetic 
and  electrostatic  components  are  nearly  equal.  Finally  for  8  small, 
the  mode  becomes  the  well  known  modified  two-stream  instability  (Wu  et  al., 
1983). 


In  summary  we  have  studied  the  plasma  instability  resulting  from  a 
cross-field  ion  drift  and  an  electron  temperature  anisotropy,  T  /T  >  1  . 
We  find  that  both  effects  significantly  affect  the  instability.  In  parti¬ 
cular,  the  present  study  has  led  to  the  following  major  conclusions: 

1.  A  relative  electron-ion  cross-field  drift  can  cause  instability, 

even  if  Vn/vA  >  1  . 

0  A 

2.  An  electron  temperature  anisotropy  T  /T  >  1  significantly  en— 

ex  ej| 

hances  the  peak  growth  rate, 

3.  The  unstable  waves  have  both  electrostatic  and  electromagnetic  com¬ 
ponents  in  general.  For  8^1  and  nearly  parallel  propagation, 


the  electrostatic  and  electromagnetic  components  are  comparable  and  the 
polarization  is  mixed. 

4.  The  instability  is  highly  kinetic  for  high  6^1.  From  numerical 
solutions  of  the  dispersion  equation  we  find  that 


Cd 

- — L_  <  i 
k»ve„  ' 


and 


kv. 

i 


>  1 


over  the  entire  range  of  angles  of  propagation.  In  addition,  electro¬ 
magnetic  effects,  high  order  electron  cyclotron  harmonics,  and  ion 
terms  all  affect  the  dispersion  equation  significantly  and  therefore  it 
is  difficult  to  treat  the  dispersion  equation  accurately  by  analytic 
techniques . 

Since  the  instability  is  driven  by  an  ion  —  wave  resonance,  the  * 

growth  rate  is  proportional  to  the  number  density  of  the  resonant  ions. 

Thus,  if  one  wants  to  apply  the  results  to  the  instability  associated 
reflected  ions,  he  must  bear  in  mind  that  the  growth  rate 
should  be  reduced  by  a  factor  n./n0.  If  n./n0  s  0.2,  such  an  approxi¬ 
mate  estimate  should  be  sufficient.  Correspondingly,  for  transmitted 
ions  the  growth  is  modified  by  a  factor  nfc/n0  which  is  of  order  unity. 
However,  in  other  situations,  the  real  frequency  as  well  as  the  growth 
rate  may  be  significantly  modified  by  the  densities  of  the  two  ion 
species  (Gaffey,  da  Jornada,  Winske,  and  Wu,  1984). 

In  addition  to  the  obliquely  propagating  kinetic  cross-field  streaming 
instability,  there  is  also  a  quasi-parallel  propagating  whistler  instabil¬ 
ity  which  can  be  excited  by  an  electron  temperature  anisotropy  T  /T  >1 

e-  e|| 

(Kennel  and  Petschek,  1966;  Sharer  and  Trivelpiece,  1967;  Mikhailovskii, 
1975).  The  maximum  growth  rate  ocoura  for  exactly  parallel  propagation,  and 
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the  growth  rate  decreases  very  rapidly  as  the  angle  of  propagation  in¬ 
creases.  This  mode  is  insensitive  to  ion  streaming,  and  it  does  not 
couple  to  the  cross-field  streaming  instability.  The  application  of 
the  parallel  whistler  to  the  earth's  shock  is  discussed  in  a  recent  article 
(Wu  et  al . ,  1984) , 
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APPENDIX  A 


In  this  appendix  we  derive  an  alternative  form  of  the  dispersion 
tensor,  including  the  ion  terms  and  the  electron  temperature  anisotropy, 
in  a  representation  based  on  the  electrostatic  and  electromagnetic  poten¬ 
tials  54>,6A  ,  rather  than  the  electric  field  6E  .  We  do  this  for  two 
reasons;  first,  to  make  contact  with  and  generalize  our  previous  analysis 
of  the  kinetic  cross-field  streaming  instability  (Wu  et  al . ,  1983),  which 
is  based  on  an  earlier  paper  (Hsia  et_  al . ,  1979),  and  second,  for  the 
greater  algebraic  simplicity  inherent  in  the  6<j>,5A  representation. 

This  point  will  be  discussed  further  below. 

Since  6E  =  —  ik6<J)  +  i(to/c)6A  and  k*6A  =  0  (Coulomb  gauge),  we  can 
write  * 


where  we  have  taken  k  =  (k^ , 0  ,k|| )  and  N  =  ck/to  .  In  the  following  we 
define  the  dispersion  matrix  |  |d  |  |  in  the  6E  and  represen¬ 

tations  such  that  Maxwell's  equations  become 
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One  can  readily  show  from  (Al)  and  (A2)  that 
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The  expressions  for  D  (i,j  =  x,y,z)  in  the  6E  representation  can 

be  written  in  terms  of  the  dielectric  tensor  such  that 
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The  electron  contribution  to  £_^_.  may  be  found  in  the  literature.  For 

the  b i-Maxwe 1 1 ian  distribution  function  given  in  Eq.  (5)  with  T  ^  T 

ei  e|| 

(e) 

£  „  has  been  derived  by  Stix  (1962)  and  recently  by  Tsai  et_  cil .  (1981)  , 
including  relativistic  effects.  Making  use  of  the  approximations-  that  the 
ions  may  be  treated  as  unmagnetized  and  neglecting  the  electromagnetic  ion 
terms  one  finds 
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I  is  a  modified  Bessel  function,  Z  is  the  plasma  dispersion  function, 
and  the  following  definitions  have  been  introduced: 
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Substituting  (A5)  into  (A3)  yields 
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Obviously  in  the  6<j),6A  representation,  the  dispersion  equation 
reduces  to  D  =  0  if  electromagnetic  effects  are  ignored.  In  other  words, 
the  advantage  of  writing  the  dispersion  equation  in  terms  of  *  * ' * 

(in  the  6(f),  6A  representation)  is  that  it  greatly  facilitates  the  dis¬ 
cussion  of  electromagnetic  effects  because  it  simplifies  the  separation  r 
of  the  electrostatic  and  electromagnetic  perturbations. 
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APPENDIX  B 


In  this  appendix  we  present  an  approximate  analytical  treatment 
of  the  kinetic  cross-field  streaming  instability. 

To  facilitate  our  discussion  we  introduce  a  quantity  which  is 

the  usual  whistler  frequency  for  6  g  «  1  ,  without  including  the  ion 
dynamics.  That  is. 
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Obviously,  it  is  convenient  to  normalize  the  frequency  u>  to  uQ  . 

It  not  only  provides  a  direct  comparison,  but  also  results  in  consider¬ 
able  algebraic  simplification,  as  will  be  seen  later.  In  the  following 
we  list  the  expressions  for  the  elements  of  the  dispersion  matrix  D 
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in  which  the  index  of  refraction  N(  =  ck/aj)  and  the  propagation  angle 

9[  =  cos  (k*B0/kB0)]  are  introduced.  If  we  retain  only  the  terms  of  order 
unity,  it  is  found  that 
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and  the  index  j  =  1,2,3.  In  these  expressions  we  have  defined: 
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and  Z(rj)  is  the  usual  plasma  dispersion  function. 
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Equation  (B2)  is  obtained  from  Eq.  (A5)  by  assuming  u  =  kfv  / 2fT  «  L 
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and  keeping  only  the  terms  associated  with  n  =  0  and  +  1  .  However,  we  should 
bear  in  mind  that  since 
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the  postulate  becomes  unjustified  when  8  .  >  ck/to  and  sin©  are  all 

el  e 


of  order  unity.  Moreover,  in  Eq.  (B2)  we  have  also  made  use  of  the 
assumption 
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We  remind  the  reader  that  Eq.  (B2)  was  derived  for  0(1).  For 

the  case  0  «  1  ,  the  term  (2N2ft»2/B  o7)  [1  +  z.Z(zn)]  in  D 

0  U  U  zz 

becomes  large  compared  to  the  other  elements.  As  a  result,  the  dis¬ 
persion  equation  may  be  reduced  to  a  2x2  determinant  (Mikhailovakii , 

1975).  In  short,  the  discussion  in  Sec.  lib  shows  that  for  (3^1 

e 

we  must  treat  the  3x3  determinant  in  order  to  derive  the  appropriate 
dispersion  equation.  This  is  the  principal  task  of  this  Appendix. 

In  the  following  we  write  w  =  oo^  +  iy  and  postulate  that 
>:>  Y  •  Furthermore,  since  we  assume  (jJ  «  (2  and 
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as  usual,  the  term  co  c  k  /u)  co  which  appears  in  D  and  D  may  be 

U  e  xx  yy  J 

considered  to  be  negligible.  Such  an  approximation  is  commonly  made 
in  the  literature  when  the  dispersion  equation  for  whistlers  is  dis¬ 
cussed  (Mikhailovskii,  1975).  Making  this  approximation,  we  find  by 
setting  the  real  part  of  the  determinant  of  D_  to  zero  that 
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where 
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From  this  expression  it  can  be  shown  that  increasing  the  electron 
temperature  anisotropy  Te±  >  decreases  the  real  frequency  when  the 

angle  of  propagation  is  large.  This  effect  may  be  seen  by  comparing 
Fig.  1  to  Fig.  2  and  is  also  shown  in  Fig.  3. 

For  the  case  of  isotropic  electrons,  Eq.  (B6)  may  be  simplified. 
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The  ion  dynamics  influence  the  wave  dispersion  even  when  the  relative 
electron-ion  cross-field  drift  is  small.  The  effect  of  the  ions  is 
is  greater  when  the  angle  is  larger. 

Next,  the  imaginary  part  of  u)  is  calculated  by  the  standard  approach. 
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Dyy*  However,  we  have  found  that  this  approximation  should  be  avoided 
in  deriving  the  imaginary  part  of  D,  the  determinant  of  D  ,  because  a  care¬ 
ful  study  has  shown  that  these  terms  can  qualitatively  modify  the  growth  rate 
at  sufficiently  large  0.  The  reason  is  explained  above  and  in  Wu  et  al.  (1983).. 
In  the  following  we  write 
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where  y  denotes  the  ion  contribution  to  the  growth  rate,  y  is 
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the  electron  Landau  damping  decrement,  and  Y  is  the  part  of  the 
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and  the  electron  Landau  term  for  T  =  T 
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Taking  f  =  f  and  g-  =  g  =  g  for  the  case  T  =  T  M  and  suppres- 
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sing  the  index  j  ,  the  electron  cyclotron  resonance  term  is 
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The  expression  for  0)^  is  given  in  Eq.  (B6)  and  the  function  A(k,(jO,0) 
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Ic  is  easy  to  see  that  one  of  the  reasons  the  electron  temperature  ani¬ 
sotropy  substantially  increases  the  growth  rate  y  is  because  the  functions 

f .  and  g  given  in  Eq.  (B3)  are  much  greater  for  T  >  T  than  for 
J  J  el  ell 

^ei  ^ell  ’  This  conclusion  is  in  agreement  with  the  numerical  solutions 
of  the  complete  dispersion  equation,  which  are  shown  in  Figs.  1,  2,  and  3. 


FIGURE  CAPTIONS 


Figure  1.  Growth  rate  (dashed  curve)  maximized  OTOr 

maximized  over  wavenumber  and  real 

part  of  the  frequency  (solid  curve)  as  a  function  of  propagation 


angle  0  for  several  values  of  the  relative  electron-ion  drift 
velocity  VQ/vA  =  5,10,15,20,  T^/T^  =  1.0  ,  ^  -  Q±  =  1, 

“e^e  =  150  *  1116  lower-hybrid  frequency  ui^  =  (JM^)*5. 


Figure  2.  Growth  rate  (dashed  curve)  maximized  over  wavenumber  and  real 
part  of  the  frequency  (solid  curve)  as  a  function  of  pro¬ 
pagation  angle  0  for  several  values  of  the  relative  electron- 
ion  drift  velocity  VQ/vA  =  5,10,20,  T@f/Te||  =  1.5,  6^=  8±  =  1, 

and  a)  /fl  =150.  £ 

e  i  T 


Figure  3.  Maximum  growth  rate  (dashed  curve)  and  real  part  of  the  fre¬ 
quency  (solid  curve)  versus  propagation  angle  for  several 
values  of  the  electron  temperature  anisotropy  T  (/T  =1,  1.5 

2  and  the  parameters  ft  j*R  =10  t,\  /O  -  „  / 

ett  Di  ±*u’  “  150 >  and  Vo^vA  =  10 • 

Figure  4.  Maximum  growth  rate  (dashed  curve)  and  real  part  of  the  frequency 

(solid  curve)  versus  relative  electron-ion  drift  velocity  V  /v 

0  A 

f°r  Tei/Te„  =  1'°  and  !-5.  ^  =  6.  =  1  ,  and  =  150. 

Figure  5.  Real  part  of  the  frequency  versus  wavenumber  ck/we  cor¬ 
responding  to  maximum  growth  for  several  values  of  V  /v 

0  A 

Tej/Tetl  =  1  (solid  curve),  T  /T  =1.5  (dashed  curve) 

6[|  9 

Be||  =  =  1  ,  and  =  150. 


Figure  6 


Figure  7 


Values  of  | |  /kv^  and  0>r /knve^  corresponding  to 

maximum  growth  versus  angle  of  propagation  0  for  several 

values  of  Vn/v  ,  T  /T  =1  (solid  curve),  T  /T  =  1.5 
u  a  ej.  e  e±  eu 

(dashed  curve),  8  .  =8.  =1  ,  and  a J  /ft  =  150. 
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ABSTRACT 


The  acceleration  of  electrons  along  the  auroral  magnetic  field  by 
large  amplitude  lower-hybrid  turbulence  and  a'  quasistatic  parallel 
electric  field  is  discussed.  Perturbations  are  excited  by  a  cross-field 
current  due  to  a  relative  electron-ion  drift.  The  fluctuating  electric 
field  is  primarily  electrostatic  and  propagates  nearly  perpendicularly 
to  the  auroral  magnetic  field.  It  is  assumed  that  the  instability 
saturates  rapidly.  The  time  evolution  of  the  energy  of  an  individual 
test  electron  is  calculated.  The  quasilinear  diffusion  equation,  with 
a  power  law  spectral  energy  density,  is  solved  using  a  Green's  function 
method.  The  time  evolution  of  an  initial  Maxwellian  distribution  is 
given,  and  the  average  electron  energy  is  obtained  for  various  cases. 

It  is  found  that  turbulent  acceleration  and  electric  field  acceleration 
are  comparable  in  some  cases. 
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I .  INTRODUCTION 


One  of  the  most  important  problems  in  auroral  physics  is  the  acce¬ 
leration  of  auroral  electrons  to  high  energies.  Experimental  measurements 
of  the  energy  spectra  of  the  precipitating  auroral  electrons  have  been 
made  by  the  ISIS— 1  and  S3— 3  satellites  which  have  orbits  passing  through 
the  auroral  region  at  2000  and  8000  km,  respectively.  These  data  suggest 
that  electrons  from  the  plasma  sheet  are  injected  into  the  auroral  region 
with  energies  of  order  1  keV  and  are  subsequently  accelerated  to  energies 
of  order  10  keV  before  precipitating.  The  results  from  the  S3-3  satel¬ 
lite  indicate  that  the  energization  of  auroral  electrons  along  the  auroral 
magnetic  field  lines  seems  to  occur  at  altitudes  between  4000  km  and  8000 
km  (Sharp  et  al.,  1980).  Furthermore,  the  presence  of  precipitating  elec¬ 
trons  with  energies  in  excess  of  10  keV  at  low  altitudes  in  the  auroral 
zones  during  disturbed  times  has  been  known  long  before  the  S3-3  measure¬ 
ments. 

Recent  electric  field  measurements  by  the  S3-3  satellite  indicate 
the  presence  of  a  quasistatic  electric  field  parallel  to  the  auroral  mag¬ 
netic  field  (Mozer  et  al. ,  1977,  1979,  1980;  Torbert  and  Mozer,  1978;  Cladis 
et  al.,  1979;  Sharp  et  al.,  1979;  Mizera  et  al.,  198Q) .  Such  an  electric 
field  may  provide  the  field  aligned  acceleration  process  required  to  explain 
the  observations  of  energetic  precipitating  auroral  electrons.  Although 
the  origin  of  the  parallel  electric  field  is  still  not  completely  under¬ 
stood,  theoretical  discussions  have  given  by  several  authors  (Chiu  and 
Schulz,  1978;  Kan  et  al.,  1979;  Shawhan  et  al .  ,  1979). 

In  addition,  wave  measurements  by  the  instruments  onboard  the  Hawkeye 
1MP6  satellites  have  observed  that  electrostatic  plasma  turbulence 
is  generated  in  the  region  where  the  quasistatic  parallel  electric  fields 
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occur  (Gurnett  and  Frank,  1972;  1977).  Observations  by  the  S3-3  satel¬ 
lite  also  confirm  the  existence  of  ion  cyclotron  and  lower-hybrid  waves 
in  this  same  region  (Kintner  et  al.,  1978;  1979;  Temerin,  1978,  1979;  Mozer 
et  al . ,  1979,  1980;  Temerin  et  al.,  1979).  The  lower-hybrid  waves  are  ge¬ 
nerally  broad  band  and  the  wave  intensity  is  relatively  high. 

In  the  literature,  several  schools  of  thought  have  been  proposed  to 
explain  the  acceleration  of  auroral  electrons .  A  brief  summary  of  these 
models  is  given  below. 

(1)  The  acceleration  may  result  from  the  presence  of  a  parallel  electric 
field  (Hoffman  and  Evans,  1968;  Gurnett,  1972;  Paschmann  et  al>>  1972. 
Wh-3-lsn  and  McDiarmid ,  19  72;  Ha.ll  and  Bryant,  19  74  5  Kaufmann  at  al.  ,  1976  5 
Bryant  et  al.,  1978;  Whalen  and  Daly,  1979;  Fridman  and  Lemair,  1980). 

(2)  The  acceleration  may  be  associated  with  electrostatic  shock  waves 
(Swift,  1970;  Kan,  1975;  Hudson  and  Mozer,  1978). 

(3)  The  acceleration  may  result  from  electrostatic  plasma  waves  (Gary 
et  al.,  1968;  Laval  and  Pellat,  1970;  Swift,  1970). 

The  fundamental  question  is  what  is  the  principal  mechanism  re¬ 
sponsible  for  the  acceleration  of  auroral  electrons.  Most  of  the  dicus- 
sions  described  above  are  qualitative  and  emphasize  only  the  basic  physical 
concepts.  No  detailed  analysis  of  the  characteristic  time  scale  of  a 
specific  acceleration  process  has  been  performed. 

The  purpose  of  the  present  paper  is  to  give  a  detailed  discussion  of 
the  acceleration  of  auroral  electrons  due  to  the  presence  of  both  a 
quasistatic  parallel  electric  field  and  enhanced  electrostatic  plasma 
turbulence.  We  are  primarily  interested  in  the  time  scale  which  cha¬ 
racterizes  the  rate  of  change  of  the  electron  energy  in  the  aforementioned 
situation.  The  following  discussion  is  primarily  based  on  an  individual 
test  particle  approach,  although  a  kinetic-theoretical  analysis  is  also 
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considered.  It  is  found  that  the  results  are  qualitatively  the  same. 

The  rest  of  the  paper  is  organized  as  follows:  In  section  -I >  the 
basic  assumptions  and  the  physical  model  are  described  and  the  mathe¬ 
matical  problem  is  formulated.  In  section  III,'  the  electron  velocity  and 
kinetic  energy  are  calculated  for  an  individual  test  particle.  A  Green's 
function  technique  is  used  in  section  IV  to  obtain  the  electron  velocity 
and  kinetic  energy,  and  the  time  evolution  of  an  initial  Maxwellian  dis¬ 
tribution  is  given.  Finally,  the  discussion  and  conclusions  are  given 
in  section  V.  The  Green's  functions  used  in  section  IVC  are  derived  in 
the  Appendix. 
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II.  GENERAL  CONSIDERATIONS  AND  THE  PHYSICAL  MODEL 

A.  The  Quasistatic  Parallel  Electric  Field 

The  measurements  by  the  S3-3  satellite  at  altitudes  of  approximately 
one  earth  radius  indicate  that  occasionally  the  maximum  value  of  the  pa¬ 
rallel  electric  field  can  be  as  large  as  200  mV/m  to  400  mV/m  (Mozer 
et  al.,  1977).  The  average  value  of  the  parallel  electric  field  is 
typically  less  than  10  mV/m  (Temerin  et  al.,1982).  However,  the  subject 
remains  controversial . 

A  number  of  authors  have  used  the  S3-3  measurements  of  the  electron 
energy  spectrum  and  pitch  angle  distribution  to  estimate  the  potential 
difference  along  the  auroral  magnetic  field  lines  (Sharp  et  al.,  1979; 
Cladis  et  al.,  1979).  It  is  found  from  the  plasma  data  that  the  signa¬ 
ture  of  an  electrostatic  potential  drop  exists  over  a  broad  altitude 
range,  from  approximately  4000  km  to  8000  km.  The  peak  value  of  the 
total  potential  difference  is  estimated  to  be  between  4  kilovolts  and 
5  kilovolts. 

B .  The  Quasilinear  Diffusion  Equation 

The  electric  field  and  plasma  turbulence  measurements  outlined 
above  suggest  that  in  the  acceleration  region,  the  electrons  are  energ- 
gized  under  the  combined  influence  of  both  a  quasistatic  parallel  elec¬ 
tric  field  and  large  amplitude  electrostatic  plasma  turbulence.  A 
theoretical  study  of  the  collective  instability  mechanism  of  electron 
acceleration  is  essential  not  only  for  the  understanding  of  the  ano¬ 
malous  energization  of  auroral  electrons,  but  is  also  Important  in  the 
context  of  fundamental  plasma  physics.  The  statistical  acceleration 
of  electrons  and  ions  in  turbulent  plasmas  has  been  discussed  by  many 
authors  and  comprehensive  reviews  have  been  given  in  the  literature 
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(Tsytovich  1966,  1977).  Recently  the  linear  and  nonlinear  evolution 
of  plasma  instabilities  driven  by  cross-field  currents  and  the  resulting 
anomalous  heating  of  electrons  and  ions  have  been  studied  analytically 
and  numerically  (Davidson  and  Krall,  1977).  These  works  show  that  sta¬ 
tistical  acceleration  processes  can  play  an  important  role  in  space 
plasmas  and  laboratory  plasmas. 

The  physical  situation  which  is  treated  in  this  paper  may  be  de¬ 
scribed  briefly  as  follows.  We  consider  a  plasma  that  can  support  an 
instability  that  amplifies  perturbations  excited  by  a  cross-field  current, 
due  to  a  relative  electron- ion  drift.  Among  the  instabilities  driven  by 
a  cross-field  current,  the  modified  two-stream  instability  and  the  lower- 
hybrid-drift  instability  are  the  most  important.  A  rather  complete  de¬ 
scription  of  the  literature  of  cross-field  instabilities  is  contained 
in  a  recent  paper  (Hsia,  et  al.,  1979),  which  gives  a  unified  theory 
of  the  lower-hybrid-drift  and  modified  two-stream  instabilities.  These 
modes  with  frequencies  in  the  range  of  the  lower-hybrid  frequency  can 
interact  simultaneously  with  both  the  electrons  and  the  ions,  and  thus 
can  transfer  energy  from  the  ions  to  the  electrons.  The  fluctuating 
electric  field  of  the  lower-hybrid  waves  is  primarily  electrostatic  and 
propagates  nearly  perpendicularly  to  the  auroral  magnetic  field.  It  is 
assumed  that  the  instability  saturates  rapidly.  In  this  paper  we  are 
not  concerned  with  the  details  of  the  saturation  mechanism  but  rather  with 
the  acceleration  of  the  electrons  by  plasma  turbulence  that  has  already 
reached  the  saturation  level.  Recently,  some  laboratory  experiments 
have  demonstrated  that  the  combined  effects  of  a  parallel  electric 
field  and  plasma  turbulence  is  particularly  effective  in  producing 
energetic  electrons  (Yamamoto  et  al.,  1980).  This  phenomenon  has  been 
observed  in  a  tokamak,  and  we  expect  that  the  physical  mechanism  should 
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be  operative  in  the  auroral  region  also.  It  is  based  on  this  notion 

that  we  are  motivated  to  proceed  with  the  present  study.  The  approach 

consists  of  two  parts;  one  is  the  study  of  a  test  particle  problem,  and 

the  other  is  the  investigation  of  the  time  evolution  of  the  distribution 

function.  The  latter  is  based  on  a  quasilinear  kinetic  equation. 

Our  discussion  starts  with  the  quasilinear  kinetic  equation  for  a 

magnetized  plasma  (Davidson,  1972;  Kennel  and- Engelmann,  1966),  including 

the  quasistatic  parallel  electric  field.  Since  the  plasma  turbulence  has 

a  characteristic  frequency,  U)  £5  ,  and  since  we  only  consider 

long  wavelengths,  kt  v  ,  <  ft  (v  .  is  the  electron  thermal  velocity), 

-  eth  e  eth 

the  most  important  diffusion  process  is  associated  with  the  n  =  0  term. 
In  this  case,  the  quasilinear  equation  reduces  to 


3F  ,  |  3F  *  3F 

—  e _ Li  1  E.i  _ —  =  —  d  ■  D  i(v  )  _ — 

<5t  m  H  3v  3v  N  M  3v. 

e  |i  11  il 


(1) 


where 


D||(v,i)  = 


o  2  2 

8tt  e 


m 


d3kek45(“-v») 

k 


(2) 


In  Eqs .  (1)  and  (2)  F^  is  the  electron  distribution  function,  and  the 
spectral  energy  density  ££•  =  <|e^|^>/8tt  ,  where  E£  is  the  fluctuating 

field . 

To  facilitate  our  discussion  we  adopt  a  coordinate  system  in  which 
the  z  axis  is  parallel  to  the  ambient  magnetic  field  and  the  turbulence 
is  excited  by  a  cross-field  drift  vd  parallel  to  the  y  axis.  For  a 
given  unstable  mode,  the  frequency  0)  and  the  wavevector  k  are  related  by 


a)  -  k*v ,  =  k  v ,  . 
d  yd 


(3) 
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We  further  assume  spectral  energy  density  has  the  form 


££  =  iKkjHCkjj) 


(4; 


where  4> (kjj)  is  the  power-law  spectrum  discussed  recently  by  Wu,  Gaffey 
and  Liberman  (1981) , 
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The  spectral  index  ct  is  considered  to  be  a  parameter,  and  d)  is 

00  * 

malized  such  that  /  dk  *  (ku)  -  1  ,  for  a  >  1  .  The  saturation  wave 

—00 

energy 


w  = 

s 


d  k  e?  = 
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d  k,  iKkt) 


(6) 


is  discussed  in  Section  IIC.  In  the  following  we  further  assume  that 


9  to  2  2  2 

<<  V,<  and  ky  >>  K  >  kt  * 


The  diffusion  coefficient  (2)  can  be  simplified  using  Eqs.(4),  (5) 
and  (6) . 


Dil(v(j)  =  D 
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The  electron  plasma  frequency  0)^  =  /m^  and  =  (2T  /m  )** 

is  the  electron  thermal  velocity. 

It  is  convenient  to  introduce  the  dimensionless  parameter, 


E  = 


e|Ellvthe 

m  D 
e 


(9) 


and  the  dimensionless  variables, 


u  =  7T~  ’  do) 

the 

t 

- - 1 —  f  dt*D(t»>  .  (U) 

Vthe  0 


If  vd  and  Ws  can  be  taken  be  approximately  time-independent,  Eq. 
(ll)  simplifies  to 


T 


(Ha) 


In  terms  of  these  variables,  the  diffusion  equation  (^  can  be  written  as 


9_ 

3u 
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E 


(12) 
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C .  The  Physical  Model  of  the  Acceleration  Region 

Since  the  present  study  is  primarily  motivated  by  the  observations 
of  precipitating  energetic  auroral  electrons,  we  now  discuss  the  model  of 
the  acceleration  region,  which  was  mentioned  previously  in  Section  IIB. 

In  this  region,  where  the  quasi-static  parallel  electric  fields  occur, 
the  Hawkeye  and  IMP  6  satellites  have  measured  large  amplitude  electro¬ 
static  plasma  turbulence  at  the  lower-hybrid  and  ion  cyclotron  frequencies 
(Gurnett  and  Frank,  1972,  1977).  We  expect  that  the  turbulence  at  the 
lower-hybrid  frequency  is  a  manifestation  of  the  lower-hybrid-drift  in¬ 
stability,  which  is  excited  by  a  cross-field  current  due  to  a  relative 
electron-ion  drift  (Hsia  et  al.,  1979;  and  references  therein).  The 
sources  of  free  energy  available  in  the  acceleration  region  to  drive  the 
lower-hybrid-drift  instability  are  the  E  x  B  drift  and  the  plasma  density 
inhomogeneity.  In  this  region  there  is  a  rather  strong  electric  field 
perpendicular  to  the  ambient  magnetic  field,  Et  72  400  mV/m.  This  large 
perpendicular  electric  field  can  generate  quite  a  large  relative  cross¬ 
field  drift  between  the  electrons  and  the  ions.  In  the  large  electric 
field  region,  the  low  energy  particles  are  expelled  and  there  is  a  strong 
density  gradient  at  the  boundary  of  this  region.  These  ideas  are  sup¬ 
ported  by  the  observations  of  Whalen  and  Daley  (1979)  that  the  most  in¬ 
tense  energetic  electron  precipitation  occurs  at  the  edge  of  the  auroral 
arcs.  The  lower-hybrid  waves,  which  are  generated  at  the  edge  can  pro¬ 
pagate  into  the  center  of  the  regions,  and  thus  the  waves  are  able  to  ac¬ 
celerate  electrons  in  regions  of  large  spatial  extent. 


The  cross— field  electron  drift  velocity 
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and  the  cross-field  ion  drift  velocity  is 
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where  v_  is  the  E  x  B  drift  velocity,  v^ 

E  De 

electron  and  ion  diamagnetic  drift  velocities, 

the  electron  and  ion  temperatures  in  the  acceleration  region  anc}  the 

electron  and  ion  cyclotron  frequencies  are  Q  =  ' e’B^/m  c  and  Q  =  !e  R 

e  1  1  0  e  i  1  1  0 

/ELc  ,  The  ions  are  not  able  to  t  x  t  drift  because  the  ion  gyroradius  is 
comparable  to  the  length  scale  of  the  acceleration  region,  i.  e.-> 

r.  =  v  / Q .  =  (3£n  n/3x)  ^  . 

x  Di  l  n 

Recent  research  results  indicate  that  ion  diamagnetic  drift  energy  provides 
the  free  energy  reservoir  which  drives  the  lower-hybrid-drift  instability 
(Hsia  et  al.,  1979),  The  saturation  electrostatic  turbulent  energy  den¬ 
sity,  given  in  Eq.  (8),  may  be  expressed  as. 


and 

VDi 
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where  y  is  the  energy  conversion  efficiency  coefficient. 
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III.  test  particle  problem 

The  time  evolution  of  the  velocity  and  of  the  energy  may  be  obtained 
by  taking  moments  of  the  diffusion  equation  (12)  and  integrating  by  parts. 
For  a  test  particle  we  take  the  distribution  function  to  be 

Fe(u,x)  =  6[u  -  u(t) ]  ,  (16) 

and  obtain 

Cot-3)  *  u  |  sgnu-E  ,  (17) 

and 

2 Cct-2)  |u|a~3  -  2Eu  .  (18) 

2 

In  Eq.  (18)  the  dimensionless  energy  £  =  u  .  This  expression  is  similar 
to  a  criterion  for  statistical  acceleration  that  has  been  used  frequently 
(Melrose  1968,  1969)  but  includes  the  effect  of  the  parallel  electric  field. 
We  caution  the  reader  that  some  care  must  be  taken  when  using  the  expres¬ 
sions  (17)  and  (18).  As  was  pointed  out  in  a  recent  paper  (Wu,  Gaffey  and 
Liberman,  1981) ,  these  equations  are  valid  only  when  the  terms  that  were 
integrated  by  parts  are  well  behaved  functions  that  are  integrable  over 
the  entire  range.  For  the  problem  considered  in  this  paper,  a  >2  is 
required. 

In  the  following  subsections  we  give  solutions  to  Eqs.  (17)  and 
(18)  for  various  cases  so  that  a  comparison  may  be  made. 

A.  Electric  Field  Only 

Under  the  influence  of  the  electric  field  alone,  without  diffusion, 
the  velocity  is 


u(T)  =  -  ET  +  u(0)  , 


(19) 
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where  the  initial  velocity  at  time  x  =  0  is  u(0)  ,  and  the  energy  is 

£ (x)  =  E2T2  -  2u(0) Ex  +  £(0)  ,  (20) 

2 

where  the  initial  energy  at  time  x  =  0  is  e(0)  =  u  (0)  . 

B.  a  =  3  case 

The  velocity  is 

u(x)  =  -  EX  +  u(0)  .  (21) 

If  the  electric  field  is  zero,  the  energy  is 

£(x)  =  2x  +  e(0)  ,  for  E  =  0  .  (22) 

If  the  electric  field  is  nonzero,  Eq.  (18)  may  be  integrated  to  obtain 

/eCO  -  -  Et  sgn  u  -  +  ^(0)  ,  £or  E  +  0  ,  (23) 

Ve/£(0)  sgn  u-1' 

where  sgn  u  is  the  sign  of  the  velocity.  For  E  small,  the  logarithm 
may  be  expanded,  and  to  second  order  Eq .  (22)  is  recovered.  For  E  large, 
the  logarithm  term  is  small,  and  we  iterate  to  obtain 

e(T)  «E2t2  +  2T  in  f— "  \  ,  (23a) 

V/e(0T/ 

where  sgn  u  =  -  1  has  been  used. 

C.  a  =  4  case 

The  velocity  is 

u(x)  =  -  (E-sgn  u  )t  +  u(0)  ,  (24a) 


and  the  energy  is 
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e(t)  =  (2-E  sgn  u)2t2  +  2v/e(0T(2-E  sgn  u)  x  +  e(0 ) 


(25a) 


For  E  >  1  ,  all  particles  are  eventually  accelerated  to  negative  velo¬ 
cities, 


u(T)  =  -  (E+1)T  +  u(0)  , 


(24b) 


and  the  energy  becomes 


2  2 


e(T)  »  (E+2)  x^  +  2/iToT(E+2)x  +  e(0)  . 

D.  Averaged  moments 

For  many  purposes  it  is  useful  to  have  the  moments  averaged  over 
the  distribution  of  initial  values,  which  we  take  to  be  a  Maxwellian. 
We  define  the  average  value  of  a  quantity  A  to  be 


(25b) 


<A>  =  —  [  Ae~U  (0)du(0) 


(26) 


For  the  case  with  only  the  electric  field,  the  average  velocity  is 


<u(t)>  =  -  EX  . 


(27) 


The  average  energy  is 


<£(T)>  =  E2X2  +  i 


(28) 


The  mean  velocity  squared  is 


<(Au)2>  =  <u2>  -  <u>2  =  j 


For  the  case  a  =  3  ,  the  average  velocity  is 


(29) 


<u(x)>  =  -  EX 


(30) 


If  the  electric  is  zero,  the  average  energy  from  Eq.  (22)  is 


<e(x)>  =  2x  + 


(31) 


and  the  mean  velocity  squared  is 


<(Au)2>  =  2t  +  y  . 


(32) 


If  the  electric  field  is  large,  the  average  energy  from  Eq.  (23a)  is 


<£(T)>  ar  e2T2  +  2x£n(Ex)  +  (2£n2  +  c) 


(33) 


where  Euler's  constant  c  is  defined  as 


_  lim 
c  =  n-x=° 


n  1 

(  )  -  -  Jin  n)  =  0.577... 

,  m 


(34) 


m=l 


(Gradsteyn  and  Ryzhik,  1965) .  The  mean  velocity  squared  is 


<(Au)  >  =  2x£n(Ex)  +  x(2£n2  +  c)  . 


(35) 


For  the  case  a  =  4  ,  for  E  <  1  ,  the  average  velocity  is 


<u(x)>  =  -  Ex  , 


(36a) 


the  average  energy  from  Eq.  (25a)  is 


<£(x)>  =  (E2+4)x  +  —  x  + 

/f  2 


(37a) 


and  the  mean  velocity  squared  is 


<  (Au)  2>  =  4x2  +  —  x  +  1  . 

2 


(38a) 


For  E  >  1  ,  the  average  velocity  is 


<u(x) >  =  -  (E+1)X  , 


(36b) 
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the  average  energy  from  Eq.  (25b)  is 

<£(!)>  =  (E+2)2t2  +  —  (E+2)x  +  j  , 

/rr 


(37b) 


and  the  mean  velocity  squared  is 


<(Au) 2>  =  (2E+3)t2  +  — (E+2)x  +  i  . 

/rr  2 


(38b) 
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IV.  KINETIC  THEORY 


In  this  section  we  obtain  solutions  to  the  diffusion  equation 
(12)  for  the  electron  Green's  function  and  distribution  function,  and  for 
moments  of  the  Green's  function  and  distribution  function.  Before  pro¬ 
ceeding  we  mention  three  points. 

First,  in  order  to  conserve  particle  number  density  the  functions 

1  I “-3  8Fe  „ 

lul  and  Fe  must  be  continuous  (39) 

and  vanish  at  the  endpoints.  This  can  be  seen  by  integrating  the 
diffusion  equation  ( 1 2 )  over  the  range  of  velocity. 

Second,  the  wave-particle  resonance,  k,,v.(  =  to  ,  is  not  expected 
to  be  satisfied  for  very  small  velocities  because  lower— hybrid  waves 
usually  have  k„  <  k(me/lO  2.  However,  since  the  diffusion  coefficient 
is  proportional  to  the  spectral  energy  density  ej  which  varies  as  k~a  , 
the  resonant  electrons  with  small  parallel  velocities  will  diffuse  very 
little. 

Third,  in  Section  III  it  was  pointed  out  that  the  spectral  index 
a  must  be  >  2  .  Furthermore,  a  must  be  <  5  so  that  the  diffusion 
process  does  not  "runaway"  to  infinite  velocities.  In  the  following 
we  shall  find  solutions  for  a  =  3  and  4  so  that  a  comparison  can  be 
made . 

We  solve  the  diffusion  equation  (12)  by  using  the  method  of  Laplace 
transforms.  The  Laplace  transform  of  (12)  is 


where  the  initial  value  fg(u)  —  ^  Q  (u50) .  We  first  construct  a 
Green’s  function  g(u,uf;s)  which  statisfies 


(41) 


(42) 


Then 


G(u, u 1 ;t)  ,  the  inverse  Laplace  transform  of  g(u,uT;s), 


is  calculated,  and  the  time  evolution  of  the  distribution  function  is 
given  by 


CO 


(43) 


The  Green's  function  G(u,u';T)  is  somewhat  analogous  to  the  test  particle 
distribution  (16),  and  the  velocity  and  energy  moments  of  G(u,u';T) 
are  calculated  and  compared  to  the  velocity  and  energy  moments  calculated 
for  the  test  particle  problem  in  section  III. 

In  the  next  subsections  we  carry  out  the  solution  procedure  for 
various  cases,  as  we  did  in  section  III,  so  that  a  comparison  may  be  made. 
A  Electric  Field  Only 

Under  the  influence  of  the  electric  field  alone,  without  diffusion, 
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the  Vlasov  equation  is 


3F  9F 

ir-E3Tr-°  • 


(44) 


which  is  satified  by  any  function  F  (u,x)  =  F  (u  +  Ex).  The  Green’s 

e  e 

function  is 


G(u,uT;x)  =  6(u-uT  4*  Ex) 


(45) 


For  an  initial  Maxwellian, 


F  (u,0)  =  —  e~U 
6 


(46) 


the  distribution  function  is 


^  1  -(u+ET)‘ 

F  (u,x)  =  —  e 

/if 


(47) 


The  velocity  moment  of  the  Green's  function  (45)  is 


u(t)  =  -  Et  +  u' 


(48) 


and  the  energy  moment  of  Eq.  (45)  is 


2  2 

£(r)  =  E  T  -  2u'Et  +  e 


(49) 


The  velocity  moment  of  the  distribution  function  (47)  is 


<u(T)>  =  -  ET  , 


(50) 


the  energy  moment  of  (47)  is 


<£(T)>  =  E2t2  +  j  , 


(51) 
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and  the  mean  velocity  squared  is 


< (Au) 2>  =  i 


(52) 


B .  Solution  for  the  case  cl  =  3 

For  the  case  of  the  spectral  index  &  =  3,  the  diffusion  equation 
(42)  is  particularly  simple. 


Su 


sg  =  -  6 (u-uT) 


(53) 


The  solutions  of  the  associated  homgeneous  differential  equation  are 


1  + 

g+(u)  =  e  , 

(54) 

where 

O  1^ 

E  +  (E  +  4s) 2 

Y±  =  2 

(55) 

We  note  that  y  <  0  <  y+  for  E  >  0 

Because  g  (u  -  <»)  =0  and  g+(u  -*■<»)  =  0  ,  we  look  for  a  solution 
fo  Eq.  (53)  of  the  form. 


-(Y  u+y  u') 
g(u,uT;S)  =Ae  “  + 


tt/  ,  ,  ~(Y  u+y  u') 

U(u'-u)  +  Ae  + 


U(u-u ' ) 


(56) 


where  A  is  a  function  of  u'  and  s  ,  U  is  a  Heaviside  step  function, 
and  we  have  required  that  g  be  continuous  at  u  =  u'  . 

The  coefficient  A  is  determined  by  substituting  Eq.  (56)  into 
Eq.  (53).  The  desired  solution  has 


A 


(Y++Y_)uT 

e 


J 


(57) 
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and  thus. 


g(u,u';s)  = 


1  -^E(u-u')  -/s  +  isE2|u-u' 


l/s  4-  JjE2 


(58) 


The  inverse  Laplace  transform  is  well  known  (Abramowitz  and  Stegun 
1964,  Erdelyi  1954).  Thus, 


,  .  1  -(u-u'+Et)  / 4t 

G(u,u  ,t)  =  — —  e 


(59) 


2v^T 


The  distribution  function  is  obtained  from  Eq.  (43),  using 
function  (59)  and  the  initial  condition  (46). 


F  (u,t)  =  - i -  e  -(u+Et).^(4x4-1) 

6  /  (4x4-1)  tt 


If  the  electric  field  is  zero,  then  the  Green's  function 


is 


2 

G(u,u'  ;  x)  =  — —  e  /4t  » 

2/Jtt 


the  Green's 


(60) 


(59a) 


which  is  the  well  known  solution  to  the  heat  equation  (Morse  and  Feshbach, 
1953),  and  the  distribution  function  is 

F-(u,t)  -  — ^ -  e-2/(4 T-U)  _  (60a) 

A  (4x4-1) 

which  was  given  by  Wu,  Gaffey  and  Liberman  (1981). Note  that  the  results  in 
this  paper  differ  by  a  normalization  factor  from  those  of  Wu,  Gaffey  and 
Liberman. 

The  velocity  moment  of  the  Green's  function  (59)  is 


u(x)  =  -  EX  4-  u’ 

and  the  energy  moment  of  Eq.  (59)  is 


(61) 
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£(t)  =  E2T2  +  2 (I-u'E)t  +  e' . 


(62) 


The  velocity  moment  of  the  distribution  function  (60)  is 


<u(t)>  =  -  Ex 


(63) 


the  energy  moment  of  Eq.  (60)  is 


<£(x)>  =  E2X2  +,  2x  +  |  , 


(64) 


and  the  mean  velocity  squared  is 


< (Au) 2>  =  2x  +  y  .  (65) 

C .  Solution  fox  the  case  q  =  4 

For  the  case  of  the  spectral  index  a  =  4  ,  the  diffusion  equation 
(42.)  is 


h  ,u|  If +  E 15  ■ sg  =  -  6(u-u,) 


(66) 


Changing  variables  to  z  =  2/[uJ"  yields, 


3  g+  (1+2E)  3g 


3-2 


3z 


-  -  sg  ±  =  -  —  6(z-z ' )  , 


(67) 


where  the  j-  signs  are  for  u  0  .  Without  loss  of  generality  we 
have  taken  E  >.  0  . 

C*1  Solution  for  u  >  0 

The  Green Ts  function  for  u  >  0  is  derived  in  Appendix  A.l.  From  Eq - 

(A5) 


G+(u,u>;x)  =  ^(^) 


,  E/2  - 

e 


(u+u') 


E  x  ' 


(68) 
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This  expression  is  rather  complicated.  In  many  situations,  the  long  time 
asymptotic  behavior  is  of  more  interest.  Using  the  small  argument  ex¬ 
pansion  of  the  Bessel  function, 


G+(u,u' ;t) 


E  _  (u+U1) 

u'  T 

—  —  0 

2r(l+E)x(1+E) 


(6?) 


The  distribution  function  is  obtained  from  Eq.  (43),  using  Eqs .  (69) 
and  (46) .  Using  the  small  argument  expansion  of  the  parabolic  cylinder 
function  and  the  duplication  formula  for  the  gamma  function  (Abramowitz 
and  Stegun  1964),  r(l+E)  =  Tr" V  r (!+£) r ( 2+E}  >  we  find 


F+(u,x) 


e 


ti 

T 


2  4-  E 

re-^)(2x) 


(l+E) 


(70) 


We  note  that  Eq .  (70)  satisfies  (12)  exactly  for  a  =  4  . 

C .  2 .  Solution  for  u  <  0 

The  Green fs  function  for  u  <  0  is  derived  in  Appendix  A. 2. 
From  Eqs.  (A13)  and  (A15) 


f  _  ( M-Hu?  1) 

T 


1 

u 

2  “  i— 

T 

uu T 

2t 

u' 

e 

I 

XE( 

T 

+ 


uu 


T 

f 


2F (E)F(1+E) 


dA 


[X(x-X)] 


(l+E) 


(71) 


0 
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The  long-time  asymptotic  behavior,  for  T  >  2/\  uu ' |  , 
is  found  by  using  the  small  argument  expansion  of  the  Bessel  function  in 
the  first  term. 


G  (u,u',x)  = 


2r(l+E)T(1+E) 


(|u|  +  lur  1) 

T 


+ 


UU 


»|E 


dX 


_  (M  +  inli) 

d-A  a  ' 


2r(E)r(1+E)  J  [A(t-A)](1+E>' 


(72) 


The  distribution  function  is  obtained  from  Eq.  (43)  using  Eqs. 

(72)  and  (46).  x 

.  i  E  -  M  T 


F  <»■*>  *  — e  T  e4T  “fc  + 


2r(i+E)x 
i  E 


+ 


u 


^2(1^E)r(E) 


dA 


.M 

t-A  8A~ 


[A(t-A) ] 


(1+E) 


(73) 


where  erfc  is  a  complementary  error  function  and  D  .  is  a  parabolic 

-(.1+E) 

cylinder  function. 

For  T  >  %  ,  the  first  term  can  be  expanded  to  give 


F  (u,x)  = 


_  M 

T 


2r(l+E)T 


(1+E) 


+ 


v/TT2(^)r(E) 


dA 


_  M 

t-A  8A 


0 


[A(t-A) ] 


(1+E) 


M 
2 

e  D 


-(1+E) 


(— ) 
v/2A 


(74) 


We  note  that  Eq.  (74)  satisfies  Eq.  (12)  exactly  for  a  =  4  . 

If  the  electric  field  is  zero,  then  the  Green's  function  is 


du  +  ui  h 


G(u,u '  ;x)  =  ^  e 


IQ(-  v'luu'  f) 


(75) 
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Thfi  long- time  asymptotic  behavior  is  found  by  talcing  the  small  argument 
expansion  of  the  Bessel  function. 


(lul+lu'l) 


G(u,u' ;x) ~  2r  e 


(76) 


The  long-time  asymptotic  form  for  the  distribution  function  is 


1  _  M 

Fe(u,T)  =  ^  e  x 


(77) 


These  results  were  obtained  recently  by  Wu,  Gaffey  and  Liberman  (1981) . 
Note  that  the  functions  in  this  paper  differ  by  a  normalization  factor 
from  those  in  the  previous  paper. 

The  velocity  moment  of  the  exact  Green's  functions  (71)  and  (68)  is 


u(T)  =  -  i  T{  (E+l)+ 


+  SW)  [r(E’  T  }  -^1)] 


—Ljy  (JulL)  *(E-i,Efl;  -  ^ 


-)} 


(78) 


and  the  energy  moment  is 


e(T)  -  T2  ij  (E+2)  (E+l)  +  (E+2)  -^i  +«1 1  +  .(E+2)(E+l)[r  JVL  _ 
*■  t  .  l  2i  (E)  ’  x 


2x 


-  2 1L1  r<E-i,LLi>+  it!  r(E-2,lti)) $(e-2,e+i;-  1h5L)} 


(79) 


where  T(a)  is  the  gamma  function,  T(a,x)  is  the  incomplete  gamma  func¬ 
tion  and  $(a,b ;x)  is  the  confluent  hypergeometric  function  (Abramowitz 
and  Stegun,  1964;  Gradsteyn  and  Ryzhik,  1965).  The  first  several  terms 
of  the  long-time  expansions  of  Eqs.  (78)  and  (79)  are 
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u(T)  =  -  T  [ (E+l)  - 


i  u 


(E-l) 


+  .  .  •  ] 


for  E  >  1 


(78a) 


and 


e(T)  ~t  [(E+2)  (E+l)  - 


2 (E+2) 
(E-l) 


u 


(E2+2) 


,2 


(E-2) (E-l) 


+...]» for  E  >  2 


(79a) 


where  E  >  1  has  been  assumed  in  Eq.  (78a)  and  E  >  2  in  Eq.  (79a) 

I  I  E 

to  neglect  terms  of  order  (|u,|/t) 

The  dominant  terms  of  the  velocity  and  energy  moments  of  the  approxi¬ 
mate  Green Ts  functions  (72)  and  (69)  are 


u(t)  —  -  (E+1)T 

(78b) 

and 

e(T)  —  (E+2) (E+l) t2  , 

(79b) 

which  agree  with  the  expansions  (78a)  and  (79a)  to  leading  order. 

The  first  few  terms  of  the  velocity,  energy  and  mean  velocity  squared 
moments  of  the  exact  distribution  functions  obtained  by  averaging  Eqs. 
(78a)  and  (79a)  over  an  initial  Maxwellian  distribution  are 


<u(x)>  =  -  (E+l) x  + 


/F*  (E-l) 


+  .  .  -  , 


for  E  >  1  (80) 


<£(T)>  -  (E+2) (E+l) T2  -  -(E+2)T  +  ^  +...,  for  E  >  2, 

/F  (E-l)  2 (E-2) (E-l) 


(81) 


and  <  (Au)  2>  =  (E+l)t2 - — -  +  E 2+2 


/fr(E-l)  2 (E-2)  (E-l)  2 


+  .  . ,f  or  E  >  2.  (82) 


The  dominant  terms  of  the  velocity,  energy  and  mean  velocity  squared 
moments  of  time  asymptotic  distribution  functions  (74)  and  (70)  are 


<u(x)>  «  -  (E+l)x  , 


(80a) 
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<e(x)>  =  (E+2)(E+1)t2 

9 

(31a) 

and 

<(Au)2>  -  (E+Dt2 

(82a) 

which  agree 

with  Eqs .  (80),  (81)  and  (82)  to 

leading  order. 

If  the 

electric  field  is  zero,  the  moments  (78)  and  (79)  may  be 

easily  evaluated  to  obtain 

u(t)  =  (T  +  | u ' | )sgn  u  , 

for  E  =  0  , 

(78c) 

and 

£(t)=2t2  +  4 |u' | T  +  u'2  , 

for  E  =  0  . 

(79c) 

The  moments 

averaged  over  a  Maxwellian  are 

<u(t)>  =  0 

,  for  E  =  0 

(80b) 

<£  (T)>  =  2T2  +  —  +  \ 

/n  2 

,  for  E  =  0  , 

(81b) 

2  2  4  1 

<(Au)Z>  =  2x  +  —  +  ^ 

/? f  2 

,  for  E  =  0  . 

(82b) 
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V.  DISCUSSION  AND  SUMMARY 


A  summary  of  the  results  of  the  test  particle  moments  calculated  in 

Section  HI  is  given  in  Table  1.  A  summary  of  the  kinetic  theory  moments 
calculated  in  Section  IV  lis  given  in  Table  2.  The  averaged  test  particle 
moments  are  given  in  Table  3  and  the  averaged  kinetic  theory  moments 
are  given  in  Table  4. 

We  remark  that  the  energy  £(x)  varies  linearly  with  time  for  a  dif¬ 
fusion  process,  whereas  e(x)  varies  quadratically  with  time  for  an  ac¬ 
celeration  process. 

When  the  electric  field  is  zero,  under  the  influence  of  quasilinear 

effects  alone,  for  the  case  a  ~  3  ,  the  energy  e(x)  ^  T  ,  which  is 

2 

characteristic  of  a  diffusion  process,  and  for  a  -  4  ,  £(t)  'V  t  because 

the  quasilinear  kinetic  equation  (12)  contains  an  acceleration  term  as  well 
as  a  diffusion  term  for  this  case.  Thus,  the  larger  the  spectral  index 
a  (the  narrower  the  spectrum)  the  more  effective  the  energization  process. 
We  further  note  that  in  the  absence  of  an  electric  field,  the  quasilinear 
effects  increase  the  electron  energy  symmetrically  for  electrons  with  po¬ 
sitive  and  negative  velocities,  since  in  this  case  Eq.  (12)  depends  only 
on  the  absolute  value  of  the  velocity.  This  results  is  in  agreement  with 
satellite  observations  of  Sharp,  Shelley,  Johnson  and  Ghielmetti  (1980). 

When  the  electric  field  is  large,  E  >  1  ,  the  long  time  asymptotic 

2  2 

velocity  u(t)  £2  -  Ex  and  the  energy  e(x)  S3  E  X  in  all  cases.  This  is 
characteristic  of  an  acceleration  process.  For  a  =  3  ,  the  fact  that  the 
test  particle  moments  and  the  kinetic  theory  moments  are  approximately 

comparable  if  £n(Ex)  ^  0(1)  may  be  attributed  to  the  diffusion  being 
constant,  independent  of  the  velocity.  For  a  =  4  ,  the  more  rapid 
increase  with  time  of  the  test  particle  moments  compared  to  the  kinetic 
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theory  moments  may  be  attributed  to  the  fact  that  the  test  particle  dis¬ 
tribution  function  (16)  was  assumed  to  remain  a  delta  function  whereas 
the  Green Ts  function  spreads  with  time*  The  Green’s  functions  for  a  =  3 
and  4  both  have  denominators  in  the  exponentials  that  vary  linearly  with 

time,  which  gives  the  velocity  spread  which  increases  with  time.  However,  for 
the  case  of  electric  field  acceleration  alone,  the  Green’s  function 

G(u,uT;x)  =  6(u  -  uT  4*  Ex)  ,  because  the  diffusion  has  been  neglected. 

An  important  difference  between  the  various  cases  is  the  spread  in 

velocity.  For  electric  field  acceleration  only,  the  mean  velocity  squared 
2 

<(Au)  >  =  h  >  which  is  the  thermal  spread  of  the  initial  Maxwellian.  For 
2 

ct  =  3  ,  <(Au)  >  ^  t  ,  which  is  characteristic  of  diffusion  with  a  constant 
diffusion  coefficient.  For  a  =  4  ,  <(Au)2>  <\,  x2  ,  which  is  character¬ 

istic  of  an  acceleration  process. 

Another  feature  is  that  in  all  cases  the  leading  term  of  the  long 
time  asymptotic  moments  is  independent  of  the  initial  distribution.  This 
may  be  attributed  to  the  nature  of  the  diffusion  and  acceleration  processes. 

The  observations  of  precipitating  energetic  auroral  electrons  by 
rockets  at  low  altitudes  (Kaufmann  et  al.,  1976;  Lundin,  1976;  Bryant  et  al., 
1978)  and  by  satellites  at  higher  altitudes  (Whalen  and  Daly,  1979;  Sharp 
et  al . ,  1980)  have  several  features  that  cannot  be  accounted  for  simply 
by  the  model  of  laminar  acceleration  by  a  quasistatic  parallel  electric  field 
alone.  As  discussed  in  Section  IIA,  the  peak  value  of  the  potential  dif¬ 
ference  along  the  auroral  magnetic  field  lines  is  estimated  to  be  about 
4  or  5  kV,  and  therefore  can  provide  only  about  one  half  of  the  observed 
energy, ^  10  keV, of  the  precipitating  electrons.  The  additional  energy 
must  be  the  result  of  some  other  energization  process,  such  as  turbulent 
acceleration  by  electrostatic  plasma  waves.  If  the  electrons  fall  through 
a  potential  difference  $  ,  the  resulting  energy  gain  will  be  e$  ,  but 
the  velocity  spread  will  remain  the  same  as  that  before  the  acceleration. 
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Also,  in  this  case  the  electron  disrtibution  function  will  fall  off  rather 
abruptly  with  increasing  energy  above  e<J>  ,  whereas  the  observed  electron 
distribution  functions  tend  to  be  rather  broad  and  fall  off  more  slowly. 

Our  model,  which  includes  the  combined  effects  of  velocity  diffusion  and 
electric  field  acceleration  agrees  better  with  the  observed  electron  dis¬ 
tribution  functions . 

Because  the  electron  acceleration  occurs  near  the  open  auroral  magnetic 

3 

field  lines,  we  assume  a  dipole  magnetic  field  B  =  B  (1  +  h/R  )  ,  where 

0  E 

Bq  is  the  magnetic  field  strength  at  the  earth,  h  is  the  altitude  above 

the  earth,  and  is  the  radius  of  the  earth,  to  calculate  the  electron  and 

ion  cyclotron  frequencies  as  a  function  of  altitude.  The  numerical  values 

are  listed  in  Table  5.  The  values  for  the  electron  number  density  used  to 

calculate  the  plasma  frequencies,  are  also  given  in  Table  5.  The  values  taken 

for  the  physical  quantities  are  B„  =  0.40  Gauss,  T  =1  keV  and  T  =  200 

0  e  i 

9  ^  6 

eV,  which  give  v  .  =  1.88.10  cm/ sec  and  v,.^  %  v  .  .  =  9.79.10  cm/ sec. 

the  di  thi 

As  an  example,  we  take  the  intermediate  altitude  h  =  5000  km, 

Er|  =  1  mV/m  and  5  mV/m,  E,  =  400  mV/m,  y  «  10~3  ,  oo/k  v  ,  'V  0(1), 

=  1  keV  and  T\  =  200  eV .  In  Figure  1,  we  plot  the  electron  energy  as  a 
function  of  time  from  Eq.  (20)  for  acceleration  by  the  parallel  electric 
field  alone,  for  Ejj  =  1  mV/m  and  5  mV/m  and  an  initial  energy  of  1 
keV.  The  larger  electric  field,  of  course,  gives  rise  to  a  greater  acce¬ 
leration.  In  Figure  2,  the  electron  energy  is  plotted  as  a  function  of 
time  for  acceleration  by  the  plasma  turbulence  alone  from  Eq .  (22)  for 
a  =  3  and  from  Eq.  (25a)  for  a  =  4  for  an  initial  energy  of  1  keV. 

The  larger  the  spectral  index  a  (the  narrower  the  spectrum)  the  greater 
is  the  turbulent  acceleration  of  the  electrons.  Comparing  Figures  1  and  2, 
we  find  that  the  electric  field  acceleration  for  Ej|  =  1  mV/m  and  the 
turbulent  acceleration  for  a  =  4  are  very  nearly  comparable  at  an  alti- 
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tude  of  5000  km.  At  higher  altitudes  the  turbulent  acceleration  is  greater, 
whereas  the  electric  field  acceleration  is  independent  of  altitude.  In 
Figure  3,  the  electron  energy  is  plotted  as  a  function  of  time  for 
acceleration  under  the  combined  influence  of  plasma  turbulence  and  a 
parallel  electric  field  of  1  mV/v  and  5  mV/m  using  Eq.  (23)  for  a  =  3 
and  Eq.  (25a)  for  a  =  4  for  an  initial  energy  of  1  keV.  As  in  Figure 
2,  the  acceleration  is  greater  for  the  larger  value  of  a  .  Also,  the 
difference  between  a  =  4  and  a  =  3  is  more  pronounced  for  the  smaller 
value  of  the  electric  field.  Of  course,  when  both  acceleration  mechanisms 
are  operable,  the  acceleration  is  greater  than  when  only  one  acceleration 
mechanism  is  operable. 

In  this  paper  we  have  investigated  the  acceleration  of  auroral 
electron  by  parallel  electric  fields  and  electrostatic  plasma  turbulence 
near  the  lower-hybrid  frequency.  We  find  that  the  combined  effects 
of  a  parallel  electric  field  and  electrostatic  plasma  turbulence  is 
especially  effective  in  producing  energetic  electrons.  In  particular, 

the  present  study  has  led  to  the  following  major  conclusions: 

1.  A  larger  electric  field  produces  greater  acceleration;  however, 
acceleration  by  an  electric  field  alone  does  not  increase  the  velocity 
spread  of  the  electron  distribution  function. 

2.  A  larger  spectral  index,  i.e.  a  narrower  spectrum  of  waves,  pro¬ 
duces  greater  turbulent  acceleration.  Turbulent  acceleration  increases 
both  the  velocity  spread  and  the  average  energy  of  the  electron  distri¬ 
bution  function. 

3.  The  turbulent  acceleration  mechanism  is  more  effective  at  higher 
altitudes,  whereas  the  electric  field  acceleration  is  independent  of 


altitude . 


4. 


For  some  values  of  the  parameters  which  are  typical  of  the  auroral 
zones  the  turbulent  acceleration  and  the  electric  field  acceleration 
are  comparable. 

^  When  both  acceleration  mechanisms  are  operable  the  acceleration  is 
greater  than  when  only  one  acceleration  mechanism  is  operable. 
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APPENDIX 


In  this  appendix  we  derive  the  Green's  function  used  in  Section  IVC. 


A. 1 .  Derivation  of  the  Green's  Function  for  u  >  0 

The  homogeneous  differential  equation  associated  with  Eq.  (67)  for 
u  >  0  is 


9z 


+ 


(1+2E)  9g+ 
Z  9z 


-  sg+  =  0 


(Al) 


The  solutions  are 


g+(z) 


z  I+gCzv^)  and  z 


(A2) 


where 
bounded  as 
more, 


and  K£  are  modified  Bessel  functions.  Only  z'^Cz/s)  is 

-*E 

z  -*■  00  and  only  z  I  (zi/s)  is  bounded  as  z  0  .  Further- 


I  I  +  jL  + 

|U!  9 1 u|  ~  2  2  sT"  =  0  at  u  -  0  and  oo  ,  so  that  Eq.  (39)  is 
satisfied. 

We  look  for  a  solution  of  the  inhomogeneous  differential  equation 
(67)  of  the  form 


g_l_(z ,  z  ' ;  s)  = 


C+(zz ') 


-E 


[KE(z'y7)IE(zys')U(z'-z)  +1  (z'*^s)K  (z^)U(z-z')  ].  (A3) 


where  the  coefficient  C+  is  a  function  of  z'  ,  U  is  a  Heaviside  step 

function,  and  we  have  reauired  that  g+  be  continuous  at  z  =  z'  . 

The  coefficient  C+  is  determined  by  substituting  Eq.  (A3)  into  Eq . 

(67).  The  particular  solution  of  Eq.(70),  with  C  =  z'2E  ,  is 

4*  9 
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.  E 


g+(z,z';s)  =  (— )  [KE(z'/s)IE(z/s)U(z,-z)+IE(z'ys)KE(zv/s)U(z-z')  ]  .  (A4) 


The  inverse  Laplace  transform  of  g  is  (Erdelyi  1954) 


E  fz2+zr2v 

G+(I>*'!T)  -  4T  \,*z\ 


o 

Returning  to  the  original  velocity  variable  u  =  hz  ,  gives  Eq.  (68). 

A. 2  Derivation  of  the  Green's  Function  for  u  <  0 

The  homogeneous  equation  associated  with  Eq.  (67)  for  u  <  0  is 


(A5) 


82g-  ,  (1  -  2E)  3s- 

a  2  z 

a  Z 


3z 


-  Sg  =  0  . 


(A6) 


The  solutions  are 


g_(z)  -  zEI±E(2/i)  and  z^Czv's)  < 


(A7) 


E  __  TJ» 

Only  z  Kg ( z v s )  is  bounded  as  z  -»■  °°  .  Both  z^I  (z/s)  are  well  be¬ 
haved  as  z  -*■  0 .  Furthermore , 


,u|  Sg-  =  i  K. 

U  3R  2  2  8z  0  at  z  -  0  and  co  ,  so  that  Eq.  (39)  is  satisfied. 


We  look  for  a  solution  of  the  inhomogeneous 


equation  (67)  of  the  form 


g_ (z ,  z  '  ;s)  =  (zz')E{[B_Ie(z'v^)  +  C_I_E(z  ’  ^*)  ]Ke(z/F)U(z-z ’) 
+  KE(z'y^)[B_  IE(zvS>  +  C_I_E(z/i‘)]U(z'-z)}  . 


(A8) 
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We  will  need  two  conditions  to  determine  the  two  coefficients  B  and 
C  .  Continuity  at  z  =  z'  has  been  imposed. 


Requiring  that  g_(z,z';s) 


lim  ,  , 

z->0  S+\z>z  ;s)  gives 


C 


rq-E)  ,s,E 

rq+E)  V 


(A9) 


The  coefficient  B_  is  determined  by  substituting  Eq.  (A8)  into  Eq.  (67), 


B 


_,-2e  rq-E)  rs  E 

z  -  rq+iT  ( 4 } 


(A10) 


With  the  coefficients  given  in  Eqs.  (A10)  and  (A9) ,  the  Green's  function 
(A8)  becomes 


8_(Z,2';s)  -  (fr)  {IE(«’^  +  (5f^)E[I_E(z'^)-IE(2.^)]}KE(2^)U(2-2’ 

+  <7r)  +  fa+gj  : )[i_E(z^)-i 


(All) 

Expressing  I_g  -  I  in  terms  of  K^,  and  using  the  reflection  formula 
for  the  gamma  function  (Abramowitz  and  Stegun  1964),  T(l-E)  =  TT/r(E)sin  ttE  , 
we  separate  the  inhomogeneous  and  homogeneous  terms  in  Eq.  (All). 

g-<z’z';s>  '  ^ 

j_  2  szz ' .  E  _ 

r(E.»rq+E)  (~4_)  . 


Z  )  ]  + 


(A12) 


The  inverse  Laplace  transform  of  the  first  term  g,  (z,z';s)  is 


E  -  (z2jrz'2\ 

G1-(2’2,’t)  -  Jr  (fr)  e  4t  °  z  («\ 
T  2  EV2x  ; 


(A13) 


41 


The  inverse  Laplace  transform  of  the  second  term 


g2_(z,z ' ;s) 


is 


obtained  from  the  formula  (Erdelyi  1954) 


E  a 

L-'Kf)2  V2V=>1  -^(We’T 


(A14) 


Using  the  convolution  formula,  we  find 


g2_(z,z,;t)  = 


^T~) 


2E 


2r(E>ra+E) 


X 

f 


J 


dA 


-[ 


J  [X(t-A)](1+E) 


2  z.2 

4(x-A) +  4A  ] 


(A15) 


Returning  to  the  original  velocity  variable  |u| 


gives  Eq.  (71) 


TABLE  1 


Test  particle  moments  for  various  cases 


u(T)  =  -  Et  =  u(0) 

electric  field  only 

e(t)  =  E2t2  -  2u(0)Ex  +  e(O) 

u(x)  =  u(0) 

a  =  3 

£(x)  =  2x  +  £ (0) 

diffusion  only 

u(T)  =  x  sgn  u(t)  +  u(0) 

a  =  4 

£(T)  =  4t2  +  4t/£(0)  +  £ (0) 

diffusion  only 

u(t)  =  -  Et  +  u(0) 

a  =  3  ,  E  >  1 

0  0  T7T 

e(t)  «  ET  +  2x  £n(— — ) 

v'e(O) 

diffusion  and  electric  field 

u(T)  =  -  (E+1)t  +  u(0) 

a  =  4  ,  E  >  1 

£ (x)  =  (E+2)2t2  +  2/e(0T(E+2)t  +  £ (0) 

diffusion  and  electric  field 

Table  2 


Kinetic  theory  moments  for  various  cases 


u(t)  =  -  Et  +  u? 

e(t)  =  E2t2  -  2u'Ex  +  e* 

u(t)  =  u' 

£(t)  =  2t  +  £’ 

u(t)  =  (x  +  |u’|)  Sgn  u 
£ (t)  =  2t2  +  4 |u'  |t  +  £  * 

u(t)  =  -  Et  +  u' 

£(t)  =  E2t2  -  2(u'E-1)t  +  £  * 

u(x)  «  -  (E+1)T 
£ (x)  «  (E+2) (E+l)x2 


electric  field  only 

a  =  3 

diffusion  only 

a  =  4  , 

diffusion  only 

a  =  3 

diffusion  and  electric  field 

a  =  4,  E  >  2 

diffusion  and  electric  field 


TABLE  3 


Averaged  test  particle  moments  for  various  cases 


<u(x)>  =  -  Et 

<£(T)>  =  E2T2  +  j 

<(Au)2>  =  | 

electric  field  only 

<u(T)>  =  0 

a  =  3 

<e(T)>  =  2r  +  j 
<(Au)2>  =  2t  +  | 

diffusion  only 

<u(t)>  =  0 

<e(t)>  =  4t2  +  —  t  +  4 

/if  2 

<(Au)2>  =  4t2  +  —  t  +  -ir 

/if  2 

a  =  4  . 

diffusion  only 

<u(t)>  =  -  Et 

a  =  3  ,  E  >  1 

<£(!)>  «  E2T2  +  2t  £n(Ex) 

diffusion  and  electric  field 

<(Au)2>  «  2t  £n(Ex) 

<u(x)>  =  -  (E+1)T 

a= 4  ,  E  >  1 

<£(T)>  =  (E+2)2t2  +  —  (E+2)t  +  \ 

/n  2 

<(Au )2>  =  (2E+3)x2  +  —  (E+2)t  +  \ 

/if  2 

diffusion  and  electric  field 

TABLE  4 


Averaged  kinetic  theory  moments  for  various  cases 

<u(t)>  =  -  Et 


<e(x)>  =  E2t2  +  j 
<(Au)2>  =  | 

electric  field  only 

<u(x) >  =  0 

a  =  3 

<e(T)>  =  2x  +  | 

<(Au)2>  =  2t  +  j 

diffusion  only 

<u(t)>  =  0 

a  =  4 

<£(T)>  =  2x2  +  —  T  +  4 
/F  2 

<(Au)  2>  =  2t2  +  —  t  +  4 
/if  2 

diffusion  only 

<u(t)>  =  -  Et 

a  =  3 

<£ (x) >  =  E2t2  +  2t  +  -| 

diffusion  and  electric  field 

<(Au) 2>  =  2x  +  j 

<u(t)>  a  -  (E+l)x 

a  =  4  ,  E  >  2 

<£(T)>  «  (E+2) (E+1)T2 

diffusion  and  electric  field 

<(Au)2>«  (E+1)t2 

diffusion  and  electric  field 
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Plasma  parameters  in  the  acceleration  region 


FIGURE  CAPTIONS 


Figure  1  Electron  energy  (in  keV)  as  a  function  of  time  (in  sec) 
from  Eq.  (20)  for  acceleration  by  a  parallel  electric 
field  only,  for  =  1  mV/m  and  5  mV/m. 

Figure  2  Electron  energy  (in  keV)  as  a  function  of  time  (in  sec) 
from  Eq.  (22)  for  a  =  3  (dash  dot  curve)  and  Eq.  (25a) 
for  a  =  4  (solid  curve)  for  acceleration  by  plasma  tur¬ 
bulence  only. 


Figure  3  Electron  energy  (in  keV)  as  a  function  of  time  (in  sec) 
from  Eq.  (23)  for  a  =  3  (dash  dot  curve)  and  Eq.  (25a) 
for  a  =  4  (solid  curve)  for  acceleration  by  plasma  tur¬ 
bulence  and  a  parallel  electric  field,  for  E|j  =  1  mV/m 
and  5  mV/m. 
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Comments  on  "Energization  of  Ionospheric  Ions  by  Electrostatic  Hydrogen 

Cyclotron  Waves"  by  N.  Singh,  R.W.  Schunk  and  J.J.  Sojka,  Geophvs.  Res. 
Lett.;  8,  1249  (1981). 


J.D.  Gaffey,  Jr. 

Institute  for  Physical  Science  and  Technology 
University  of  Maryland 
College  Park,  Maryland 
20742 


In  a  recent  paper,  Singh  et  al.(1981)  make  a  comparison  between 
the  stochastic  acceleration  of  ionospheric  ions  by  electrostatic  hydro¬ 
gen  cyclotron  waves  observed  in  their  numerical  simulations  and  the 
predicted  acceleration  of  these  ions  based  on  the  theory  of  stochastic 
acceleration  (Papadopoulos  et  al.,  1980).  The  assertion,  made  by  Singh 
_et__al_,  that  the  heating  of  the  helium  ions  observed  in  their  simulations 
is  appreciably  larger  than  the  heating  predicted  by  the  theory  of  stochastic 
acceleration  (Papadopoulos  et  al.)  is  the  result  of  an  algebraic  error. 

In  fact,  upon  correcting  the  error,  we  find  the  agreement  between  the 
simulations  and  the  theory  to  be  quite  good. 


The  error  occurs  in  the  second  line  of  Eq.  (5)  of  Singh  et  al. ,  where 


the  dimensionless  wave  amplitude  is  defined  as 


a. 


■  ^1*0 
Mini 


(la) 


where  k±  is  the  perpendicular  wavenumber,  q  and  M^  are  the  ion 

charge  and  mass,  the  cyclotron  frequency  ft.  =  q  B/M  c  ,  and  $  i3  the 

I  i  i  o 

amplitude  of  the  electrostatic  potential  of  the  wave.  In  terms  of 
the  dimensionless  parameters  of  the  background  hydrogen  plasma, 


(lb) 


®H  = 

(2T^/M^)  .  Equation  (lb)  is  given  in  Papadopoulos  et  al.  on  the  line 

below  their  Eqs.  (3)  and  (9).  In  proceeding  from  Eq.  (la)  to  Eq.  (lb) 
Singh  et  al.  omitted  the  factor  1/2,  and  thus  obtained  a  value  of 
a  factor  of  2  too  large. 

The  values  of  the  dimensionless  parameters  of  the  background  hydrogen 
plasma  given  by  Papadopoulos  et  al.  on  the  line  below  their  Eq.  (8)  and 
by  Singh  et  al.  in  their  Eq.  (3)  are 


where  the  Larmor  radius  r^  =  and  the  thermal  velocity 


to  .  _ 

—  =  1.2 


Vh 


1.73  , 


H 


*  0.4 


(2) 


Substituting  the  numerical  values  of  Eq.  (2)  into  Eq.  (lb)  gives 


a. 

i 


n 

n  ,  ns  _  i 

0,6  a±  ~  2 


V, 


(3) 


where  the  frequency  ratio  .  The  numerical  values  of  W 

used  by  Singh  et  al.  in  their  numerical  simulations  is  a  consequence 
of  omitting  the  factor  of  1/2  in  Eq.  (lb)  and  Eq.  (3),  and  therefore 
is  twice  the  value  of  the  wave  amplitude  used  by  Papadopoulos  et  al . 

for  the  same  values  of  the  parameters  given  in  Eq.  (2). 

The  maximum  velocity  achieved  by  stochastic  acceleration,  given 
by  Singh  et  al.  in  Eq.  (6)  and  by  Papadopoulos  et  al.  Eq.  (7)  is 


VJ.max,i 


(4aiV1) 


2/3 


(4a) 


which  in  terms  of  the  dimensionless  background  plasma  quantities  becomes 
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v  .A^^Vil2'3  i/Vsf, 

-max,i  it'  k_jTH  e  \  M^  /  i 


(4b) 


The  minimum  velocity  of  the  stochastic  region,  given  by  Singh  etal. 
Eq.  (8)  and  by  Papadopoulos  et  al.  Eq.  (6)  is 


03 


-imin ,  i  k. 


a: 


1  - 


(5a) 


which  in  terms  of  the  dimensionless  background  plasma  quantities  becomes 


H 


.imin.i  k,rH 


,  !±\/2A\ 
"vlvi/  a± 


(5b) 


In  Table  1,  we  have  calculated  the  minimum  and  maximum  velocities  of 
the  stochastic  region  from  Eqs.  (4b)  and  (5b)  for  the  parameters  of  the 
simulations  of  Singlet  al. ,  shown  in  their  Figures  la,  lb,  2,  and  4, 
and  also  for  the  parameters  of  the  simulations  of  Gaffey  (1981)  and 
Gaffey  and  Mulbrandon  (1982) . 

The  correct  numerical  values  of  v,  .  calculated  from  Eq.  (4b) 

2max, i  n  ' 

2/3 

are  a  factor  of  2  1.6  larger  than  Singh  et  al.  have  in  their  Figures 

la,  lb,  2,  and  4  because  of  the  error  they  made  in  calculating  . 

The  correct  numerical  values  of  v>min  ±  calculated  from  Eq.  (5b) 

are  a  factor  of  2  ~1.4  larger  than  Singh  et  al.  have  in-  their  Figures 

la  and  lb.  Singh  et  al.did  not  calculate  v  .  .  for  their  Figures  2  and  4, 

i.min ,  i  ° 

In  Table  1,  we  have  also  summarized  the  results  of  the  numerical 
simulations  of  Singh  et  al .  and  the  numerical  simulations  of  Gaffey  (1981) 
and  Gaffey  and  Mulbrandon  (1982)  for  singly  charged  ions. 

The  agreement  between  the  theoretically  predicted  values  and  the  nu¬ 
merically  calculated  values  is  quite  good  in  most  cases.  In  general,  it 
is  seen  that  the  numerically  calculated  v 

y  —  IT13K  ^  i 

mi  a 


is  slightly  less  than  the 


theoretically  predicted  value,  and  the  numerically  calculated  v 

■tmax,  i 

is  slightly  larger  than  the  theoretically  predicted  value. 

By  comparing  the  values  of  v.  .  .  and  v  ,  for  a  ~  1/2  v 

i.min,i  -imax,i  i  i 

with  those  for  ^  ~  vi  »  is  seen  that  the  effect  of  increasing 

is  to  decrease  v  .  .  and  to  increase  v  .  ,  thus  broadening  the 

2min , i  amax , i  6 

stochastic  region. 

We  note  that  the  maximum  and  minimum  velocities  given  by  Eqs.  (4a) 
and  (5a)  depend  only  on  the  ion  charge  and  mass,  but  not  on  the  ion  tem¬ 
perature.  This  is  to  be  expected  because  they  were  derived  from  the  equa¬ 
tions  of  motion  of  an  individual  test  particle.  Changing  the  ion  tem¬ 
perature,  changes  the  width  of  the  initial  Maxwellian  distribution,  which 
changes  the  number  of  particles  able  to  participate  in  the  stochastic  ac¬ 
celeration  process,  but  not  the  maximum  and  minimum  velocities  of  the 
stochastic  region. 

The  theoretical  minimum  velocities  for  the  0+  ions  agree  very 
well  with  the  lower  boundaries  of  the  energized  regions  in  Figures  la  and 
lb  of  Singh  et  al.  The  theoretical  maximum  velocities  for  the  0+  ions 
are  about  a  factor  of  two  greater  than  the  largest  velocities  in  Figures 
la  and  lb  of  Singh  et  al.  This  may  be  due  to  the  fact  that  v 

!max,0 

| , 

is  large  and  Singh  et  al.  have  shown  plots  after  only  160  0  cyclotron 
periods.  In  the  simulations  of  Gaffey  (1981)  and  Gaffey  and  Mulbrandon 
(1982)  the  0  ions  with  a  —  10,  continued  to  gain  energy  for  a  very  long  time 
and  eventually  exceeded  the  theoretical  maximum  velocity  given  in  Table  1 
and  exceeded  the  largest  velocities  shown  in  Figures  la  and  lb  of  Singh 

et_al.  for  ct  ,  —  18.5.  The  surface  of  section  plots  of  perpendicular 
0  + 

velocity  versus  gyrophase  for  the  0  ions  contain  a  large  number  of 
islands  of  various  orders  which  slow  down  the  stochastic  diffusion  pro- 
cess.  This  is  another  reason  why  the  0  ions  require  a  very  long  time 


to  reach  the  maximum  velocity. 

For  the  He+  ions,  the  theoretical  minimum  velocity  agrees  very  well 

with  the  lower  boundary  of  the  energized  region  in  Figure  2  of  Singh  et  al. 

The  theoretical  maximum  velocity  is  about  25%  greater  than  the  largest 

velocities  in  Figure  2  of  Singh  et  al.  Again,  this  may  be  because 

v  is  large  and  Singh  et  al.  have  plotted  the  distribution  fun- 

i.max,He  + 

ction  after  only  300  He  cyclotron  periods.  In  the  simulations  of  Gaffey 

(1981)  and  Gaffey  and  Mulbrandon  (1982)  the  He+  ions  with  a  =2.4 

He+ 

required  a  long  time  to  gain  energy  and  eventually  exceeded  slightly  the 

theoretical  maximum  velocity  for  a  =2.4. 

He+ 

For  the  He  ions,  the  theoretical  minimum  and  maximum  velocities 
agree  quite  well  with  the  distribution  function  in  Figure  4  of  Singh 
et  al.  Only  a  few  ions  have  exceeded  the  theoretical  maximum  velocity. 

In  conclusion,  we  emphasize  that  the  values  of  the  minimum  and  maximum 
velocities  of  the  stochastic  region  calculated  from  the  theory  of  Kamey 
(1978)  and  Papadopoulos  et  al.  (1980)  agree  quite  well  with  the  numerical 
simulations  of  Singh  et  al . ,  Gaffey,  and  Gaffey  and  Mulbrandon  in  most 
cases.  The  assertion  made  by  Singh  et  al.  that  an  appreciable  number  of 
helium  ions  are  heated  to  velocities  greatly  exceeding  the  theoretical 
maximum  velocity  is  the  result  of  an  algebraic  error.  Upon  correcting  the 
error  we  find  that  in  all  cases  almost  all  of  the  energized  ions  are  in  the 
stochastic  region  between  the  minimum  and  maximum  velocities  predicted 
by  the  theory  of  stochastic  acceleration. 
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Theoretical  and  numerical 

minimum  and 

maximum  velocities  of 

the  s 

tochastic  region 

in  units 

of  the  ion  thermal  velocity 

ion  species 

ai 

vi 

theoretical 

T,  v  ,  , 

i  JLmin,i 

th  ai 

numerical 

V  V 

imax,i i.min,i 

a .  a. 

i  i 

v 

-max, i 

ai 

Reference 

0+ 

18.5 
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16.1 
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10 
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STATISTICAL  ACCELERATION  OF  ELECTRONS  BY  LOWER- 
HYBRID  TURBULENCE  IN  THE  PRESENCE  OF  A  PARALLEL 
ELECTRIC  FIELD* 

J.D.  Gaffey , Jr. ,  Z.X.  Liu+and  C.S.  Wu  (insti¬ 
tute  for  Physical  Science  and  Technology, 
University  of  Maryland,  College  Park,  Md. 
20742) 

The  statistical  acceleration  of  electrons 
along  ambient  magnetic  field  by  large  ampli¬ 
tude  lower-hybrid  turbulence  in  the  presence  of 
a  parallel  electric  field  is  discussed.  Per¬ 
turbations  driven  by  a  cross-field  current  and 
propagating  nearly  perpendicular  to  the  applied 
magnetic  field  are  considered.  It  is  assumed 
that  the  instability  saturates  rapidly  and  that 
the  fluctuating  electric  field  is  predominantly 
electrostatic.  If  the  turbulence  is  character¬ 
ized  by  a  spectrum  of  small  parallel  wavenum¬ 
bers,  such  that  the  parallel  phase  velocity  of 
the  waves  is  greater  than  the  electron  thermal 
velocity,  then  the  turbulence  can  only  accele¬ 
rate  electrons  moving  with  large  velocities 
along  the  magnetic  field.  The  quasilinear  dif¬ 
fusion  equation  is  solved  using  a  Green’s  fun¬ 
ction  technique,  assuming  a  power  law  spectral 
energy  density.  The  time  evolution  of  an  ini¬ 
tial  Maxwellian  distribution  is  given  and  the 
time  rate  of  change  of  the  average  electron  ve¬ 
locity  and  energy  are  calculated  for  various 
cases.  The  acceleration  process  is  also  dis¬ 
cussed  in  terms  of  a  test  particle  problem. 
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Category  Number  and  Subject  6.1.2,  Magnetosphere 

Theory  CH  Experiment 

Acceleration  of  Auroral  Electrons  by  Lower- 
Hybrid  Turbulence  and  a  Quasistatic  Parallel  Elec¬ 

tric  Field . *  J.D.  GAFFEY,  JR.,  Z.X.  LIU**  and  S.C. 
GUO***,  Univ.  of  Maryland. — The  acceleration  of  elec¬ 
trons  along  the  auroral  magnetic  field  by  large  am¬ 
plitude  lower-hybrid  turbulence  and  a  quasistatic 
parallel  electric  field  is  discussed.  Perturbations 
are  excited  by  a  cross-field  current,  due  to  a  rela¬ 
tive  electron-ion  drift.  The  fluctuating  electric 
field  is  primarily  electrostatic  and  propagates  near¬ 
ly  perpendicularly  to  the  auroral  magnetic  field.  It 
is  assumed  that  the  instability  saturates  rapidly. 

The  time  evolution  of  the  energy  of  an  individual 
test  electron  is  calculated.  The  quasilinear  dif¬ 
fusion  equation ,  with  a  power  law  spectral  energy 
density,  is  solved  using  a  Green’s  function  method. 

The  time  evolution  of  an  initial  Maxwellian  distri¬ 
bution  is  given,  and  the  average  electron  energy  is 
obtained  for  various  cases.  It  is  found  that  tur¬ 
bulent  acceleration  and  electric  field  acceleration 
are  comparable  in  some  cases. 

*Supported  by  ONR  and  NASA. 
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Bistabilities  Due  to  Transmitted  and  Reflected  Ions  and 

Electron  Temperature  Anisotropy  in  the  Quasiperpen- 
dicular  Bow  Shock* 


I 


J.  D.  GAFFEY,  JR,  and  E.  H,  DA  JORNADA  (Institute  for 
^Physical  Science  and  Technology,  University  of  Mary- 
Bland,  College  Park,  MD  20742) 

u.  WINSKE  (Los  Alamos  National  Lab.,  Los  Alamos,  NM 
87544) 


i 


The  theory  of  cross-field  ion  streaming  instabilities 
is  generalized  to  treat  two  ion  species ,  transmitted 

Id  reflected,  and  applied  to  the  parameter  regime  of  a 
gh  Mach  number  quasiperpendicular  shock  wave.  Elec¬ 
tron  temperature  anisotropy  is  also  included  in  the 

»alysis  because  the  electrons  are  perferentially  heat- 
in  the  direction  transverse  to  the  ambient  magnetic 
eld  by  adiabatic  compression  at  the  shock  surface. 
Analytical  and  numerical  results  for  the  linear  growth 

Bte  are  presented.  Three  waves  are  discussed:  a 
arly  perpendicular  propagating  modified  two-stream 
like  mode,  an  obliquely  propagating  ion-ion  streaming 

Kde,  and  an  obliquely  propagating  kinetic  cross-field 
reaming  mode.  The  first  two  modes  are  predominately 
ectrostatic,  whereas  the  third  mode  is  a  whistler 
rave  with  mixed  electromagnetic  and  electrostatic  polar- 
Bation.  The  free  energy  which  drives  the  instabili- 
Res  is  provided  by  the  cross-field  ion  drift  and  the 
electron  temperature  anisotropy,  T^  >  T^jj. 
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Instabilities  at  the  Bow  Shock  I:  Kinetic  Cross-Field 


Streaming  Instability^ 

.  TANAKA,  J.  D.  CAFFKY,  .JR.,  D.  WINSKE,  C.  S.  WIJ, 
Y.  M.  ZHOU,  S.  T.  TSAI  AND  K.  PAPADOPOULOS 
(University  of  Maryland,  College  Park,  MD  20742) 


I 
I 

The  kinetic  cross-field  streaming  Instability  is  the 
fcigh  beta  generalization  of  the  modi fi ed-two-st ream 
liability.  It  propagates  at  oblique  angles  (9  << 

90  )  to  the  magnetic  field  and  is  not  stabilized  by 
electromagnetic  effects  when  the  drift  speed  exceeds 
»ie  Alfven  speed.  It  is  shown  that  the  mode  in  fact 
i.s  an  obliquely  propagating  whistler.  Particle 
simulations  of  the  instability  are  presented  and 

tmpared  with  detailed  linear,  quasi  linear  and 
nlinear  theory.  The  results  are  applied  to  the 
earth's  bow  shock--the  instability  is  expected  to 
Bccur  in  the  foot  region  of  the  quasi-perpendicular 
Jiock — to  explain  electron  and  ion  heating  and  the 
observed  low  frequeny  100  Hz)  electromagnetic 

"ves. 
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instabilities  at  the  How  Shock  II  T:  Cross-Field 
B  Streaming  Instability  with  Tempera to re  Anisotropy* 

i.  T.  TSAI,  M.  TANAKA,  C.  S.  WII  AND  J.  D.  GAKKKY,  JR, 
(University  of  Maryland,  College  Park,  MD  20742) 


1.  Fall  Mooting 
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I  The  fact  that  whistlers  with  oblique  propagation  can 
e  amplified  by  a  cross-field  ion  stream  implies  that 
electron  Landau  damping  Is  unable  to  stabilize  growth 
due  to  the  ion  drift.  On  the  other  hand  It  is  also 

til  known  that  electron  thermal  anisotropy  T  >  T 
n  result  in  v/hlstler  Instability  with  k  «'L0.  11 

It  is  natural  to  question  what  happens  If  both  the  Ion 
^rift  and  electron  thermal  anisotropy  are  Included  in 
Ihe  analysis.  The  present  study  finds  the  thermal 
™nisotropy  enhances  the  peak  growth  rate 
significantly.  Since  In  a  perpendi rula r  shock 

Ilectrons  can  be  heated  preferentially  In  the 
ransver.se  direction,  the  inclusion  of 
temperature  anisotropy  Is  not  only  rea 
■ecessa  ry. 


the  effect 
I  I sf  Ic  hut 


of 
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A  CROSS-FIELD  STREAMING  INSTABILITY 


WITH  ELECTRON  TEMPERATURE  ANISOTROPY  IN  TANDEM  MIRRORS* 

J.D.  Gaffey,  Jr.,  M.  Tanaka,  S.T.  Tsai 

Institute'  for  Physical  Science  and  Technology 
University  of  Maryland 
College  Park,  Maryland  20742 

The  kinetic  cross-field  streaming  instability,  which  is  the  high 
beta  generalization  of  the  mod  i  f  i  ed-tvo--5;t  ream  instability,  has  been 
studied  in  detail,  including  the  effect  of  an  electron  temperature  ani¬ 
sotropy  T,e  >  T  .  The  mode  propagates  at  oblique  angles,  much  less  than 
90°,  to  the  magnetic  field,  and  it  is  not  stabilized  by  electromagnetic 
effects  when  the  drift  velocity  exceeds  the  Alfv/n  velocity*  In  fact,  it 
lias  been  shown  that  this  mode  is  an  obi  i quel y  propngat  i ng  vhi st  1  er  wave. 

The  fact  that  obliquely  propagating  whist ler  waves  can  be  amplified  by 
cross-field  ion  streaming  implies  that  the  electron  Landau  damping  is 
unable  to  stabilize  the  growth  due  to  t  lie  ion  drift.  Because  mirror 
machines  are  of t en  boated  by  eneigetic  neut ral  beams  injected  perpendicular 
to  the  magnetic  field,  it  is  natural  to  investigate  what  happens  when 
both  the  ion  drift  and  the  electron  temperature  anisotropy  are  included  in 
the  analysis.  Since  the  elections  can  be  preferentially  heated  in  the 
transverse  diiertion  in  a  minor  machine,  l  ho  inclusion  of  the  effect  of 
the  t  i-ir.pera  t  ure  anisotropy  is  not  only  i  cal  i  Stic  but  also  necessary.  The 
present  study  has  found  that  the  tempo nturc  anisotropy  significantly  en- 
li.inci'S  1  lie  peak  giowth  rale*.  The  icsults  have  been  applied  to  the  end 
cells  of  the  Tandem  Mill  or  device. 

"fu;  ported  by  ONE  and  NAB  A 


Bffci-t  of  Electron  Tln.-riii.j1  Anisotropy  on  the  Kinetic 
Cross-Field  Instability  in  the  low  "shock* 


J.D.  GAFFEY,  JR.,  MOTOHIKO  TANAKA,  C.S.  WU  AND  e.H. 

DA  JORDANA**  (all  at  Institute  for  Physical  Sci¬ 
ence  and  Technology,  Univ.  of  Md . ,  College  Park, 

Md.  20742) 

The  investigation  of  the  kinetic  cross-field  stream¬ 
ing  instability  (Wu  et  al. ,  (1983),  which  is  mo  t ivated 
by  the  research  of  collisionless  shock  waves  and  ap¬ 
plied  to  the  earth 1 s  bow  shock,  is  discussed  more  fully 
in  the  present  paper.  Since  in  the  quasiperpendicular 
bow  shock  electrons  are  preferentially  heated  in  the 
direction  transverse  to  the  ambient  magnetic  field, 
the  inclusion  of  the  effect  of  the  temperature  aniso¬ 
tropy  is  both  desirable  and  necessary.  The  present 
study  has  found  that  Tej_  >  Te||  significantly  en¬ 
hances  the  peak  growth  rate  of  the  cross-field  stream¬ 
ing  instability  when  the  electron  beta  is  of  order 
unity.  Furthermore,  the  present  analysis  also  im¬ 
proves  the  analytical  and  numerical  solutions  of  the 
dispersion  equation  obtained  by  Wu  et  al.  (1983). 

Wu,  C.S.,  Zhou ,  Y.M. ,  Tsai,  S.T. ,  Guo,  S.C.,  Winske,  D. 
Papadopoulos ,  K.,  1983,  Phvs.  Fluids,  2j6,  (in  press). 
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Category  Number  and  Subject  1.4  —  Scattering,  Emission  and  Absorption  of 
X  Theory  Q]  Experiment 


Radiation 


Spontaneous  Emission  of  Radiation  at  Cyclotron 
Harmonics  by  Electrons  With  A  Loss-Cone  or  Hollow-Beam 

Distribution  Function.*  C.  S.  WU,  H.  P.  FREUND**  AND  J 

D.  GAFFEY,  Jr.,  U.  of  Maryland  —  The  spontaneous 
emission  of  radiation  at  cyclotron  harmonics  by  elec¬ 
trons  with  a  loss— cone  or  hollow  beam  distribution 
function  is  investigated.  The  characteristic  energies 
of  these  electrons  are  assumed  to  be  from  several  keV 
to  several  tens  keV.  In  addition  to  these  energetic 
electrons  there  also  exists  a  component  of  low-energy 
background  electrons  which  has  a  higher  density.  The 
present  investigation  is  motivated  by  the  research  on 
the  auroral  kilometric  radiation  and  cyclotron  emis¬ 
sions  in  mirror  machines.  The  emissivity  is  calcula¬ 
ted  based  on  an  approach  similar  to  that  discussed  in 
a  previous  paper  by  Freund  and  Wul  in  which  the 
general  expression  for  the  emissivity  is. derived  and 
the  relativistic  effect  on  the  cyclotron  resonance  is 
included  in  the  calculation. 
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Effect  of  Electron  Temperature  Anisotrophy  on 
the  Kinetic  Cross-Field  Streaming  Instability  in 

Tandem  Mirrors*  J.  D.  GAFFEY,  JR.,  E.  H.  DA  JORNADA, + 

M.  TANAKA,  D.  WINSKE,  and  C.  S.  WU,  I.P.S.T.,  Univ.  of 
Md.  —  The  theory  of  the  kinetic  cross-field  streaming 
instability}-  which  is  the  high  beta  generalization  of 
the  modified-two-stream  instability,  is  extended  to 
include  the  effects  of  electron  temperature  anisotropy 
since  the  electrons  in  magnetic  mirrors  are  preferenti¬ 
ally  heated  perpendicular  to  the  magnetic  field. 

Because  supplemental  heating  in  mirror  machines  is 
often  supplied  by  neutral  beams  injected  across  the 
magnetic  field,  the  ion  drift  as  well  as  the  electron 
thermal  anisotropy  is  included  in  the  analysis. 
Furthermore,  the  present  study  also  improves  the  ana¬ 
lytical  and  numerical  solutions  obtained  previously  by 
Wu,  et_  al.  ^  The  results  are  applied  to  the  end  cells 
of  the  Tandem  Mirror  device.  The  unstable  mode  is 
identified  as  a  whister  wave  propagating  obliquely  to 
the  magnetic  field.  The  peak  growth  rate  is  signifi¬ 
cantly  enhanced  by  the  electron  temperature  anisotropy 
Te.x;>Te()  when  the  electron  beta  is  sufficiently  high. 
*Supported  by  CNPq,  NASA  and  ONR 
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ELECTRON-CYCLOTRON  MASER  INSTABILITY 
IN  MAGNETIC  MIRROR  MACHINES* 


H.  K.  Wong** 

Laboratory  for  Extraterrestrial  Physics 
NASA-Goddard  Space  Flight  Center 
Greenbelt,  Maryland  20771 
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University  of  Maryland 
College  Park,  Maryland  20742 


The  electron-cyclotron  maser  instability  for  energetic  electrons 
with  a  loss-cone  distribution  function  is  studied  for  magnetic  mirror 
machines.  The  instability  can  occur  at  all  angles  of  propagation  for 
a  wide  range  of  parameters.  The  growth  rate  is  significantly  reduced 
by  the  presence  of  a  population  of  cold  electrons,  and  the  instability 
can  be  suppressed  if  the  density  of  the  cold  electrons  is  sufficiently 
large  and  the  temperature  of  the  energetic  electrons  is  not  too  high. 
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